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1.  INTRODUCTION 
* 

During  the  past  few  years  there  has  been  an  upsurge  of  inter¬ 
est  in  the  effect  of  free  surfaces  on  the  lattice  dynamical  prop¬ 
erties  of  crystals.  There  are  several  reasons  for  this  increased 
activity.  From  the  point  of  view  of  fundamental  knowledge,,  the 
development  of  experimental  techniques  such  as  low  energy  electron 
diffraction  has  made  available)  experimental  information  concern¬ 
ing  the  dynamics  of  surface  atoms.  Those  experimental  results 
have  in  turn  stimulated  theoretical  investigations  of  the  phenomena 
involved.  From  the  point  of  view  of  applied  research,  the  use  of 
surface  elastic  waves  in  devices  such  as  delay  lines  has  led  to 
much  activity  concerning  the  continuum  theory  of  lattice  dynamical 
surface  waves,  particularly  in  piezoelectric  crystals.  Oa  the 
purely  theoretical  side,  the  Investigation  of  surface  effects  on 
lattice  dynanics  has  been  aided  considerably  by  the  availability 
of  high-speed  computers  and  the  development  of  aany-body  theoretic 


techniques .  ; 

In  the  present  review,  the  theory  of  surface  modes  of  vibra¬ 
tion  in  crystal  lattices  will  bo  developed  from  both  the  continuum 

and  discrete  lattice  points  of  view.  The  dispersion  of  surface 

/ 

mode  frequencies  and  the  damping  of  surface 'modes  by  both  anbarmon- 
icity  and  imperfections  will  be  discussed.  Changes  in  normal 
mode  frequencies  produced  by  free  surfacos  will  codify  thermodysiiwaic 
properties  such  as  specific  heat.  These  effects  on  theroodynrnic 
properties  will  be  reviewed. 

Particular’ attention  will  bo  given  to  offsets  associated  with 


-  v 
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'  In  tiSfe  case  of  low  energy  electron  diffraction,  the  Debyo-Waller 
factor  provides  information  about  the  nsean-squaro  displacements 
of  surface  atoms.  The  inelastic  scattering  of  electrons  or  neutral 
atoms  fro'a  crystal  surfaces  can  provide  knowledge  of  the  dispersion 
of  surface  node  frequencies.  Mean-square  velocities  of  surfaco 
atoms  may  be  studied  by  means  of  tho  second  order  Doppler  6hift  in 
the  Mttssbauor  effect.  The  interaction  of  electromagnetic  radiation 
with  crystal  surfaces  nay  ba  manifested  through  either  infrared 
absorption  or  the  Raman  effect.  Both  of  these  phenomena  will  be 
reviewed. 

After  the  present  review  was  started,  the  excellent  treatise  ■ 
on  tho  "Theory  of  Lattice  Dynamics  in  the  Harmonic  Approximation" 
by  Mai-adudin,  Montroll,  Weiss,  and  Ipatova  appeared^  Tho  reader 
is  reforrod  to  tho  chapter  on  surfaco  effects  contained  therein. 
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2. 


SURFACE  MODES  OF  VIBRATION 


A  knowledge  of  the  oorwal  vibrational  modes  of  the  atoms  in 
a  crystal  is  essential  for  the  calculation  of  the  lattice-dynamical 
properties  of  the  crystal.  Thus,  to  study  surface  effects,  it  is 
important  to  know  the  effect  of  a  free  surface  on  the  normal  mode 
frequencies  and  eigenvectors.  Generally  speaking,  the  creation  of 
a  free  surface  tends  to  lower  the  normal  mode  frequencies  and  pro¬ 
duce  a  dess  of  modes  called  surface  modes  in  which  the  displacement 
amplitudes  are  relatively  large  at  the  surface  and  decrease  in 
essentially  exponential  fashion  away  from  the  surface.  These  effects 
may  be  viewed  as  arising  in  the  following  manner.  Starting  from 
a  periodic  crystal,  one  may  create  a  pair  of  free  surfaces  by  sot¬ 
ting  to  zero  the  interactions  coupling  atoms  on  opposite  sides  of 
a  piano  lying  between  two  adjacent  lattice  planes.  This  decreaso 
ia  coupling  constants  produces  a  lowering  in  normal  code  frequencies 
and  an  increase  in  the  mean  square  displacements  of  surface  atoms. 

In  particular,  surface  mode  frequencies  are  typically  split  off  from 
the  bottom  of  a  band  of  bulk  mode  frequencies  having  the  sane  wave 
vector  components  parallel  to  the  surface.  However,  this  does  not 
moan  that  there  may  not  be  other  bulk  modes  with  lower  frequencies 
tbar.  surface  modes  of  the  same  wave  vector. 

Surface  modes  may  be  derived  from  either  the  acoustical  or 

optical  branches.  We  shall  begin  our  discussion  by  considering  the 

continuum  theory  of  surface  vibrational  nodes  derived  from  the 

acoustical  branch.  For  the  case  of  an  isotropic  aediuo,  ve  shall 

be  led  to  tho'oldcst  example  of  surface  vibrational  modes,  namely, 

(2) 

Rayleigh  waves. 


A.  Continuum  Theory  of  Surface  Vibrational  Modes. 


Within  the  framework  of  linear  elasticity  theory,  we  may  write 
Hooke's  law  iu  the  form 

du 

cpq  ’  cpqrs  Tx^ 

where  c  is  an  element  of  tho  stress  tensor,  Cpqra  is  an  element 
of  the  elastic  constant  matrix,  uf  is  tho  r-th  Cartesian  component 
of  displacement  of  the  medium  at  the  point  whose  position  vector 
has  tho  s-th  Cartesian  component  x^.  Tho  summation  convention  for 
repeated  indices  is  assumed  in  Eq.(l).  Wc  shall  take  tho  1  and  2 
Cartesian  components  to  be  parallel  to  tho  free  surface. 

In  order  to  discuss  vibrations,  wo  need  tho  equations  of  motion 
which  can  be  wri-tten  in  the  form 


a2u  ac 
r  _  rs 


,  r  -  3,2,3 


(2) 


whore  p  is  the  density  of  the  medium.  The  definition  of  tho  problem 
is  completed  by  specifying  the  boundary  conditions  at  the  freo 
surface.  If  tho  free  surface  is  specified  by  Xg-  0,  then  tho  bound- 
ary  conditions  correspond  to  vanishing  of  tho  three  components  of 
stress  across  thi3  surface: 


°r3  "  0  at  *3  “  0  i  r  "  1,2,3  (3) 


At  this  point  it  is  helpful  to  specialize  to  a  cubic  crystal 
with  principal  axes  parallel  to  tho  coordinate  axes  and  adopt  the 

t 

Voigt  N  '  notation  for  the  elastic  constants.  Combining  £qs.{l} 
and  (2)  we  obtain  the  equations  of  motion  . 


plus  two  other  equations  obtained  by  cyclically  permuting  (u,v,w) 
and  (x,y,z).  In  Eq.{4),  (u-pUg.u^)  have  been  replaced  by  (u,v,w) 
and  (x^x^Xg)  by  (x,y,z),  respectively.  The  boundary  conditions 
take  the  form 


Bw  Bu 

*S3e  +  32 


'12 


ay/ 


aw  ,  av  _ 
■Si?  +  aS 


♦  On  S'- 0 


<5) 


Surface  wave  solutions  to  Eqs.(4)  and  (5)  are  characterized 

4 

by  an  exponential  decrease  in  the  displacement  amplitudes  in  tho 
positive  z-direction.  We  accordingly  seek  solutions ‘of  the  forn 


(u,v;w)  “  (U,V,W)  exp  |x[-qz+i(tx+ay-ct)3 1 


(6) 


where  c  is  the  phase  velocity,  x  is  tho  magnitude  of  the  two-dimen¬ 
sional  wave  vector  describing  tho  propagation  parallel  to  the  sur¬ 
face,  q  is  s,  dimensionless  attenuation  constant,  and  4,m  are  the 
direction  cosines  cf  the  propagation  direction.  Substitution  of 
Eq.(6)  into  Eq.(4)  yields  the  following  secular  equation  which  must 
bo  satisfied  if  tho  assumed  solution  is  to  bo  non-trivial: 


,2  2  2  2 

+s  "p  — q 

*»(g2+l)  ‘ 

4q(&^l) 


te{g2+l) 

,2  2  2  2. 
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In  Eq.<7),  gj  *"  cu/,c’44i  S2  *’*  ^12^44  3  **ud  p  *•  pc  /c^*  In 


general,  the  solution  of  Eq.(7)  gives  throe  values  of  q'2  for  givon 


values  of  the  velocity  c  and  direction  cosines  -t.n.  Denoting  tho 


three  values  of  q  with  positive  real  parts  by  q  ^ ,  J-1,2,3,  wo  seek 


to  satisfy  the  boundary  conditions;  Eq.(5),  by  the  superposition 
3  , 

(u,v,iw)  -  E^(Sj»nj,Cj)Kjexp  |x[~qjZ+i(fcc+ray-ct)]^ 


(4) 


(8) 


'  where 


§j  -  <4.2+g1m2-p2-qi2)(p2+g1qj2-l)-aS:j2<S2+l) 


V}j  “  •toCg2+l)tqj2(g2+l-g1)+l-p2] 


(9) 


Cj  -  ■tqjCs2+l)Ca2Cg2+l-g1)-t,2+p2+qj2] 


Substitution  of  Eq.{8)  into  Eq.(5>  yields  a  set  of  linear  boiaDgoneous 
equations  in  the  aoplitudes  Kj  wbese  non-trivial  solution  requires 


that 


lfij  l  - 0 


(10) 


whoro 


hi  -  «J-  Vj 


hi  "  mh~  qJnJ 


<11) 


hi  “  *«j+  »V  iall/e\2^i 


The  quantities  q^  which  appear  in  Eq.(10)  are  known  functions  of 
p  frees  tho  solution  of  Eq.(7),  the  secular  equation.  Equation 
(10)  can  therefore  bo  vised  to  dsterraino  tho  surface  wave  velocity 
or  frequency  aud  will  bo  called  tho  frequency  equation. 


i 


-  a-“ 


The  nature  of  the  surface  wavos  in  a  given  situation  depends 


on  the  elastic  constants  of  the  material,-  the  crystallographic 
plane  to  which  the  surface  is  parallel,  and  the  direction  of  propa¬ 
gation  of  the  wave.  We  now  consider  solutions  of  Eqs.(7)  and  (10) 
for  some  special  cases. 

(i).  Isotropic  Case. 

An  isotropic  elastic  medium  is  characterized  by  the  Lame' 
constants  X,u  which  are  related  to  the  usual  elastic  constants  c^, 
Cl2,  and  c44  by 

C11  “  X  +  2(1  *  c12  "  X  •  c44  “  V  *  (12> 

Equations  (12)  imply  a  relationship  among  the  elastic  constants 
which  is  given  by 

t 

C11  -  c12+  2c44  '  *  ‘  .  (13> 

For  an  isotropic  material,  the  surface  wave  characteristics 
are  independent  of  the  direction  of  propagation.  We  therefore  take 
<t«l,  o«*Q  in  Eq.(7)  and  utilize  Eqs.(12)  and  (18)  to  obtain'  the 
solutions 


"  (*-  s') 

(14a) 

"3  '  f1'  i  f 

(14b) 

where  g  »  (X+2j;)/n  ,  We  note  that  p  is  the  ratio  of  the  wave 

validity  c  to  the  velocity  cfc  of  transverse  bulk  waves,  while 
i 

p/g3  is  tho  ratio  of  c  to  the  velocity  c^  of  longitudinal  bulk 
wavos.  In  order  to  satisfy  $he  boundary  conditions,  one  must  super- 
pose. solutions  corresponding  to  and  qg.  The  frequency  equation, 
Eq.(10),  can  be  reduced  to  the  fora 


(15) 


g(p6-8p4+24p2-16)  -  16(p2~l>  **  0 


m 


Solutions  of  £*.(15)  corresponding  to  surface  waves  oust  have 
c<ct'°l  so  ql  ^2  Siven  fay  Sqs.(14)  are  real. 

Physical  values  of  the  reduced  surface  wrve  velocity  p  obtained 
by  solving  Eq.(lft)  range  from  0  96  for  g  -  *»  ,  the  incompressible 
case,  to  0.69  for  g  *»  1.333....,  the  smallest  value'  of  g  consistent 
with  crystal  stability.  The  displacements  associated  with  the 
surface  wave  can  be  expressed  in  the  fora 

u  -  K  [o"*qlz  -(1-5  p2)  **xq22j  eAxix"J-'t;; 


v  -  0 


(16) 


*  -  IK  [l-(p2/g)]^[e"Hqli'*-(l-5p2)"Vwl22]ci>!(5£"ct) 


■ 


where  K\is  a  constant.  We  note  froti*  Eq. (16)  that  the  particle  dis¬ 
placements  executo  ellipses  in  tho  sagittal  planc-i.e.,  tnc  plane 
containing  both  tho  surface  normal  and  tho  direction  of  propagation. 
We  further  note  that  two  attenuation  constants,  qj  and  q2,  are 
required  to  characterize  tho  surface  wave  in  isotropic  material*. 
For  the  incoopr-ssible  case,  q^~  1.00  and  q^  -  3.30.  A  diagram 

illustrating  tho  displacements  for  g  -  3  is  shown  in  Fig.  1, 

(ii)  Anisotropic  Caso,  Cubic  Symmetry 

a.  Surface  waves  on  tho  (001)  surface 

•  For  the  caso  of  tho  (001)  surface  of  a  cubic  crystal, 
one  aust  solve  Eq.(10)  for  tho  reduced  velocity  p  using  the  values 
of  tho  elastic  constant  ratios1  g^  and  g2  appropriate  to  tho  crystal 


-  ■*>  V  v*n- 
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# 
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under  consideration.  For  an  arbitrary  direction  of  propagation  of 
the  surface  wavo,  it  is  difficult  to  obtain  an  equation  analogous 
to  Eq. (15)  for  the  detenainatior.  of  p.  However,  Stonelcy'  was 
able  to  derive  such  equations  for  the  two  directions  of  high  syisaetry 
lu  the  [100]  direction,  he  obtained  the  equation 

a-p2Xgxa-B22-81P2>a  -  C!P4<Ei“P2)  -  0  07)  ’ 

V  * 

while  in  the  [llO]  direction,  he  obtained  the  equation 

<l“p2)(g123~S22“g1p2)2  -  S1P4<g3*p2)  “  0  <1S) 

•shore  g^  -  * (.e-j.+g2*5-2> .  In  other  directions,  it  is  necessary  to 
solvo  &cis. (7)  and  (10)  iinultaneouoly  using  a  computer.’  Tha  dis- 

placeaent  patterns  for  surface  waves  propagating  in  the  [100]  and 

% 

[110]  direction"  are  qualitatively  similar  to  the  isotropic  case*' 
a.®.,  the  displacetseats  are  superpositions  of  twd  attenuated  terns 
and  trace  out  ellipses  lying  in  the  sagittal  plane  with  a  principal 
axis  perpendicular  to  the  surface,  FGr  other  directions  of  propa¬ 
gation,  the  displaecEionta  are  superpositions  of  three  attenuated 
terra?  and  trace  out  ellipses  which  4n  general  are  inclined  at  scae 
accusers  aC&io  to  tho  asgittal  plane. 

It  is  found  upon  solving  -or  the  .reduced  velocity  p  and  the 
attenuation  cons tr  its  q,  tha-  surface  imves  way  e&ist  having  either 
real  or  cosplex  attenuation  cons touts  depending  on  the  elastic  con¬ 
stants  of  the  r.aterial.  Surf  aca  waves  characterized  by  real  atten¬ 
uation  constants  will  be  called  ordinary  Bajloigfc  waves,  while  those 
characterised  fey.  ooaplax  attenuation  constants  will  be  called  gen¬ 
eralised  Rayleigh  waves? ».  -%  - 


*  > 


y 


*  Tha  distinction  between  aatoriala  having  the  two  different 

it 

types  of  surface  waves  can  bo  illustrated  graphically  using  a  plot 

of  tbo  clastic  constant  ratios  g-  vo  g2<  The  case  of  propagation 

in  the  [200]  direction  is  shown  in  Fig.  2.  The  region  of  crystal 

stability  is  that  to  the  right  of  the  lines  g^-gj,  and  g..+2g2"0. 

Points  corresponding  to  fi  nuiabsr  of  well-known  cubic  crystals  are 

shown.  Katerials  to  the  right  of  the  line  labeled  "bound  1"  possess 

ordinary  Rayleigh  waves,  while  those  to  the  left  of  bound  1  possess 

generalized  Rayleigh  waves.  The  line  of  isotropy  g is  shown 

in  Fig.  2  and  is  seen  to  lie  entirely  in  tho  ordinary  Rayleigh  wave 

region  as  expected.  One  aay  acts  that  whether  a  surface  wave  is  An 

ordinary  or  a  generalized  Rayleigh  wave. does  net  siaply  aorrespoad 

to  whether  the  anisstropy  paraaotor  ^c^/vc^-c^o?  xr  less  than  oV 

greater  than  unity,  respectively.  The  lino  labeled  "bound  2"  in 

*  *  _  4 

Pig.-  2  will  6e  discussed  later ,  I 

Since  cubic  crystals  in  general  are  nai-so tropic,  the  velocity 

of  surface  waves  will  depend  on  their  direction  of  propagation. 

The  results  of  solving,  the  frequency  equation,  Rq, (10) ,  for  the 

-  ( c;  v 

surface  wave  velocities  of  several  jaateria.1*  is  shown  in  Fig,'*  3 
whore  the  reduced  velocity  p  is  plotted  as  a  function  of  direction 
of  propagation  on  the  (061)  surface.  Oic  sees  that  the  velocity 
varies  very  slowly  for  the  alkali  halides  KCl  and  HaCt,  hut  shows 
a  no re  pronounced  variations  including  a  aaxi-ua,  for  Al,  dicsond, 
ciuc  rloudc,  and  Oa .  It  aao  hs  noted  that  these  alkali  halides 
exhibit  ordinary  Rayleigh-  wnvss,  whereas  the  other  eftteriaic .enhihlt 
generalized  Rayleigh,  waved,  Another  point  -whit*  say  io  asde  is  * 
that  the  alkali  halide  J4.IV  in  -eahtrast  to  JCCv  and  UaCt,  has .  geaor-  . 


14.  I 


aiized  JUyloigh  waves. 


*hb  throe  attenuation  -constant?  4^j*^2,^3  to  character¬ 

ize  a  Rayleigh  wayo  yary  snoothly  and  rather  uninterestingly  with 


direction  of  propag-ttion  lor  the  ordinary  Rayleigh  wave  aateri&ls 


KOt  tad  SaCt-.  For  &  generalized  Rayleigh  wave  natorlal  such  as 


Cu„  however,  the  attenuation  constants  exhibit  the' rather  interest**  \ 

i  n  y  *  i 


iog  behavior  *>howc  ia  Fig.  4*  Here  q^  and  \v.bo  real  and  imaginary 
parts  o*  q?  and  are  plotted  as  functions  cd  direction  of  propa¬ 
gation.  One  cotes  that  q^  goes  to  zs)o  when  tho  direction  reaches 
the  [llDJ  direction.  This  p»c&nS  that  the  flurfjioc  wave  cEen.^^erates 


to  a  hulk  wave  at  ihe  [HO]  direction.  Early  aur.orical  calculations ' 


indicated  that  q^^  reaches  zero  SsefOre  the  [HO]  direction  is  reached, 

but  recently  nore  precise  calculations^  *  have  shown  that  it  is 

% 

enly  in  the  £  110-3  direction  itself .that  q^  is  zero.  Nevertheless, 
ths  generalized  Rayleigh  waves  penetrate  deeply  into  the  n&torial 
for  propagation  near  the  [110]  direction,  and  for  crystals  of  ordin¬ 
ary  thickness  0. 1-1.0  cs,  surface  waves  will  exist  near  [110]  only 
for  very  short  wave  lengths. 

Froa.tbe  foregoing  one  concludes  -that  the  curves  in  Fig.  3  for 
the  generalized  Rayleigh  wave  materials  correspond  to  bulk  (transverse) 
Wt-vss  when  the  [ZXC]  direction  is  reached.  However,  surface  waves 
:  do  exist  for  these  saaterials  in  the  £110]  direction,  art  tho  vel- 
ocitiaB  are  higher  than  the  Hailing  velocities  in  .the  [no] 
direction  shown  a  Fig.  3,  The  true  ^surface  wave  velocities  arc* 
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solutions  of  Eq.< 18).  For  Ca,  tho  reduced' Velocity  p  of  the  surface 


ways  is  found  to  bo  0.3)  which  stay  be  ecspar?d  t&  the  linitisgro- 
daced  velocity  of  0.56  in  the  £110]  direction,  *  . 
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.  We  now  ask -which  cubic  materials  have  the  surface  wave  velocity 


in  the  (113]  direction  higher  than  the  limiting  velocity  in  that 

direction  for  the  surface  wave  branch.  Tho  answer  is  those  materials 

whose  etaAiic  constant  ratios  fall  to  the  left  of  bound  2  in  Fig.  2. 

Bound  2  itself  corresponds  to  the  equality  of  these  two  velocities 

o 

and  is*  obtained  by  eliminating  p  from  Sq ,  (18)  and  from  the  equation 


!<£;“S2> 


-  0 


(19) 


which  specifies  tX.s  liaHiihr  reduced  velocity  in  tho  [110]  direction. 
The  result  is. 


2<2~g>L+S2)<g1ss*£J;-^/ /  -  ga'ft 


i.  ft£f’ 


(20) 


which  is  the  equation  ic-i  b-5't«.4  > 

b.  Psoudosurface  waves  on  (001)  surfaces. 

Tho  surface  wave  in  the  (110}  direction  on  the  (001) 
plane  of.  Cn  is  the  limiting  case  of  so-called  pseudpsurfaco  waves. 

For  directions  near  the  (1101  direction,  these  waves  involve  dis¬ 
placements  o£  the  type  given  by  Sq,(8),  but  one  of  tho  attenuation 
constants  q^  is  pure  imaginary  rather  than  real  or  complex.  This., 
moans  that  there  is  &  component  of  the  wave  .which  docs  not  attenuate 
in  the  direction  away  from  the  surface;  hence,  tho  n-m?  pseudosuJia.ee 
wave.  The  non-atteuuat^d  component  causes  a  flow  of;  energy  away 
from  the  surface.  This  in  turn,  requires  that  tho  pray®  vector  com¬ 
ponents  parallel  to  the  surface  have  imaginary  parts,  so  that  the 
wave  decays  as  it  propagates  along  the  surface.,  in  many  cases,  tho 


decay  is  so  small  that  the  pccudesvr/amq wevc  i-*  rcadily  obsorvablc. 

-In  Fig.  S e, ,psendosur£Sc»  wave  velocities  calculated  fey  him  and 
Faras!!1  ^ ,  are  chown  Xpr.Sth&  pOi)  suriaec'of:  coj»p'or^  Fdf  assist 
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below  2b°,  the  pseudosurface  waves  are  nearly  pure  bulk  wave3, 

while  at  45°,  the  wave  is  a  pure  surface  wave.  Also  shown  in  Fig. 5 

are  the  curves  for  the  generalized  Rayleigh  surface  wavas  and  for 

(7) 

the  lower  transverse  bulk  wave,  as  well  as  experinontal  points  for 
the  surface  and  psoudosurface  waves.  The  bulk  wave  curve  approaches 
very  closely  to  the  surface  wave  curve  for  angles  greater  than  30° 
and  finally  serges  with  the  surface  wave  curve  at  45°,  As  we  have 
seen,  tho  surface  wave  degenerates  into  a  bulk  wave  at  4Ca. 

Pseudosurface  waves  a?e  typically  exhibited  by  natorials  having 
elastic  constant  ratios  to  the  left  of  bound  2  in  Fig.  2.  A  nunber 
of  specific  oaterials  are  discussed  in  detail  in  tho  recent  review 

article  by:  Parnell . ^  Tha  situation  for  Cu  is,  however,  fairly 

representative  of  that  for  other  cubic  naterials  exhibiting  paeud^- 


surface  wsvaa .  „  • 

c.  Surface  waves  end  pseudo-surface  waves. on  the  (i’Q) 
and  (111)  surfaces 

Rayleigh  waves  exist  on  the  (110)  and  (111)  surfaces 
of  cubic  crystals.  Materials  which  exhibit  ordinary  Bayloigh  waves 
on  tho  (COX)  suylsce  frequently  exhibit  generalised  Bayloigh  vav@s 
on  tho  (210)  surface  one  vice-versc.  On  the  (111)  surface,  the 
surface  ws>ves  ere  typically  of  the  geftBialisjed  Rayleigh  type. 
Psoudosurface -waves-  oa^-exist  on  both  (110)  iknd  (111)  surfaces  < 

Fr.r  further  infonsatien  coccorniny  surface  waves  an.  various  sur--  - 
faces  of  both  cubic  aon-cubic  tiStoriaXs,  the  reader  oay  con¬ 
sult  the  review  by  -Israeli.  - 
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(iii).  Non-Cubic  Materials. 

When  one  investigates  materials  of  lower  symmetry  than 
cubic,  the  calculation-  become  more  formidable  and  a  general,  sys¬ 
tematic  treatment  becomes  desirable.  Such  a  treatment  has  been 

1*0} 

given  for  the  r.on-piozoolectric  case  by  Ingubrigtsen  and  Tonning'  ' 
who  utilized  a  six-dimensional  vector  formalism. 

Specific  calculations  have  been  carried  out  for  a  number 
of  cases.  Deresiewicz  and  Mindlin^  showed  that  surface  waves 
car.  propagate  on  the  AT  cut  of  quartz.  Engan  et.al^11^  have  ex¬ 
perimentally  observed  both  surface  and  pseudosurface  waves  on 
various  surfaces  of  quartz.  Their  data  are  in  reasonable  agree¬ 
ment  with  thooroticnl  calculations. FarncXX^  has 
summarized  results  of  surface  -wave  calculations  for  a  variety  of 
tetragonal,  hex’agonal,  trigonal,  and  orthorhombic  materials.  Some 
of  the  features  revealed  by  these  calculations  will  now  be  dis¬ 
cussed  . 

Surface  waves  propagating  on  the  basal  plans  ((001)  plana) 
of  tetragonal  crystals  show  marked  similarities  tc  surface  waves 
propagating  on  the  (001)  plane  of  a  cubic  crystal.  Some  materials 
such  as  tin  and  zircon  have  ordinary  Royloigh  waves  with  little 
variation  of  surface  wave  velocity  with  direction.  Rutile,  on 
the  other  hand,  exhibits  both  generalized  Rayleigh  waves  and  pseudo¬ 
surf  ate  waves,  and  has  a  narked  variation  of  surface  wave  velocity 
with  direction.  Thu  (110)  plane  also  exhibits  surface  wave  be¬ 
havior  analogous  to  that  found  for  tho  (110)  plane  of  cubic  cvys 
tals.  The  (010>  plane,  however,  has  lower  symmetry  than  the 
corresponding  pluae  of  a  cubic,  crystal,  and  the  surface  wave 
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behavior  wore  nearly  resembles  that  of  the  (110)  plane  than  the 
(001)  plane. 

The  basal  piano  of  hexagonal  crystals  is  isotropic 
for  tl.  propagation  of  surface  waves.  For  some  materials,  such 
a3  CdS,  ZnO,  and  Cd,  the  surface  waves  are  of  the  ordinary  Rayleigh 
type,  while  for  Zn  the  surface  waves  are  of  the  generalized  Rayleigh 
type.  The  basal  plana  of  trigonal  crystals  is  normal  to  a  three¬ 
fold  axis  of  symmetry  and,  as  night  be  expected,  exhibits  surface 
wave  behavior  similar  to  that  of  the  (111)  plane  of  cubic  crystals. 
For  further  results  on  non-cubic  crystals,  the  article  by  Farnoll 
say  oe  consulted. 

(iv).  Piezoelectric  Materials, 
a.  General  discussion. 

% 

Many  piezoelectric  crystals  have  electromechanical 
coupling  constants  which  are  sufficiently  small  so  that  the  piezo¬ 
electricity  can  be  neglected  when  one  is  discussing  surface  elas¬ 
tic'  waves.  For  other  materials,  such  as  lithium  aiobate,  this  is 
no  longer  the  case,  and  the  surface  wave  problem  must  be  reformu¬ 
lated  to  include  the  piezoelectric  coupling.  This  has  been  done 
by  Tseng  and  White .and  by  Xraut.  ^  The  three  Cartesian 
components  of  displacement  must  be  ‘augmented  as  variables  by  the 
scalar  potential  <p  associated,  with  ths  electric  field  produced  by 
the  strains.  The  stress  ia  related  to  the  strain  and  to  the 
gradient  of  «?  by 
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wbara  the  are  the  components  of  the  piezoelectric  coupling 
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tonsor.  The  equations  of  motion  are  now  obtained  by  substitution 
of  £q.  {  21  )  into  Eq.  £  2  )  giving  the  result 


a  /  %  t  at?  \ 

7?  "  3^  \Grspq  ST"  +  °prs  *x“  )  . 


The  constitutive  relation  involving  the  electric  displacement 
D  and  the  dielectric  tensor  e  „  ■ 
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together  with  the  Maxwell  equation 


yields  the  equation 
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(251 


which  Is  to  be  solved  simultaneously  with  Eq.  (22). 

The  boundary  conditions  for  a  surface  cefined  by  the  plane 
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3Cj“0  have  the  .torsi 
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corresponding  to  the  vanishing  of  the  three  components  cf  stress 


at  Xj^O,  and 
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corresponding  to  the  continuity  of  tho  normal  component  of  g  at 
Xg  “  0.  in  Hq.{  27‘),  $  is  the  electric  potential  outside  the 
crystal.  It  stunt  satisfy  Laplace's  equation  and  approach  aero 
as  x3 

For  a  general  direction  of  propagation,  one  has  4x4  deter¬ 
minants  in  the  secular  equation  and  the  frequency  equation  rather 
than  3x3  determinants  as  in  the  non-piezoelectric  case.  Surface 
wave  solutions  have  been  calculated®^  for  a  cumber  of  piezo¬ 
electric  materials.  In  an  approximate  sense  ono  can  express  the 
surface  wave  velocity  in  terns  of  elastic  constants  "stiffened’’ 
by  the  piezoelectric  coupling.  Another  consequence  of  the  piezo- 

i 

electric  coupling  is  the  existence  of  a  new  type  of  surface  wave, 
tho  Blouste in-Gulyaev  wavc.^2&’2®^  . 

Specific  calculations  for  "the  interesting  material-  ■ 

(  17} 

lithiua  niobata  have  been  carried  cut  by  Cambell  and  cones.  ' 

Their  results  show  that  omission  of  the  piezoelectric  coupling 

leads  to  surface  wave  velocities  which  are  on  the  order  of  15% 

P3} 

too  low.  Tseng  and  White v_  7  have  calculated  surface  wave  veloc¬ 
ities  on  the  hasal  plane  of  CdS,  CdSe,  ZnO,  and  lead-titanate 
zirconate.  The  effect  of  various  electronagnotie  boundary  con¬ 
ditions  on  surface  wave  propagation  in  piezoelectrics  has  been 
investigated  by  Ingebrigisen,  AO' 

-  b.  Bleustein-Gulyaev  surface  waves.  . 

In  the  case  of  piezoelectric  crystals  a  surface 
wave  is  possible  which  has  no  counterpart  in  ordinary,  non- 
piezoolectric  crystals.  This  wave  was  independently  discovered 
by  Bloustoin'15^  and  Gulyaev,  It  consists  of  a  transverse 

wave  having  its  mechanical  displace neat  parallel  to  the  surface. 
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'An  ordinary  earfoco  wave,  on  the  othor  hand,  has  Its  noebanlcal 

<t 


displacement  nor.-parallel  to  tho  surface.  Considor  a  siaterisl 
such  as  CdS  belonging  to  class  Cgv  and  having  ita  six-fold  axis 
parallel  >£o  tho  direction.  If  wo  considor  a  surface  defined 
by  =  0  and  restrict  ourselves  to  displacoaents  only  in  -no 
3-direction  which  are  functions  only  of  and  then  Eqs,  <  22  ) 
and  (  25  )  reduce  to 

x2.. 
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v2  u3  -  €X1  91  <? 


(28a> 

(28b) 


where'  V2  is  tfcd  Laplaciaa  in  tho  variables  and  the  Voigt 

notation  has  been  used  fer  the  elastic  and  piezoelectric  constants. 
The  non-'Vivial  boundary  conditions  at  Jtg* 0  can  be  written  as 


Taking 


ties  that 
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(29b) 

~  (o15 /€UJ  U3  mld  c44  “  ?44*  C15/,'ll  * 

COS  Ooc^-HJtJcXpC-QXg),  X2>  0 

(30a) 

cos  (kx^*ait)exp(-kx2) ,  Xg>  0  • 

(30b) 

cos  (kx^-ojtJexpCteg)  ,  x2<  0  . 

•  (30c) 

tting  into  Zqs .  <  28  ),  wo  find  for  a  noc-trivlal  solu- 

ptii2  -  c44(fcS-(i2)  . 

(31)  . 
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Substituting  Eqs.{  30  >  into  Eqs.(  29  )  wo  gat  the  relation 
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Combining  Eqs.(31)  *  and  (32),  w®  obtain  for  the  speed  of  the 


Bloustoin-Gulyacv  wave 
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It  can  bo  seen  that  the  speed  is  somewhat  loss  than  the  transverse 


bulk  wave  speed  calculated  with  tho  piezcolectrically  stiffened - 


elastic  constant  c. 


Bleustein' 


has  also  treated  tho  case  lr.  which 


tho  ourface  is  covered  with  a  grounded  this  perfectly  conducting 
layer.  Gulyaev has  considered,  the  situation  where  free  charge 


earriors  are  present  in  tho  crystal.  Tseng1.  1  has  shown  that 
surface  waves  of  tho  Bleustein-Gulyaev  typo  can  exist  on  certain 


surfaces  of  certain  classes  of  cubic  and  orthorhombic  crystals. 


Specific  results  are  given  for  Bi12  Ge020  (cubic)-  and  Ba2  Na  Nbs  015 


(orthorhombic) . 


B.  Lattice  Dynamical  Theory  of  Surface  Vibrational  Modes. 


(  i).  Static  Displacements  Near  a  Surface. 


In  the  lattice  dynamical  theory,  one  assumes  that  the 


material  undor  consideration  is  made  up  of  discrete  atoms  or  ions 


which  vibrate  about  certain  equilibrium  positions.  For  an  infinite 
crystal  those  equilibrium  positions  are  associated  with  the  sites  of  a 
crystal  lattice  periodic,  in  three  dimensions.  The  crystal  can  be 
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considered  as  made  up  o t  unit  cells,  oach  cell  containing  a  certain 
numbor  of  atoms.  Tho  position  vector  r(ix)  of  the  x-th  :\ton  in 
the  -t-th  unit  cell  may  be  expressed  in  tho  form  r(<-x)  -  + 

u(ix)  where  (4x)  is  the  position  vector  of  the  equilibrium 

site  and  u(ix)  is  the-  displacement  vector  from  equilibrium.  2a 
terms  of  the  primitive  translation  vectors  of  the  lattice  a^,  • 

&2‘  %3‘  oae  can  WI"*t0 

S(o)(ix)  «  ^  Sl  +  a2  +  l3  a*  4-  R(x)  (34) 

where  l2>  arc  integers  and  R(k)  is  the  position  vector  of 
the  x-th  atom  relative  to  tho  origin  of  the  unit  cell. 

He  shall  consider  semi-infinite  crystals  with  a  free 
surface  parallel  to  some  lattice  plane.  The  crystal  is  periodic 
in  the  two  dimensions  parallel  to  the  surface  but  not  in  the  di- 

*  i- 

mension  which  is  normal  to  the  surface.  Two  basis  vectors  and 
a,  may  be  chosen  parallel  to  the  surface,  while  the  tpird  basis 
vector  Bq  is  chosen  not  parallel  to  the’  surface.  These  basis  vec¬ 
tors  may  or  may  not  bo  the  primitive  translation  vectors  of  the 
infinite  lattice.  The  equilibrium  position  vector  R(-tx)  for  an 
atom  in  the  semi-infinite  lattice  can  be  expressed  as 

V 

5(4x>  -  S(o)(4x)  v  S(4x)  05) 

whore  the  vector  5(ix.)  tahes  into  account  tho  ’'relaxation"  cf 
the  equilibrium  atonic  sites  in  the  semi-infinite  crystal  away 
i rdm  the  corresponding  sites^  of  the  infinite  crystal. 

f  '  * 

The  static  displacements  S(-tx)  are  specified  by  tbs 
equilibrium  conditions 
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where  i  is  the  potential  energy  of  the  crystal  as  a  function  of 
the  atonic  positions  and  c  is  the  number  of  atoms  in  a  unit  cell. 
If  the  static  displacements  are  small,  they  stay  be  determined  by 
expanding  the  forces  appearing  on  the  left  in  Eq.  (36)  in  power 
series  about  the  configuration  the  atoms  would  hove  at  equilibrium 
in  the  infinite  lattice:  ' 

Fc Wx)  -  ? <°>  (U)  -  i  ;S, p  4  <  J>  ax U'k"}  S^Cl'K’)  +.  ; . .  <3?> 


where 


Fa  'W  “  Tr^CZxT  S(tx)  -  0  , 


'  ’•4"J  "S'! 


(ixJU'K'J  - 
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S(4x),sa*K*)  -  0  ». 


and  higher  terms  in  the  series  expansion  are  neglected.  One  should 
note  that  the  forces  ®tx)  are  not  all  zero  because  4  is  the 
potential  energy  for  the  semi-infinite  lattice;  in  fact,  the 
^  are  non-zero  for  those  atoms  whose  distance  froa  the 
surface  is  less  than  the  range  of  the  interatomic  fo.  jes.  The 
•aTe  ti3°  harmonic  coupling  constants  when  the  atoms 
are  in  the  configuration  appropriate  to  the  infinite  lattice.  For 
some  assumed  set  of  interatomic.- interactions,  the  static  displace¬ 
ments  can  be  calculated  from  Eqs,  <3®  and  {  37)  by  solving  a  set 
of  inhomogeneous  algebraic  equations. 

Physically,  the  semi- infinite  lattice  may  be  regarded  as 
arising  from  the  Infinite  lattice  by  tbs  removal  of  all  atoms 
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beyond  the  bounding  plane  dofining  the  free  eurface.  Atoms  near 
the  surface  are  then  acted  upon  by  unbalanced  forces, 
duo  to  the  removal  of  the  Interactions  with  atoms  across  the 
bounding  plane  and  therefore  relax  to  new  equilibrium  positions. 
Symmetry  frequently  requires  that  the  static  displacements  bo 
normal  to  the  surface,  but  this  is  not  always  the  case.*'20^ 

A  general  discussion  of  the  relaxation  of  atoms  near  a 
free  surface  has  been  given  by  Feuchtwang. A  number  of  cal¬ 
culations  have  been  reported  based  on  numerical  solutions  to 
Eq.  (36)  using  models  involving  Leonard- Jones  or  Morse  -potential 
interactions  between  pairs  of  atoms.  For  example,  various  sur¬ 
faces  ct  face-centered  cubic  crystals  with  Lennard-Jones  inter¬ 
actions  have  been  studied  by  Shuttleworth, Alder  ot  al,^22^ 

Schmidt  and  Jura, ^23'  Benson  and  Claxton, ^ 2<1 ^  and  Allen  and 
(25) 

Do  Wette.  -  -  The  results  typically  show  an  expansion  outward 
of  the  crystal  at  the  surface.  The  increase  in  tho  interlayer* 
spacing  is  on  the  order  of  2-3%  at  the  surface  and  decreases 

it- 

rapidly  toward  tho  interior.  Allen  and  de  Wette ^25)  have  shown 
that  for  an  interatomic  potential  varying  as  l/rp,  p  >  3,  then 
the  increase  in  interlayer  spacing  varies  at  large  distances  d 
from  tne  surface  as  1/d1^3.  Thus,  for  a  6-12  Lennard-Jones 

rx 

potential,  the  spacing  varies  as  1/d  ,  a  result  previously  given 
by  Alder  et  al.^22^ 

Tho  (001)  surface  of  a  staple  cubic  lattice  has  been 
treated  by  Yun  and  Benson using  a  Lennard-Jones  potential. 
They  find  only  a  0,5%  increase  in  interlayer  spacing  at  the  sur- 

(27) 

£acs«  Vail  9  on  tho  other  handy  used  a  liorse- potential  and 


found  a  2%  incre&ce.  More  extensive  use  of  the  Mo'.ase  potentia 
has  been  reported  by  Jackson'  *  who  .s^.dieQ  the  (100)  (IIP’, .. 
and  (111)  surfaces  s.  .  variety  of  ..•.je-contered  and  body-centeree 
cubic  metals.  lie  f'Jund  some  surp.l '*n^ly  large  relaxations  of  the 
surface  layex  s.g.,  20%  for  the  (2’0)  surfaces  of  Ca  nnd  Sr, 
and  80  -  95%  for  the  (111)  surfaces  of  alkali  csetsls.  Whether 
such  largo  deviations  iron  the  bulk  interplaear  spacings  are 
physically  reasonable  remains  to  ba  seen, 

considerable  wbrk  has  been  done  on  ionic  crystals  of 

{pft  ) 

tha  alkali-halide  typo  by  Benson  and  coworkers.  For  tho 

(100)  surfaces,  they  find  typically  an  alternating  contraction 
and  expansion  of  the  surface  layers.  The  surface  layer  con¬ 
tractions  range  froa  about  seven  percent  for  Hal  to  one-half 
*  * 
percent  for  HbF. 

Benson  and  coworkers  allowed  only  the  £ irot  five  'surface 
layers  to  relax,  but  pored t ted  separate  relaxation  of  the  cations  and 
anions  within  those  layers.  Tong  and  Maradudin^30'  have  followed  a 
'different  tuck  by  considering  a  NaCl  crystal  of  15  layers  and  allowing 

all  layers  to  relax  with  respect  to  both  intralayer  and  interlayer 
spacings.  They  did  cot,  however,  allow  separate  relaxation  of 
cations  and  anions  and  required  tho  itttralayev  relaxation  to  be 
the  .ease  in  all  layers.  They  found  aiy  intralayer  contraction  of 
about  Irrelative  to  the  Infinite-crystal  nesrent-seigbbor  spacing  ' 
and  an  inf^rlsyer  expansion  of  about  1%  .relative  to  the  nsw  infra- 
layer  sonrast-neifhbor  ejscieg.  The  fnterlay«r  expansion  is 
ss&Iier  at  the  surf  sen  than  iu  tbs  .interior. 
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*  Relatively  few  analytical  results  for  the  static  dis¬ 
placements  near  a  surface  have  been  presented.  Gaais  and  Wallis 


(31) 


have  solved  Eqs.  (36)  and  (37)  for  a  semi-infinite  linear  chain 
with  nearest  and  next-nearest  neighbor  interactions.  A  similar 
calculation  has  been  carried  out  for  a  body-centered  cubic  lattice 
with  a  (001)  surface  by  Clark  et  ai^2)  us«ag  kennard-Jones  inter¬ 
actions  between  nearest  and  next-nearest  neighbors.  In  both 
cases  the  static! displacements  were  found  to  vary  as  exp(-q<j) 

for  a  surface  at  l3  -  0  whore  l 3  is  an  integer  a  0  labelling  the 
lattice  planes,  q  has  the  form  q  -  qQ  +  iuD,  qQ  is  a  real 
positive  lumber  on  the  order  of  unity,  and  8  is  0  or  1  depending  - 
on  the  force  constants.  This  exponential  drop-off  of  the  static 

.  k 

displacements  with  distance  from  the  surface  is  to  be  contrasted 
with  the 'inverse  power  law  dependence  obtained  by  Alder  et  alv“ 
and  by  Alio.;  and  Do  U6ttc^S^  for  Lehnard-Jones  interactions 
between  all  atom  pairs  rather  than  just  nearest  and  next-nearest 
neighbors.  * 

(li).  Equations  of  Motion  for  a  Crystal  with  a  Surface. 

The  vibrations  of  the  atoms  in  the  crystal  nay  be  ana¬ 
lyzed  by  expending  the  potential  energy  in  power  series  in  .the 
components  of  the  displacements  from  the  equilibrium  sites  of  the 
semi-inf ini  to  lattice,  • 
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ar.d  (u(4k)  3  stands  for  the  entire  set  of  displaceaonts  «(■£«)  • 

The  coefficients  ia^(tK,VK’  ,  etc,,  are  the  harmonic  cou¬ 

pling  constant*,  the  cubic  anbaraonic  coupling  constants,  etc. 

The  coupling  constants  must  satisfy  certain  invariance 
conditions.  The  periodicity  of  the  lattice  parallel  to  the  sur¬ 
face  (specified  by  -'0)  requires  that  the  coupling  coefficients 
be  functions  only  of  the  differences 

.-Cj- 1  Ij,  etc;  however,  they  are  in  general  functions  of  t3,  -14, 

,  •  • 

4g,  ...  separately.  Infinitesimal  translational  invariance  re¬ 
quires  that 

»  • 

T*  y>  V 

4V  *a9C4,l»4,x'>  "°»  A*  A”  “°'  etc“*  (42> 

a 

while  infinitesimal  rotational  invariance  requires  that 

A‘{*a&<4*‘4’*,)[V4,0'“  RYCt,X,)] 

■  W(t*’4’x’>[V'tR>  '  V*  **>]}  "  °  (43) 

with  corresponding  equations  for  the  anharsonic  coupling  coeffi¬ 
cients.  The  constraints  imposed  on  the  coupling  coaifdcients  by 
rotational  invariance  are  particularly  important  fer  surface 

problems.  ^33 
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The  normal  nodes  of  vibration  for  the  semi-infinite 
crystal  nay  bo  found  by  solving  the  harmonic  equations  of  motion 
for  the  displacements  u(4x).  Using  £q„  (40),  these  equations 


can  be  written  as 


a 


c4 

cu^T VtCf 


<  44) 


-  *aSU*>1  * >*e(l ‘> 

4’x'P  a*  9 

I 

where  we  have  used  the  syoma try  property  (ix.,  V  >,*)“ 

k*,4x)»  The  solution  of  £q.(44)  is  complicated  by  two  factors. 

_  • .  . 

First,  there  is  no  periodicity  of  the  atomic  sites  i.*  ths  direc¬ 
tion  normal  "to  the  surface;  second,  the  coupling  coefficients 
i=  AtK&x1)  are  not  necessarily  the  same  as  tho  corresponding 

Rp  ,  /  % 

coefficients  ‘for  the  infinite  lattice,  •  - 

A  general  discussion  of  the  coupling  constants  in  a  semi- 


infinite  crystal  has  been  given  by  Fcuchtvang. 


One  may  expand 
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tho  coupling  coefficients  in  power  series  in  the  static  displace¬ 
ments  S(4x): 

4a3Ux*4‘y'‘>‘5G?(0><i5i'4,,l,)'!'  2  „  \g°U4K,-t,x',t«x«)Sy(4"H-) 

+  ...  "  (45) 

It  is  clear  from  Eq.fes)  that  the  changes  in  the  coupling  constants 
arise  from  anharmonic  effects  and  that  they  will  be  large  aaar  the 
surface  where  the  S(4,h)  sire  large.  The  coupling  coefficients  with 
the  superscript  zero,  it  must  be  Emphasised,  are  still  cot  osci-ssarily 
equal  to  the  corresponding  coefficients  of  the  infinite  lattice, 
.because  the  atoms  near  the  surface  ca^  be  polarized  even  -*hca  thsy 
•occupy  tha  sites  th»ty  would  occupy  in  infinite  '„■= 
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Changes  in  surface  coupling  constants  have  boon  calculated  by 
Clark  et  al^32^  for  the  case  of  a  (001)  surfaco  of  iron  using  a 
model  with  Leonard- Jones  interactions  between  nearest  and  next- 
acarest  neighbors.  They  found  relatively  large  decreases  (20-30%) 
in  certain  of  the  coupling  constants  at  the  surface  compared  to  the 
bulk  values  even  though  the  static  displacements  at  the  surface  arc 
only  about  three  percent  of  the  bulk  lattice  spacing. 

Turning  now  to  the  solution  of  the  equations  of  notion,  we 
exploit  the  translational  periodicity  parallel  to  the  surface  by 
seeking  solutions  to  2q.(44)  of  the  form 

uaUR)  -  M^va(43x,kp)  ei(kl¥¥fwt)  (46  ) 


whore  u)  is  the  circular  frequency  and  fc_  is  a  dimensionless  two- 
dimensional  wave  vector  parallel  to  the  surface  with  components 

kj  and  kg*  Substitution  of  Eq.(46>  into  Eq.(44)  yiolds  a  set -of 

*  *  2 
linear  difference  equations  in  the  amplitudes  v' (4x,uj  )  , 


m  va(43n»V  D«S(<3X»*3K’  'Sp5*^*' ’  V 


<  47) 


where 

(48  ) 


are  the  elements  of  the  reduced  dynamical  matrix  of  the  semi-infinite 
crystal. 

The  solution  of  the  equations  of  notion  specified  by  Eq.(44) 
is  more  complicated  than  for  an  infinite  crystal  for  two  reasons. 
First,  a  given  atom  does  not  interact  with  atoms  beyond  the  bound¬ 
ing  plena  of  the  crystal  bo  certain  terms  are  missing  from  Eq.(44). 
Second,  the  coupling  coefficients  that  are  actually  present  may 
vary  with  the  distance  of  the  coupled  atoms  from  the  surface. 
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.  Th*e  noraal  mode  frequencies  of  the  crystal  aro  specific!  by 

ft 

setting  the  determinant  cf  the.  coefficients  of  the  amplitudes  in 
2q.(47>  equal  to  zero: 


|Dc0<<3K,<3*'»£p>  _u|2  6kk‘  6aj3 


-  o 


(  49  ) 


For  a  semi-infinite  crystal,  the  determinant  in  Eq.(4S0  is  infinite 
in  size.  A  finite  determinant  can  be  obtained,  however,  by  intro¬ 
ducing  a  second  surface  parallel  to  the  first  and  specified  by 
■Lj-Lg.  One  is  now  dealing  with  a  plate  having  L3  atonic  layers. 

The  size  of  the  determinant  in  Eq.(49)  is  then  3sJLjx3sL3  where  s 
is  the  number  of  atoms  in  a  unit  cell.  Results  that  seem  reason¬ 
ably  representative  of  real  crystals  aro  obtained  if  is  as  small 
as  20  or  30.  For  simple  crystals  one  than  must  work  with  60x60  or 
90x90  determinants,  which  can  easily  bo  handled  using  modern  high 
speed  computers. 

Ac  a.1  terns, iiva  approach,  particularly  useful  situ,  short-range 

v 

forces ;  is  to  regard  each  atosa  as  having  the  same  equations  of  jKoticsv, 

They  must  be  augmentod  by  boundary  conditions,  which  are  constraints 
that  restore  the  effect  of  the  surface  on  the  equations  of  motion. 

The  procedure  is  completely  analogous  to  that  employed  with  the 
continuum  theory  of  surface  elastic  waves.  The  determinants  which 

V  • 

occur  in  the  analogues  of  the  secular  equation  and  the  frequency 
equation  are  gonerally  on  the  order  of  3sx3s  in  size,  which  aro 
relatively  small  and  are  easily  handled  on  a  computer. 

Some  general  aspects  of  solutions  to  the  equations  of  motion, 

Eq . (44) ,  have  been  discussed  by  Feueh twang. For  a  semi-infinite 
crystal,  such  solutions  may  be  ’’bulk”  solutions  whore  tbo  displace— 

,  moots  ua<4n)  do  not  approach  zero  as  4-t  «  and  "surface”  solutions 

w 

where  they  approach  zero  as  •tj’*  «.  For  a  finite  slab,  the  surface 

and  hulk  solutions  differ  according  to  whether  the  displacements 
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do  or  do  not  approach  zero,  respectively,  at  the  center  of  the  slab 


as  tho  thickness  becomes  very  large.  In  the  following  we  shall 


consider  tho  surface  solutions  for  a  variety  of  special  cases. 


(iii).  Linear  Monatomic  Chain 


A  specific  illustration  of  soae  of  the  points  mentioned 


above  is  provided  by  the  semi- infinite  monatomic  chain  with  nearest- 


neighbor  interactions.  The  equations  of  motion  for  the  displacement 


u  of  atom  n  are 


M  «*  a’Cug-Uj^) 

M  u 2  -  aiU3-u2)+  a 


(  50a  ) 


(50b  ) 


M  un  -  tt(u.+u1-2u) ,  a  a  3 


(50c  ) 


whero  a  and  a*  are  the  coupling  constants  for  bulk  atoms  ancf  f or 


surface  atoms,  respectively,  H  is  the  atvaic  mass,  the' double  dot 


indida'es  a  second  fciae  d<<ri  votive,  and  the  surface  atom  is  speci¬ 


fied  by  o  “  0. 


For  an  infinite  chain,  every  atom  fcas  an  equation  of  ssot-ioa 


of  the  form  of  Eq.(50e).  The  normal  nodes  of  vibration  can  be 


written  as  .travelling  waves 


ua  -  U  expCiGpa-cjt)] 


<51  ) 


where  9  is  a  dimensionless  wave  number  and  tn  is  the  frequency 
specified  by 


ta2  -  (4a/U) sin2 (<p/2)  . 


(52  > 


The  wave  vector  t?  ranges  from  -it  to  +tt,  so  that  the  squared  -frer 


quencies  range  from  0  to  4o/£; 
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For  tho  serai-infinite  chain,  we  can  obtain  surface  solutions 
under  certain  conditions.  Let  us  write  the  displ&cenents  in  tho 
fora 

U1  ”  Ujcxpttwt)  (S3a  ) 

uQ  -  U  exp(-qn+imt) ,  n  >  1  .  (S3b  > 

Substitution  of  £q,(53b)  into  Eq.(SOc)  yields  tho  result 

lies2  -  2a (1-cosh  q)  .  (54  ) 

An  equation  specifying  the  attenuation  constant  q  is  obtained  by 
substituting  Eqs.(5ai)  and  (53s)  into  Eqs.(5<h)  and  <50>),  using 
Eq.(54),  and  sotting  the  determinant  of  the  coefficients  ci  !>,,  ■ 

U  equal  to  zero.  Tho  result  can  bo  written  as 

"  -c-(j£2+c)^  (55  ) 

■sfhere  c  «•  (a'-a)/a.  For  a  surface  mode,  tho  real  part  of  q  must 
be  positive.  Such  &  solution  to  Eq.(55)  exists  for  c  >  1/3  and 
has  the  fora  q  «  qQ-s-ln  where  qQ  is  real  and  positive.  The  fre¬ 
quency  is  now  specified  by 

to  -  [(2a/i5)  (1+cosh  Qd)3^  (56) 

j 

and  is  larger  than  the  maxir.ua  buXH  aode  frequency  (4o/a)-*.  Thus, 
for  c  >  1/3,  a  surface  aodo  frequency  rises  above  the  allowed  band 
of  frequencies  of  the  inflate  chain.  It  should  be  noted  that  no 
surface  mods  exists  when  tho  surface  coupling  constant  a'  is  un¬ 
changed  from  the  bull:  value  a. 

If  next-nearest  neighbor  interactions  are  included,  the  sent- 

ir  v 

infinite  aonatoaic  linear  chain  still  docs  not  possess  a  «?urf aca  node 
'  if  tho  surface  coupling  constants  arc  unchanged  froa  the  fciilfc  • 

,<31> 

I  «• 


values 


One  or  two  surface  inodes  can  bo  produced,  however,  by  suitable 

(31 1 

changes  in  the  surface  coupling  constants. 

Surfaco  modes  can  also  bo  produced  by  changing  the  mass  of 
the  atom  at  the  free  end  of  the  chain.  This  situation  corresponds 
to  an  absorbed  atom.  For  nearest- neighbor  interactions  only,  one 
finds^3^  that  a  surface  mode  exists  provided 


where  M'  in  the  mass  of  the  and  atom. 

If  the  surface  coupling  constant  cf  has  the  bulk  value  a,  a 
surface  mode  exists  if  <  11/2.  This  result  has  also  been  obtained 
by  Malsuda^3^  and  by  Asahi  and  Hori^3®  using  the  transfer-matrix. 
The  closely  related  seattoring-“3S5atrix  method  has  been  used  by 
Hori'and  Asahi*?7'  to  treat  the  ca&e  of  tvo  isotopic  impurity 
stems  at  the  free  end  of  a  r,on*tomic  chain  with  nearest-neighbor 
interactions.  If  the  two  impurity- a  fconis  have  the  same  mass  H', 
surface  modes  exist  if  Ur  <  J£(2-,/2)V4.  This  problem  has  also 
been  treated  by  Ifesri  and  Anawsd'  *  who  discussed  cne-dlacnsional 
raodols  for  argon  condensed  on  krypton  and  for  esrbon  monoxide 
adsorbed  cc  platinum.  The  results  for  the  latter  case  agree  well 
with  those  obtained  using  ft  three-diKensiohsl  acdsl, 

That  perturbations  in  masses  cr  coupling  constants  are  required 
to  produce  a.  surface  mode  in. the  scat-infinite  manstcaic  linear 
chain. is  understandable  ofl  the  basis  of  Rayleigh’s  theorem.^1/ 
Without  such  perturbations,  the  creation  of  the  surfaco  correspond? 
to  setting  ono  o?  caaplisg  constants  to  soro.  and  this  by 
Rayleigh's  theorem  must  depress  the  aovaal  isede  frequencies.  For 
»  inc-diaiaasics&l.syatem,  a  localised  eoia  such. an  a-  surface  mod? 
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cannot  have  its  frequency  in  an  allowed  band  of  frequencies  of  the 

infinite  lattice.  Since  the  allowed  band  extends  down  to  zero  Ire- 
% 

quency  and  since  the  surface  mode  cannot  have  a  negative  frequency, 
there  is  no  gap  into  which  a  depressed  frequency  can  aovo  and 
thereby  produce  a  surface  node , 

(It}.  Diatomic  Linear  Chain 

Tho  serai-infinite  diatonic  linear  chain  offers  an  inter¬ 
esting  contrast  to  the  serai-infinite  raonatoraic  linear  chain  in  that 
a  surface  node  is  possible without  changing  the  surface  force 
constants,  provided  the  atom  at  tho  free  end  is  the  lighter  of  tho 
two  atoms  in  the  chain.  If  ono  assumes  nearest-neighbor  interactions, 
the  equations  of  notion  have  the  fora 


a  ■  afu^-u^) 

(  58a) 

U  hi  *  C(u2j+l'>u2j-r2u2j>  »  1 

(  58b) 

a  u2j_a  -  a<«2j+u2j-2~2u2j-r  »  J  4  2 

(  58c) 

where  tho  end.  atoa  is  taken  to  have  the  scalier  of  the  two  masses 
a  <  U.  A  surface  solution  to  £qs.(58)  exists and  can  be 
writtoa  as 

u2j-l  “  exp{iot)  (  59a) 

u2j  ”  U(-X)^(a/H)^  exp(iint)  (  59b) 

where  the  frequency  o  is  specified  by  1 

m  -  (  60  ) 

The  atoaic  displacements  for  the  surface  code  in  a  eb&in  of  10  atocs 
are -shown  in  Fig.  6.  " 

The  normal  node-  frequencies' of  the  pariodio  diateni©  chain 
li©  i»  either  the  acoustical  branch,  Di/  £  2c/M,  or  tho. optical 
•branch.,  2a/»  c  o  s -2c{srfK)/^;  ^  Ko-noraal  code,  frequencies  lie  ia_. 
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the  gap  bctwoen  the  bottom  of  tho  optical  branch  and  tho  top  of  tho 
acoustical  branch.  The  semi-infinite  lattice  with  a  light  atom  at 
tho  free  end,  how-vor,  has  its  surface  mode,  frequency  within  the 
"forbidden'’  gap.  In  terms  of  squared  frequencies,  this  surface 
modo  lies  exactly  at  the  center  of  tho  gap  as  shown  in  Fig.  V. 

The  creation  of  free  ends  involves  reducing  a  coupling  constant 
to  zero.  In  accordance  with  Rayleigh's  theorem,  this  leads  to  a 
reduction  in  the  normal  mode  frequencies.  The  surface  mode  in 
particular  arises  from  an  optical  mode  \'nose  frequency  drops  from 
tho  bottom  of  tho  optical  branch  into  tho  "forbidden'’  gap.  . 

Surface  modes -in  the  semi-infinite  diatomic  chain  with  nearest 

and  next-nearest  neighbor  interactions  have  been  studied  by  Gazls 
,  {401 

and  Wallis.  '  When  a  light  atom  is  at  the  free  end,  the  surface 
mode,  may  be  inhibite  d  by  making  tho -next-nearest  neighbor  coupling 
constant  sufficiently  large.  When  a  heavy  atom  is  at  the  free  end, 
however,  a  surface  node  may-  be  induced  by  propor  choice  of  the 
next-nearest  neighbor  coupling  constant. 

The  effect  of  changing  tho  masn  of  the  free-end  atom  in  a  semi- 
infinite  diatonic  chain  with  nearest-neighbor  interactions  has  been 
studied  by  Kori  and  Asahi  '  using  the  scattering  matrix  method. 
Defining  Q  *  (a'-cO/iB  they  find  that  a  surface  mode  appears  with 

* 

frequency  above  the  top  of  the  optical  branch  if  Q<  -$  for  either 
m>15  or  n<H.  For  ti  <  M,  the  surface  mode  with  frequency  in  tho  for¬ 
bidden  gap  merges  into  the  optical  branch  and  disappears  for  Q<  ~J, 
while  for  a  >  H,  it  comes  into  existence  only  for  Q< 

{Vf .  Two  and  Three-dimensional  Monatomic  Lattices 
a,  '  Eoacnstcck-Kewell  model  •  - 

&  relatively  siRplo  lattico-4y«asical  aodbl  which 
admits  ©f  exact  analysis  is  that  of  ftosenstock  aa‘4  Sewell  who 
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consider*  only  ncaroat-noighbor  interactions  but  with  central  and 
non-central  components.  For  the  two-dimensional  simple  square 
lattice  with  free  boundaries,  the  equations  of  motion  have  the  form 


*5jk  “  T<uj,k-rujk>(1-6k,i.^  c<Vi.k"Va'6j.i> 


+  T(uij,k+rujk>(1-6k,N) 


(61  ). 


for  the  x-component  of  displacement  u  with  a  similar  set  of  equations 
for  the  y-component  v  in  which  u  is  replaced  by  v,  the  central  force 
constant  a  is  replaced  by  the  non-central  force  constant  T-,  and  t 
is  replacod  by  o.  As  in  2q.  (50)  the  lattice  sites  vro  identified 
by  integers  j,  k  each  ranging  from  1  to  H  and  the  free  boundaries 
are  doflnod  byj-l,N  and  k-l,N. 

The  exact  solutions  to  the  equations  cf  r< otion  have  the  form 

ujk  -  U  eosOo^^cosOs-Ds^-expCiiot) 

v.s*  - 0 

and 

Ujk  "  0 

-  V  cohvj-^q^cosCk-i)^  exp (i at) 

where  Pj-n^n/N,  P2"n2n,/’N»  and  nl,n2  ar?  *ntlS°rs  tho  range  from 

0  to  N-l,  inclusive.  The  normal  mode  frequencies  are  specified  by 

*  - 

a2  -  (4o/K)[sio2(91/2)+('r/o)sin2(p2/2)3  .  C  «3  ) 

-4»  » 

The  normal  modes  of  vibration  are  clearly  all  wave-like  in  character 
and  ha*©  frequencies  which  lie  in  the  range  0  to  2£(ct+'f)/y]^.  No 
surface  modes,  ‘in  particular,  no  Rayleigh  waves,  exist  for  the  two- 
dimensional  Rosenstpck-Newell  mods.  &  completely  similar  situation 

exists  for  the  three-dimensional  Rosens fccsk-Kewell  model. 

' 


(  62a) 


<  62b) 


I 


*58. 


Tho  question  arise*?  as  to  why  the  Rosens tock^Newell  model 


possesses  ao  Rayleigh  waves.  The  answer  seems  to  be  associated 

(33) 

with  the  fact  that  this  raodol  is  not  rotationally  invariant. 


Tho  lack  of  rotational  invariance  shows  up  when  one  applies  Eq.(43) 


to  the  coupling  constants  involving  a  surface  atom. 


Physically, 


one  can  see  rather  easily  that  a  rigid  rotation  of  the  crystal 


about  a  surface'  atom  will  produce  a  net  force  acting  on  that  atom 


as  a  result  of  the  non-central  forces.  This  is  a  clear  violation 


of  rotational  invariance.  Alternatively,  one  nay  view  the  lack  of  , 


Rayleigh  waves  as  due  to  the  decoupling  of  the  x,  y,  and  z  equations 


of  motion. 


Ludwig  and  Lengoler'*  '  have  shown  how  to  modify  the  Rosens tock- 


Hewell  model  to  make  it  rotationally  invariant.  They  idtrodupe 


additional  non-central  interactions  involving  tho  surface  atoms 


which  restore  the  rotational  invariance.  Thus,  for  the  two-dimen¬ 


sional  case  with  a  free  boundary  dofined  by  k-0,  the  coupling  co¬ 
efficients  involving  the  atom  (00)  are  taken  to  have  tho 'fora 


5(00,  ±10)-  -  (®6  )  ,  5(00,01)-  '(j  j)  (  64  ) 


whore  the  new  coupling  constant  6  must  have  the  value  t/2  in  erder 


for  E<$.(  43  )  to  bo  satisfied.  It  should  bo  noted  that  the  addi¬ 


tion  to  the  coupling  constant  matrix  involving  &  is  antisymmetric 


in  contrast  to  the  symmetric  form  of  the  hulk  coupling  constant 
matrix.  We  also  observe  that  violating  the  condition  of  rotational 


invariance  corresponds  to  the  use  of  incorrect  boundary 'conditions. 
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Although  the  Benge ler- Ludwig  model  makes  possible  the  existence  of 
Rayleigh  waves,  the  simplicity  of  the  Rosens  to.ck-Hewell  model 
associated  with  the  decoupling  of  the  x~,  y~,  and  z-coaponeats  of 
.displacement  in  the  equations  of  motion  is  lost. 

Surface  nodes  nay  rise  in  the  Rosenstoek-Kewell  nodei  without 
the  Leagcier-Ludwig  modification  if  isotopic  impurities  with  suffi¬ 
ciently  -light  mass  are  present  on  the  surface.  The  two-dimensional 
aonatonLc  case  has  been  investigated  by.  ..sahi  and  Hori^36^  who  have 
considered  the  effect  of  replacement  of  either  one  or  all  of  the 
host  atoms  in  the  surface  row  by  an  impurity  atom  cr  atoms.  These 
authors  present  the  conditions  for  the  appearance  of  a  surface  mode 
above  the  bulk  continuum  in  terms  ef  the  mass  defect  parameter 

f  *** 

Q-{H  -U)/M  and  the  force  constant  ratio  t/o.  This  work  has  been 
extended  by  Pukushiaa . 

t  * 

b.  Uodels  with  Rotationally  Invariant  Interactions 

The  simplest  way  to  insure  rotational  Invariance  is 
to  chooso  individual  interatomic  interactions  which  are  rotatiooally 
invariant.  This  can  be  accomplished  by  associating  the  forces  act¬ 
ing  on  the  displaced  atoms  with  changes  in  the  distance  between 
atoms  or  the  angle  between  atomic  bonds  or  other  quantities  which 
reaaia  invariant  under  a  rigid  body  rotation.  A  simple  application 
of  thla  idea  wag  made  by  Gazis,  Berman,  and  Wallis^  who  treatod 
the  monatomic  sisple  cubic  lattice  with  nearest  and  next-nearest 
neighbor  central  forces  plus  angle-bending  interactions  involving 
pairs  of  nearest  neighbors.  For  this,  case,  the  equations  of  motion 
ot  a  bull;  atom  can  be  written  as  ' 
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(65  ) 


where  a  and  p  are  the  central  coupling  constants,  y  is  the  angle- 
bending  coupling  constant,  and  the  other  aquations  arise  by  cycl¬ 
ically  permuting  uTvfw  and  the  increments  $,c  on  the  Indices  4,,s,n. 

The  equatipns  of  motion  must  be  augmented  by  boundary  conditions 
in  order  to  bring  in  the  effect  of  the  surface.  For  a  .{001)  surface 
characterized  by  n-0,  the  boundary  conditions  have  the  fora  . 


L±i  rCu^5,a,-ru4,m,b“teM,iiJ-i) 


*¥  ^\arr^1^o>S{sw,*)»^wiai.^.:j-0 

«£L  [ e^Vt,»+6,-l~v4,ia;t.r^^a^ ,-.1) 


+Y  ^Cvt.r.,-rvt,m,o)‘4^fm^rp^tt^6,-l>3  }  -*  ;:V: 


+?6?-ut+6  ,s,o^^iiB,~l’5'“4.,5s+4,o"vi,a;^i,-l^  J 


C  65  ) 


•-  ~  *  V  > 

3*? 


ig^gUP 

PUP 

iT$ . 

MxL 

mmm 


:  :--^gjc  **S&£g 


JjftAdsi'i  ££**  vir?:^v>  X' 


"••  v- 


’M&m& 


■  V’C. 


.  .'/if 


r  -Ml 


0i 


1 


&V 


aa;. 


•  J.  ...  -  ...,'..  ■  ....  .••  ..•  ...  -  •  :  v.,\.  .  •  •  ..•••  *-.,.  irvv^7  .--V--- -i*  '■•4vv\r'-~  '■  • ■ 

iv.  a,  ^-v-viy-  ,;•  •  •••■•.q.v  ^ ^ Va * ■ 


lllSlillt> 
-jpplilf ; 
'  IgSillli'. 

•  lllilf  i 
'fH&i 

SslslspeS; i 

,  '3SB§i'f 

llSIlPl 

lllillt I 


..'  fn* 

JpillS.: 


-'  m 


'  -  J  x" 


'  ‘  ;  V  v. 

a-*  .  _  *^a|5c£S&  * , 

iHSI*. 


^JSS 


~i~:  •  <r" ' .-  •■ ! 


-  '  ✓;  ,‘vv  i^v 


..  .  .  ,  .  „ 

a"  ..-<V-^:W  t>'^- 


one  of  tho.  attenuation  constants  at  a  critical  wave  number  <p  -  1.75. 

c 


At  <p  this  attenuation  constant  switches  iron  a  real  value  to  a 


complex  value  of  tho  fora  qo+in. 


Another  case  which  has  been  trow ted  with  angle-bonding  inter- 

.no  4  (?  Mia  ttAaafa«n<a  Ka  .1  f.aanf  avarf  A  Cd4  i  PfcSSG 


actions  is  the  monatomic  bouy-eentered  cubic  lattice. 


velocities  and  attenuation  constants  were  calculated  for  surface 


waves  propagating  in  the  [IQOj  direction  of  the  <001)  surface  of 


vanadium,  which  exhibits  ordinary  Rayleigh  waves,  and  iron,  which 


exhibits  generalised  Rayleigh  waves.  Dispersion  of  tho  phase  vel¬ 


ocities  was  established,  but-  so  divergence  of  an  attenuation  constant 


was  found  in  cither  case . 


Using  the  procedure  described,  above,  Kaplan^45  ^  has  studied 


the-  semi-infinite  monatomic  square  lattice  with  nearesfe-and  next- 


nearest  neighbor  central  .forces  for  the  Situation  where  a  homogeneous 
isotopic  impurity  layer  is  preheat  ou  the  surface  edge.  He  found 
that  the  impurity  layer  has  little  effect  on  long-wavelength  Rayloigh 
modes .  This  is  to  ba  expected  because  long-wavelength  Rayleigh  modes 
penetrate  deeply  into  the  crystal.  At  the  Brillouin  zone  boundary, 


however,  the  Rayloigh  modes  are  localized  within  an  atomic  layer  vv 

v-  ^ 

so  of  the  surface,  asd  they  should  bs  significantly  affected  by  the 


impurity  layer.  This  is  just  what  Saplac  found.  "For  a  sufficiently' 


heavy  impurity-atom  tsass  li' ,  there  are  two  low-frequency  surface 


modes,  one  of  which  is  tho  Rayloigh  mode.  As  the  impurity  mass 

? 


decrease's,  ia  succesaioa  the  <sxtra  low-frequency  surface  »edo  dis¬ 


appears,  a  high-frequbadj;  surface  modi  appears  above  tho  allowed 


hand  of  frequencies,  tho -Rayleigh  aeda  disappears,  and  finally  a 

1  -  ■  - 


second  high  frequency  soda  apj^ars.  .< 
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Surface  nodes  of  vibration  in  monatomic  faco-centercd  cubic 
lattices  have  been  investigated  recently  by  Alloa,- Alldredge,  and 


Detfette' 


who  assumed  central  interactions  of  the  Leonard- Jones 


6-12  typo  botweea  all  pairs  of  atoms  in  plhtes  up  to  21  atomic  layers 
thick.  The  parameters  characterizing  the  pair-wise  interactions 
wore  taken  to  be  the  same  at  the  surface  as  in  the  interior  of  the 
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crystal.  The  model  nay  be  considered  appropriate  for  a  slab  of  a 
rare-gas  solid . 

Allen  et  al^4®^  first  considered  surfaco  modes  on  a  (111) 
surface  with  propagation  in  the  rX(£llO])  and  ns(Lll23)  directions 


as  well  as  in  the  MK  direction  along  the  edge  of  the  Brillouin  zone. 
By  solving  £q.(49)  on  a  computer,  they  found  the  expected  dispersion 
of  the  Rayleigh  nodes.  They  also  found  additional  surface  nodes  of 
non-Saylcigh  typo  -whose  frequencies  lie  in  gaps  in  the  spectrum  of 
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the  bulk  modes  as  shown  in  Pig-.  9.  The  atomic  displacements  of 
tho  non-Rayleigh  surface  nodes  are  primarily  parallel  to  the 

surface,  whereas  those  of  Bayioigh  nodes  ar®  prinarily  pex-pon- 

*  •  *  (471 

diculer  to  tho  surface.  In  a  second  publication,  Allen  et  al1  J 

report  tho  existence  of  surface  modes  on  a  (100)  surface  which 

have  the  same  frequencies  as  bulk  nodes  with  tho  some  value  of 

the  two-dimensional  wave  vector  parallel  to  tho  surface.  Wo  have 

already  encountered  surfaco  nodes  of  this  type  in  the  continuum 

theory  of  materials  exhibiting  generalized  Rayleigh  waves. 

In  a  recent  publication.  Alien,  Alldr  edge,  and  BcWetto^®  ^ 
have  extended  these  calculations  to  the  hexagonal  close-packed 
crystals.  Her®  again,  they  find  surfaco  modes  in  gaps  in  the 
hulk  frequency  spectrum. 
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45. 

*  *  (4Q\ 

Those  some  authors  havo  recently  given  *  a  general  forau- 
lation  of  the  surface  vibrational  mode  probloa.  They  discuss  tho 
circuastancos  under  which  the  displacement  ellipse  of  the  particlos 
lies  in  the  sagittal  plane  (surfaces  with  "axial-invorsion  symmetry") 
and  when  surface  nodes  can  havo  frequenci.es  lying  within  bulk  node 
subbands. 

Detailed  calculations  havo  boon  reported  by  Allon  ct 
for  monatonic  face-conterod  cubic  crystals  with  (111) ,  (100) ,  and 
(110)  surfaces.  They  find  a  number  of  surfaco  mode  branches  for 
each  surface-five  for  the  (111)  surface,  at  least  nineteen  for  the 
(100)  surface,  and  ton  for  tho  (110)  surface.  These  results  aro 
probably  dependent  to  some  extent  on  the  model  employed  —  inter¬ 
actions  between  all  pairs  of  atoms  through  a  Lcnnard-Jones  potential. 
Sono  of  the  surface  modes  found  by  Allen  ot  al^50^  are  primarily 
localized  in 'the  second  layer  from  the  surface  or  even  a  deeper 
layor,  rathor  than  in  tho  surfaco  layoj:  itself.  These  various 

surface  nodes  can  be  regarded  as  "peeling  off"  in  succession  from 

**  • 

the  bulk  branches.  The  effects  of  surface  Xorco-conutant  changes, 
were  also  considered. 

Soso  •  calculations  emphasizing  the  long  wavolength  region 
have  boon  carried  out  by  Alldrodga,  Allen,  and  DeWatte*53^  for 
the  (100)  surface  of  a  face-centered  cubic  crystal.  In  addition 

ip 

to  the  usual  Rayleigh  waves  and  pseudosurfaco  waves,  they  find 

an  additional  surface  wave  propagating  in  tho  £110]  direction 

which  is -polaris s J  shear  horizontal.  Alldredgo*525  has  shown, 

however,  that  the  attenuation  constant  .q  -  k2  and  hence  vanishes 

^  P 

in  the  long-wavelsagtfc  lisit,  so  this  mode  is  not  a  Bayleigb«typ» 
surfaco  wave.  f 


46. 


An  investigation  of  size  and  surface  effects  on  the  lattice 

dynanical  properties  of  small  particlos  has  been  reported  by 

Oickoy  and  Paskin^  ^  who  used  the  Lennard-Jonos  potential 

« 

and  various  crystal  shapes.  The  frequency  distributions  for 
bulk,  surface,  and  odge  modes  were  calculated  and  woro  found  to 
exhibit  peaks  occurring  at  successively  lower  frequencies  in  the 
order  given. 

The  case  of  an  isotopic  impurity  layer  on  the  (111)  surface 

of  a  f&ce-ccntcrcd  cubic  crystal  has  also  been  investigated  by 

Allen  et  al.^4^47^  As  in  the  work  of  Kaplan^45  *  they  found  for  a 

light  impurity  layer  that  high  frequency  surface  modes  nay  appear. 

above  the  bulk  continuum  and  that  the  Rayleigh  node  nay  disappear 

at  sufficiently  large  values  of  the  wave  vector.  This  work  has 

(  541 

been  extended  by  Alldredge,  Allen,  and  Pe^tte  .  to  the  (100)' 
and  (110)  surfaces.  For  very  light  or  very  heavy  adsorbed 
particlos,  there  are  three* principal  surface-mode  branches 
associated  with  the  adsorbed  layer,  but  for  non-extreme  nass 
ratios,  the  separation  of  the  principal  surface-mode  branches 
is  not  complete.  They  find  that  heavy  adsorbed  layers  enhance 
the  localization  of  Rayleigh  waves  and  that  light  layers  have 
the  opposite  effect. 

Phonon  frequency  distributions  for  layered  structures 

f  rs"i 

r  havo  been  studied  by  Allen,  Alldredge,  and  EsTFstte''  *  who 
find  that  interface  effects  are  small  compared  to  surface  effects. 
Soae  remarks  are  made  concerning  superconducting  transition 
temperatures  in  thin  films. 
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A  treatment  of  an  adsorbed  layer  on  the  (001)  surface  of  a 
simple  cubic  crystal  has  been  given  by  Dobrzynski  and  Mi  11s ^56) 
using  tho  Roscnstock-Nowoll  model.  When  the  mass, ratio  of  ab¬ 
sorbed  and  bulk  atoms  is  less  than  unity,  they  find  localized 
nodes  above  the  bulk  continuum  whereas  if  the  ratio  is  greater 
than  unity,  they  find  localized  modes  bolow  the  bulk  continuum. 

Tho  latter  are  not  Rayleigh  waves,  however,  since  the  Rosonstock- 
Hewoll  model  docs  not  exhibit  Rayleigh  waves. 

Models  *0»ieb,  involve  central  forces  only  must  satisfy  Cauchy's 
relation;  c..g-  (cubic  crystal,  for  example).  Many  metals  do  not 
satisfy  Cauchy's  relation,  so  non-central  forces  are  required  for  a 
proper  ip ''Aon  of  their  lattice-dynamical  properties.  Recently, 

thiRfits  and  Rieder  1  have  studied  surface  waves  propagating  in  tho 
£00-0  di<i'«t;Vion  oa  fch«  (100)  surface  of  a  monatomic  face-centered 

cubic  crystal  using  neaxest-and  next-nearest  neighbor  central  forces 

'  n* 

plus  non-central  forces  between  nearest  neighbors.  Rotational  ia- 

.variance  wa£  preserved  using  the  procedure  of  Ludwig  and  Lengeler. 

Musser  and  Riisder  found  that  the  frequencies  of  long-wavelength 

•*  * 

Rayleigh  waves  wre  significantly  affected  if  rotational  invariance 

A 

* 

was  ignored. 

The  effect  of  changing  the  value  of  tho  coupling  constant  for 

* 

the  central  interaction  between  a  surface  atom  and  its  nearest 
neighbors  in  the  second  layer  has  been-  studied  by  liuacer  and  Riedor1'^7^ 
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using  a  Green's  function  technique.  If  the  surface  coupling  constant 

is  smaller  than  the  corresponding  bulk  value,  a  second  lev- frequency 

surface  mode  may  appear  for  wave  vectors  near  the  Brlllouin  zone 

boundary.  The  Rayleigh  mode  frequency  is  also  depressed.  If  tho 

surface  coupling  constant  is  larger  than  the  bilk  value,  a  high- 

frequency  surface  mode  may  appear, above  the  bulk  continuum.  These 

C  11 

resul1-  are  entirely  consistent  with  Rayleigh's  theorem.  ' 

Uusser  and  Rieder  also  find  surface  modes  with  frequencies  within  .a 
gap  in  the  bulk  continuum. 

Soveral  workers  have  studied  localized  modes  associated 


with  isolated  atoms  or  molecules  adsorbed  on  crystal  surfaces 

.(53) 


using  throe-dimensional  models.  Grimley considered  an 
ndatem  on  a  (100)  surface  of  c  Rosens tock-Newoli  model  of  a 
simple-cubic  lattice  and  obtained  an  expression  for  the  localized 
mode  frequencies  using  the  Green's  function  technique  of  Montroll 
and  Potted An  adsorbed  diatomic  molecule  and  the  inter¬ 
action  between  a  pair  of  adatoms  were  also  treated.  Specific 
application  was  made  to  hydrogen  and  to  carbon  monoxide  adsorbed 
on  platinum.  The  interaction  between  two  adatoms  has  also  boen 


discussed  recently  by  Cunningham,  Dobrzynski,  and  Uaradudin 


(60) 


A  problem  related  to  those  just  described  -  namely,  that  of 
an  isotopic  impurity  atom  located  in  or  near  a  (100)  surface  of 
a  simple-cubic  lattice  -has  baea  investigated  by  Ashkin^  61  ^ 
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The  sodol  employed  consists  of  nearest  and  next-nearest  neighbor 
contral  forces  with  the  force  constants  chosen  to  give  elastic 
isotropy  in  .tho  long-waveiength  limit.  Utilizing. tho  Montroll- 
Potta  Green's  function  technique  and  by  assuming  that  tho  local 
jsodo  frequency  is  large  compared  to  the  maximum  frequency 
of  tho  perfect  lattice,  Ashkin  obtained  the  following  expression 
for  the  square  of  tho  localized  node  frequency 


-  ,  s  2  a^(» ,t)  *2^°*^)  Xa3(e?,t)  9 

- +  ^Anryr  +  ~3 - 2  +  o<x2) 


i{(a,iy  +  i&^V'Dy 


<  73  > 


Where  or  danotos  the  polarization  of  the  localized  node,  l  opacifies 
the  site  of  tho  impurity  atom,  X  «  M’/K  is  the  ratio  of  tho  impurity 
and  host  atom  masses;  and 
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(74b  ) 
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(74c  ) . 

In  Eqs.  (74  ),  ^o9 ( ■W.  0  is  an  element  of  the  dynamical  matrix  *  for 
tho  crystal  with  *  surface  but  without  |ho  impurity  givon  by 

"  C-/H)4fraUt!)  where  *ffg(-M,')  is  ea  element  of  the  force 


constant  matrix..  Specific  results  aro  gives  by  Ashkin  for  the 

f 

modal  mentioned  above.  Qualitatively,  the  localized  mode,  which 

-+  ^ 

is  triply  degenerate  if  the  impurity  is;  in  the  bulk,  is  split  by 

tto  surface  into  a  singly  degenerate  mode  polarized  perpendicular 

.  *  •*  >>  -  •  ^ 
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to  tho  surface  and  a  doubly  dogenouoto  mode  polarised  parnllol  to 


the  surface.  As  wight  be  expected,  the  localized  mode  frequency 


is  decreased  when  the  impurity  is  in  the  surface  layer.  Thus, 


for  X  •*  0.2,  the  porpendicular  mode  has  a  frequency  29  por  cent 


lower  and  the  parallel  nodes  §  per  cent  lover  than  the  localized 


node  frequency  when  the  inpurity  atom  is  in  the  bulk.  For  this 


case  the  bulk  localized  node  frequency  is?  essentially  achieved 


when  the  inpurity  atoa  is  coved  just  one  layer  inward  iron  the 


surface. 


.The  work  of  Ashkin  has  been  extended  to  body "centered  cubic; 
lattices  by  Dobrzynski/62^  who  used  the  nodol  of  Clark.  Cazio,  and' 

■»■  /r*)N 

Wallis  to  consider  both  single  adatonn  and  an  adlsyer  on  a 
(100)  surface. 


A  theoretical  investigation  of  vibrational  edge  modes  has 
recently  been  carried  out  by  Haradudin,  Wallis,  Hills,  and  Ballard^64^ 
who  treated  an  edge  of  a  cubic  elastic  medium  bounded  by  two  (100) 
faces.  Jfor  the  case  of  elastic  isotropy  with  tho  Laisd  constants 
X,  u  equal,  they  found  a  localized  edge  mode  with  speed  0.S013  C* 


compared  to  tho  Rayleigh  wavo  speed  of  0.9194  C. . 

t 


(Vi),  Two-  and  three-dimensional  Diatomic  Lattices  '*1th  Short- 


range  Interactions, 


Ordinarily,  diatomic  crystals  have  ionic  character, 
and  the  long-range  Coulomb  interaction  between  the  iocs  must  bo 
taken  into  account  in  a  lattice-dynamical  mode.  Some  qualitative 
features,  however,  aay  be  studied  using  a  model  with  only  short- 
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tango  interactions,  A  general  formulation  of  the  problem  based  on 
a  Green's  function  method  was  given  a  number  of  years  ago  by 

LifShitZ  and  BOSGnZWOig^^  »hn  rin-tnrt  fh*  Av4af/4t\Aa  -n-9  hnfrk  n^nniw 


who  noted  the  existence  -of  both  acous¬ 


tical  and  optical  surface  jnodos.  Wallis^ 


later  studied  the 


Rosenstock-Newell  model  for  finite  two-  and  three-dimensional 
diatomic  lattices  of  the  rock-salt  type  using  both  a  perturbative 
procedure  and  a  Green's  function  procedure.  It  was  found  that  8 


4,/  i  -/i*’’-  -  ‘‘""t 

'jwL£  | 


•  "■  *  • 


band  of  surface  mode  frequencies  occurs  in  the  gap  between  the 
acoustical  and  optical  branches.  These  modes  are  closely  re¬ 
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lated  to  the  surface  mode  of  the  linear  diatomic  chain  with 
nearest  neighbor  interactions.  Furthermore,  for  the  finite 
lattices  considered,  edge  and  corner  modos  were  found  in  which 

e> 

the  atomic  displacements  are  localised  near  edges  and  corners, 

*  f 

respectively.  Typically,  the  frequencies  of  edge  modes  lio  colow 

«£ 

those  of  the  surface  modes  mentioned  above  and  the  frequencies  of 

*  ~i 

corner  modes  lie  r  below  those  of  the  edge  modos.  This  is  physically 

r 

reasonable  in  view  of  the  decreasing  number  of  interactions  with 

s  .  » 

neighbors  possessed  by  an  atom  as  it  changes  from  a  surface  atom 

*  « 

to  an  odgo  atom  to  a  corner  atom. 

A  surface  mode  calculation  using  s  rotationally  invariant 

£  * 

model  for  a  rock-salt  type  crystal  was parried  out  by  fallis, 

Hills,  and  Jiaiadudin®7^  who  found  a  surface  branch  in  the  gap 

between  the  acoustical  and  optical, branches, 

*  .  .  * 
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(vii).  Three-dimensional  Diatomic  Lattices  with  Long-range 
Interactions. 

Sines  diatomic  crystals  typically  have  ionic  character, 
a  realistic  lattice-dynamical  model  must  include  the  Coulomb  inter¬ 
actions  between  the  ions  as  well  as  interactions  associated  with  the 
polarizability  of  tho  ions.  The  long-range  character  of  the  Coulomb 
interactions  makes  the  calculation  of  the  surface  mode  frequencies 
somewhat  difficult  because  the  effects  of  tbe  surface  are  not  local¬ 
ized  within  a  few  atomic  layers  of  the  surface.  Furthermore,  tho 
possibility  of  retardation  effects  associated  with  the  .electromagnetic 
field  introduces  an  aspect  to  tbe  problem  not  encountered  when  only, 
short-range  interactions  are  considered. 

Wo  can  formulate  the  surface  mode  problem  with  Coulomb  Inter¬ 
actions  by  simply  noting  that  the  coupling  constant  matrix 
in<)(4Kf-l>x')  now  can  be  written  as  tbe  sum  of  a  Coulomb*' part  and  a 
short-range  part; 


*  The  Coulomb  part  has  the  fora 

-.2 


(75  ) 


lafl(4«,4‘a’}«— 2 - 

W  &r0UK)arpa 


_  / - 


£uUx) }  «  o 
(76  ) 


where  are  the  electrical  charges  on  ions  of  type  x,H‘ , 

respectively.  The  short-range  part  can  be  treated  in  terms  of  central 
forces  and  angle-beading  forces  as  was  done  earlier.  The  calculation 


o t  tho  surfaco  mode  iroquencios  then  requires  '.ho  solution  of  tho 
seedier  equation,  Eq.(  49),  using  tho  coupling  constants  specif iod 
by  £q.(75  )'.  A  major  task  involved  is  tbj  evaluation  of  tho  lattice 
suns  in  the  Coulcab  contribution  to  ihu  reduced  dynamical  matrix 
indicated  in  Eq.(  48). 

Upon  solving  the  secular  equation,  oco  obtains  acoustical  sur¬ 
face  modes  which  reduce  to  the  Rayleigh  waves  in  the  long-wavelength 
limit.  Of  more  Interest  in  the  present  context,  however,  are  opti¬ 
cal  surface  modes  which  show  significant  effects  duo  to  tho  Coulomb 
interactions. 

Fuchs  and  Kliewer^8^  carried  out  a  calculation  for  a  slab  of 

a  NaCt-type  crystal  with  a  (001)  surface  for  both  rigid  ions  and 

polarizable  ions.  The  short-range  interactior.3  were  taken  to  be 

nearest-neighbor  central  forces.  Fuchs  and  Kliewer  neglected  the 

changes  in  the  short-range  coupling  constants  and  the  relaxation 

of  the  ions  near  the  surfaces.  In  the  limit  of  long  wavelengths 

(kp  small  compared  to  the  reciprocal  of  tho  lattice  constant  a), 

tho  lattice  sums  can  be  replaced  by  integrals  and  the  equations  6f 

motion  become  a  set  of  integral  equations  which  wore  solved  by 

Fuchs  and  Kliewor.  However,  a  simpler  procedure,  as  pointed  out 

by  Kliewor  and  Fuchs^69^,  is  to  treat  the  p:  ->blera  aacroscopically 

in  terms  of  Maxwell's  equations  and  the  appropriate  boundary  con- 

^  • 

ditions.  The  effects  of  retardativn  can  be  Included  in  a  straight¬ 
forward  fashion. 

Consider  a  slab  of  thickness  L  -  2d  with  surfaces  normal  to 

*  *. 

the  x-o*rection  at  z  *  id.  The  material  is  assumed  to  bo  optically 

l  - 

isotropic  with  a  dielectric  constant  €(gj)  taken  to  no  that  of  tho  • 

bulk  medium.  The  electric  and  magnetic  vectors  aye  chosen  to  have 

tho  iota  of  plan??  waves  propagating  along  the  surface  with  vr 

-•  y,  . ..  -  -■  V 

vector  k  and  frequency  e;  ~  4 
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Kliewer  codes  are  transverse  In  the  sense  that  v-u  -  0  and  that 
at  k  *  ui^/c,  the  nodes  are  no  longer  loctflined  at  the  surfaces. 

The  surface  nodes  just  described  are  non-radiative  aodesf 
The  fields  drop  off  erponeatiully  outside  the  slab  and  do  not 
interact  directly  with  external  radiation.  Other  surface  nodes 

ifiO) 

exist,  so-called  radiative  nodes 


,  with  iace&ir.g  ov  outgoing 


waves  outside  the  crystal.  These  nodes  are  not  tree  noraal  codes 
and  arc- highly  caaped  even  when  €{&>  is  real.  The  frequencies  o? 
the  radiative  codes  lie  above  ths  light  line  ®  -  ck-»  in  contrast 
to  the  nor^  radiative  codes  which  lie  belies  this  line.  . 

Ac  independent  treatment  the  optical  asode's  of  vibration  of 
ionic  crystals  finite  in  one  or  core  dimensions  Isis  bees  given  by 
EngXnan  and  Huppin^^in  a  series  of  papers.-;  They  genaralixe 

the  results  of  Puchs  and  JQ  lower  to  crystals  of  diverse  shapes 

•  s':-.,-  -  -  -  -  -  ,-  - 

using  the  sacs  assueptions  as  Puobs  and  Xliewbr.  In  'the  cas*  id 

a  club,  their  results  agree  with  those  of  2? uchs  and  Xli^gssL  X« 
the  oas o  of.  a  ephora,  the  surface  codes  fall  inte.  tyb'-typei.  Tin 
so-called  ‘'electric  codes have  sg  radial  cenpaueat  -bf/tfea-sas^ V_ 


netic  vector  and  have  frequeixtien  specif  X£dh?  :. 

-  RJlk^Cy&S'  *-  _> 

•shore  j  t  and  b,  arc  spherical  Eessol  Hap-kCi  f u.sj  t  ions , 

v  v  .  ■  *  _  -  -  -  ;v  ~  ^ 

kQ  •  £j,“ie/c ,  k,  ”  €^s/cr£  and  are  the  deles  trie  ecnstur-.tS  of 
the  sphere  ai»4  surreuading  ,»ediUb7  *ad-E  is  tlia  yadits.  The  ^ahg*-  ' 
no  tic  codes,’*  oh  the  other  hand,  have  no  radial  eosiposeai  c-f  tbb 
electric  vectSr  sta  navs  trequejKieb-'specifleS  ,by 

^(kjBHk^ayoV  -  ..  (mb} 
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In  the  non-retarded  Halt,  the  surface  node  frequencies  of  the 


electric  type  are 


2  ,  .  2  I  | 

wso  t»T  ^  €o)+€y  C?flT7t  1 


r  4  “*  I  }2  ;3  )  «  •  «  i 


<  89 : 


whore  £  and  €„  are  the  static  and  high  frequency  dielectric  con¬ 


stants  of  the  sphere.  The  amplitudes  of  vibration  viry  with  dis- 


tancs  r  from  the  center  of  the  sphere  as  rv  .  The  node  t“l  is 

(721 

**  nWfta  *■>  W  ^  v  In  Vine*  4  a/1  Kv#  1?n*HKT  ^  aK' 


therefore  a  Ufcifdrn  sso4e  vhich  has  bears  studied  by  FrBhlich; 


Sines  there  is  a  macroscopic  electric  fiald  associated  with 


surface  optical  phonons  of  the  Fuchs~5CTiewar  type,  ono  may  expect 
that  these  modes  will  interact  with  any  free  charge  carriers 


*  C*7*l 

present  in  the  crystal.  Indeed,  it  has  been  shewn  by  Kheifots'  ** 


Suppin^7^ ,  by.  Chiu  and  Quinn^7® ,  and  by  Kallls  and  Brion^7^ 


that  the  surface- optical  phonons  coupio  to  thG  surface  plesaon*. 


?er  carriers  in  a  simple  parabolic  energy  band,  one  can  utilize 
t£e  procedure  presented,  above  r.cspt  that  the  diolectric  constant 


.given  by  £q.  (84)  is  roplacod  by 


/  2.  2 
{  y?***u/  — ,  uj  * 

s{u»)  »  |  -g~z  ”  rg  | 


2  -  5s  *  ■'Tj6  ~  *  »  • 

whstc  -*  vr  ^nno  '/^E  is  the  ,pl*sea  frequency ,  ^  ia  tho  carrier 


concha tratioa  und  e'  their  effective  jsase.  Xa  the  isa^-rotardod 


"liait  *  relatively  s^sple  cxpr.es^isa  cast  bt  gjbSsK  fsjr  the  coupled 
.sus'iaco' freqtsaecies^f  -a1  very  thi&s-'sCU^,-'-  '  3'  ' 


sSoyh  - givea  byvSql:  ( •  ^5-’J  /and  givan  hy ;. 


■M  '.-i  '33  -3  :•=•  3  H;  -V  -'"3 

ct.v 7-X>'.: J; '-  .  3;  ; 


is  the  surface  plasmon  frequency.  Experimental  verification  oi 


Eq.  £  91  )  has  recently  been  reported  by  Reshlna,  Gerbshtein,  and 

«  «  TnCh  Anc  n#  + 


for  n  -  InSb.  Theoretical  discussions  of  the  effect 


of  an  external  magnetic  field  on  the  coupled  surface  modes  has 
been  given  by  Jluppin^74^ ,  by  Brion,  Wallis,  Hsrtetein,  and 


Burstoin^^end  by  Cbiu  and  ^juian.^*^ 


There  are  a  number  of  rather  interesting  physical  phenomena 


associated  with  Fuchs-Kl lower  surface  modes.  Mahan 


has  shown 


that  the  faniliar  problem  of  image  charges  can  be  discussed  in 


terms  of  surface  .polaritons.  if  a  source  charge  is  located  near 


the  interface  between  two  insulating  dielectrics  or  a  dielectric 


and  a  vacuum,  the  resulting  potential  field  can  be  represented  by 


&  sot  of  image  charges.  According  to  Mahan,  surface  polaritons 


make  a  notable  contribution  to  the  iaat,e  charges';  indeed,  they 


give  rise  to  the  entire  polarization  field  if  the  source  charge- 


is  outside  the  dielectric. 


A  second  phenomenon  associated  with  surface  polaritons  is 


the  Van  dor  Kaala  interaction  between  solid  bodies.  It  has  been 


shown  by  Van  Xampen,  Kijboer,  and  Scbram,,Si^  (noa-retarded  case)  and 


by  Gerlnchr  ’  (retarded  case)  that  tfescVan  der  Saal£  interaction 


energy  is  simply  th£  suro-poSot  energy  of  the  surface  polar! toss. 


That  only  the  surface  modes  and  E3S  the  bulk  modes-  arov involves 
is  clear  because  only  the  surface  modes  h&vs.  fields-  t-atenaiag  o<££~ 


side  the  media. 
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.  will  bo  recalled  that  Fuchs  and  Xliowor  made  two  important 

approximations  in  thoir  treatment  -  namely,  they  neglected  tho 

effect  of  the  surface  on  the  short-range  interactions  and  they 

assumed  the  wavelength  large  compared  to  the  lattice  spacing. 

Lucas^  ^  has  avoided  the  first  of  those  approximations  in  an 

investigation  of  optical  suriaco  modes  in  thin  slabs  of  an  NaCl-typo 

crystal  with  (001)  surfaces.  Koarest-noighbor  central  forces  wore 

assumed  for  She  short-range  interactions,  and  tho  absence  of  neighbors 

beyond  the  boundaries  was  taken  Into  account  in  the  equations  of 

motion  for  the  surface  atoms.  However,  relaxation  of  the  ions  noar 

tbo  surfaces  and  retardation  were  neglected.  The  ions  were,  assumed 

to  be  rigid,  and  the  appropriate  Coulomb  sums  were  evaluated  in 
/  84  ^ 

plhne-wise-feshiea:  '  Since  each  lattice  plane  contains  equal 

cumbers  of  positive  and  r.ogativo  ions  and  is  overall  neutral,  the 
interactions  between  planes  decrease  exponentially  with  separation, 
-s»d  tha  convergewse  of  the  intorplsner  interactions  is  therefore 
very  rapid.  r 

Luesa  carried  out  specific  calculations  only  for  k  -  0  where 
k  is  the  two-dimensional  wave  vector  parallel  to  tho  surface.  Ho 
found  a  transverse  optical  surface  stodo  (riispl3cooants  parallel  to 
the  surface)  with  *  frequency  very  slightly  below  that  of  the  trees- 
versa  optical  bulk  mode  for  k  «*  SJ.  This  result  is  in  -on trast-  to 
that  of  Fuchs  and  Slfdwor  who  found  that  tho  transverse  optical 
msdc  is  not  Jocciimod  at  tbs  surface  for  k  -  0.  Thin 
difference  is  attributable  id  the  inclusion  by  laicaa  of  corrections 
-  tc  ths-  short-range  forces  'actiBc..oc  surfcaa  Ions  due  to  their 


saalXer  cusbbr  of  neighbors,  "...  However*  Eucrs  found  so  .Jea^iisdlsiX 


60. 


optical  surface  nodes  (displacements  perpendicular  to  the  surface) 

for  k  *  0  in  agreement  with  Fuchs  and  Kxiower. 

~P 

A  wore  refined  calculation  for  the  NaCl  slab  with  (001)  sur¬ 
faces  has  been  carried  out  by  Tong  and  Maradudin^  30  ^  who 
assuaod  nearest-neighbor  central  forces  and  Coulomb  forces  between 
rigid  Ions.  Proper  account  was  taken  of  the  number  of  short-ranf/i 
interactions  of  surface  atoms,  and  the  possibility  of  changes  i,i 
tho  intorplanar  spacing  near  the  surface  was  included  in  tha  formu¬ 
lation,  but  retardation  was  neglected.  Calculations  were  made  fer 

wav©  vectors  k  covering  tbo  entire  two-dimensional  first  Brillouin 
~P  . 

zone. 


If  wo  neglect  tho  doubling  of  the  surface  codos  duo  to  the 
presence  of  two  surfaces,  then  we  can  characterize  the  results 
of  Tong  and  Maradudin  as  follows.  For  propagation  in  tho  [100] 
direction,  there  are  four  surface  branches.  One  branch  is. 

4 

acoustical  in  character  and  corresponds  to  Rayleigh  waves.  There 

are  three  optical  surface  branches.  A  high  frequency  branch  lies 

bo tween  h>L  and  «T  and  approaches  a  frequency  somewhat  bolow  s«L 

as  k  -  0.  This  branch  would  be  significantly  affected  by  re- 
P 

tardasion.  A  low  frequency  branch  lies  below  the  bulk  optical 
mode  frequencies  for  small  kp,  approaches  a  frequency  just  below 
aij  as  kp  -  0,  end  is  polarized  in  the  sagittal  plane.  An  inter¬ 
mediate  frequency  branch  has  the  same  frequency  at  kp  -  0  as  tbo 
lew  frequency  branch,  but  rises-  into  tho  z-egion  between  and  a^, 
and  is  polarised  parallel  to  the  surface.  Tong  and  ilaradudin  cal¬ 
culated  tha  deviation  of  the  slab  frequency  distribution  from  the 
bulk  distribution  and  found  pronounced  peaks  corresponding  to  tho 
surface  branches. 
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The  investigations  of  Fuchs  imd  Klfsrcr  and  of  Ton?  and  Maradudin 


do  not  lesd  to  on  entirely  consistent  picture  of  the  surface  optical 
modes  of  an  ionic  cryzls 1  sleb.~  The  problem  of  reconciling  these 


results  has  been  Und^rtakon  by  Tories  asd.?ech3‘  °a  J  end  by  Chsn. 


Alien,  Ail dredge,  sbd  Do  Ifetie^?8''  ,  The  following  picture 
onerges?  cs  shown.  in  Fig-,  11c  Tito  high  frequency  branch  of  Tong 


and  Saradudirs  splits  as  h^»  o  into  ivo  breaches,  one  approaching 


53id  the  other  s/^.  These  are  the  Tschs-Klioudr  modes.  (Tong  and 


.Karadudin  had  not  u„cd  a  sufficiently  email  kp  to  reveal  tho  split* 


tisg.)  Hear  k^  "  Q,-  tho  ?uchs=-Klien-i,r  eodas  onto?  ths  bulk  IX)  end 


TO  ©patiaus  end  bsscdae  ”sixa4“  o?  pseudoeurfaca  waves,  The  Jo* 
frequency  brsach  enters,  the-  sogfisticai  bulk  continuum  about  midway 
to  the  Brillouir.  so as  boundary  and  becomes  a  alsud  soda  branch. 

The  islerasuia>g'.-braneh.  does  not  couple  to  hulk  modes  end  resales 

a  pure  surface  breach  oat  to  tfes>  seise  boundary, 

~-\  -  -  -  {§5  871  ----  - 

vhsn  oi  al  ’  have  investigated  several  directions  of 
gropegstioe  and  kayo  found  a  ^aria iy  oi  surface  modes  in  gaps  in 


the  bulk  frequency  spectrum,  They  have  also  utilised  this  shell 


jh tidal'  fox  in  ??aCX  ol&b  and  fiud  that  the  Fa'ebs-Elisser  sc-des  are 


pSguddsurfsca  **¥©#-  for  all  jsav©  vectors.  -  The  other  surface  optical 
KOdes  ars  etill-t^s-murfaco  spdaa  hssr  4s^  «  Q.  T|ie  ahel.l-^odel 

epeu^tions -3»ys  «s»i^ed(8$\  to  a&y  snd  SfeCt,  rock-salt 

structure  sulerisls-hn^iBg  tfco  largest  gaps,  bbt^esjy  acou^tit*! .  s&d 
optical  branches.  Veil  defined  '‘sir.roseopie’'  s«rface4u?de  fesisches^3?) 
with  frequencies  lyinz  witbia  th^  gap  worn  found  for  .the  <$Gii  ./ 
surfaces.  Surface  sodes  associated  *ith  the  (601) surface  of  Id? 


have  also  been  studied'89'  using  th©  shell  model.. 


Tho  problem  of  surface  modes  in  a  oeai-infinito  rocksalt 

(  qr>\ 

crystal  has  bee;;  tackled  by  Acbsr  end  Sarsch'  '  using  the 
sea fetering  matrix  method*  However,  their  results  appear  to  bo 
in  error  * 

Another  theoretical  investigation  of  the  ionic  slab  is  'that 
of  Srybsifc  and  Fir gov ^  '  who  found  24  surface  Redos  in  the 

Ic-ng  wavelength  ufu'eiardcd  limit,  of  which  12  dec%y  slowly  away 
free  the  surface  (like  Fuche-KHewer  aodss)  and  12  decay  rapidly 
within  s  few  lattice  spa-sings  -of  the  surface*  The  relationship 

a 

between  those  results  gsd  those  of  other  workers  has  not  been 
fully  established.. 

(vii.fi ,  Surface  Waves  in  the  Presence  of  Second  Sound 

Ssaonu  sound  -  x.e.,  propagating  te»psrature  waves  - 
hes  long  bash  known  in  liquid  holiua  ind  has  recently  boon  ob¬ 
served  is  solid  h&lU-sS  and  in ^  '  c  hn  investigation 

of  surface  effects  has  been  carried /out  by  Ioffe'  who  con- 

cidsred  tbs  cEouplicg  of  particle  waves  end  tsjaporature  waves 
through  the  thermal  expansion.  Ioffe  assumed  a  soni-infinite 
isotropic  el?. stic  continuum  and  obtained  two  surface  waves,  one- 
corresponding  to  the  Sayleigh  *?Ves  and  the  othor  having  a  very 
slow  velocity  given  by 

;  *  (crT>  *  < ps^/CT)  1/8st  (93 

where  .erift  the  thermal  expansion  coefficient,  p  is  the  density, 

T  is  tho  teapersturo,  C  is  the  specific  heat,  s^,  a^,  and  sT  are 
the  rpeads  of  bulk  longitudinal,  transverse,  and  second  sound 
waves,  and  CQ  is  a  constant  on  the  order  of  unity.  Ioffe's  slow 
wave  is  not  simply  a  surface  second  sound  wave,  but  a  coupled 
particls-teoperatur*  wave. 
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3L  SURFACE  SPECIFIC  HEAT 

Within  tl'e  harmonic  approximation,  tho  specific  hoat  o f  a 
crystal  can  bo  written  as  the  sum.  of  contributions  from  oach  of 
tho  s  normal  modaE  of  vibration  in  the  for® 


Cv-  kB  S  (J9ftu.'s)2  esch2^*^) 


(  94  ) 


where  6  »  (kgT)  “  and  wg  is  tho  frequency  of  the  a  th  normal  mode.  For 
a  bulk  crystal,  4he  specific  heat  bus  the  familiar  Debyo  form  at 
low  tomperatures.  When  surfaces  are  present,  the  normal  mode  fre~ 
quencies  are  modified  from  those  for  a  perfectly  periodic  crystal, 
and  the  specific  heat  acquires  a  surface  contribution  whose  form 
we  now  investigate. 


A.  Continuum  Theory. 

'»  • 

Early  treatments  of  the  surface  specific  beat  are  based  on  the 
Do bye  elastic  continuum  model  of  a  crystal.  Tho  Debye  theory  util¬ 
izes  a  method  of  counting  normal  modes  which,  in  its- original ‘form, 
yields  the  correct  result  in  tho  limit  of  an  infinite  crystal,  but 
an  Incorrect  result  for  ’a  finite  crystal,  lion  troll '  96  ^  improved 

the  Daoyo  theory  by  correcting  the  counting  procedure  for  the 

« 

normal  modes  and  found  that  the  vibrational  specific  hoat  should 
have  an  additional  term  proportional  to  the  surfaco  area.  In  tho 
low  temperature  limit,  this  correction  torn  is  proportional  to  T2 
in  contrast  to  tha  T3  dependence  of  the  bulk  term.  However,  Kontroll 
used  clamped  boundaries  rather  than  froo  boundaries. 

The  modifications  of  the  normal  modes  by  free  surfaces  havo 
been  taken  into  account  in  the  isotropic  elastic  continuum  theory 
by  Brager  and  Schuckowitzky ,  .  f  '  by  Stratton, 1  ■  and  by 


Dupuis,  Mazo,  and  Onsager. ^  ^  Wo  shall  utilizo  Stratton's 

approach  in  tho  discussion  to  follow. 

In  the  continuum  approximation  the  mean  vibrational  energy  can 
be  expressed  as  an  integral  over  the  normal  node  frequencies 


<E>  -/  g(u,)  €  (»)  d* 


(  95  ) 


where  g(w)  d®  is  tho  number  of  normal  node  frequencies  in  the  range  ui 
to  oj  ♦  db  and  ?(cb)  is  the  mean  energy  of  a  normal  mode  given  by 


i(w)  -  fta)tn(io)  +  $] 
n(ui)  «  Co,)tu/kBT-l]-1 


(  96  ) 

(  9?  ) 


The  specific  heat  is  then  obtained  from 


*  ‘  C  -  3<E>/3T 
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Our  basic  problem  is  to  calculate  the  frequency  distribution  func¬ 
tion  g(u>) . 

The  enumeration  of  the  normal  modes  Is  facilitated  if  we  con¬ 


sider  a  slab  of  material  bounded  by  parallel  pianos  of -area  4L  L 

*  y 

and  separated  by  a  distance  JiL  .  The  normal  modes  for  the  isotropic 
elastic  slab  were  discussed  by  Lamb. ^10°  '  They  may  be  classified 
as  either  symmetric  or  antisymmetric  with  respect  to  reflection  in 
tho ’ said-piano  of  tho  slab.  A  further  classification  can  be  achieved 
in  terms  of  the  values  of  the  decay  constants  q, ,qa  of  Sec.  2A(i). 
For  qltq2  both  real,  we  have  Rayleigh  waves;  for  q£  real  and 
imaginary,  we  have  mixed  waves;  for  both  imaginary,  we  havo 

pure  bulk  waves  with  a  longitudinal  component.  In  addition  to 
those  three  types  of  modes,  there  are  pure  transverse  bulk  waves. 


-  -'*'t ; 

ter 
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'  Tho  wave  vector  rsagnitudo  x.  the  attenuation  constants  q^.qg 
and  the  frequency  te  -  xc  are  related  by  the  equations 

1.  -  t  Cl+2K}q12  -  Kq22}/(UK)  (  99  ) 

u>2  -  X2Ct2(q12-q22)(l+2K)/a+K)  (100  > 

where  K  «*  ji/X  and  X,ji  ore  the  Lam6  constants.  The  wave 'vector  com¬ 
ponents  q «  xi  and  q  »  xa  are  specified  by  periodic  boundary  condi- 
•  y 


tior.e 

qx~  »j.n/Xix  ,  nx-  0,£l,i2 .  •  (101»  ) 

V"  n7n/,h7  '  ny"  0,±1j42 .  (101b  ) 


*  Tho  spec if ic at ion  of  the  Rayleigh  waves  is  completed  by  the  boundary 
condition 


tanh  xq2Lz 
'  tanh  xq^ 


4  CUX)  q+2K-(/K)a 


kX 


(1+2K+J2)2 


«-[gK(o)] 


il 


(  102  ) 


•where  a  m  q2^i  and  the- exponents  +1,-1  refer  to  the  symmetric  and 

antisymmetric  modes,  respectively.  The  number  of  Rayleigh  modes 

with  frequencies  loss  than  or  equal  to  »  is  given,  to  an  adequate 

« 

approximation,  by 

Gr(«>  -  Va>2Y02/4t7Ls,ct2  (  X03  ) 

where  V  is  the  voluao,  Y^2-  (l+2K-s2K}/(l+2K)  (1-or2)  ,  and  «0  is  the 
solution  of  'Eq.  (  103)  for  x  r*  •  , 

•  "a 

The  mixed  sodas  have  one  attenuation  constant  replaced  by  a 
pure  imaginary  quantity  q^*-  iq^,  *or  largo  x,  tba  solution  of 


Eq.  <  102  >  is 

*2  **  <nV2Lz>+  6a 
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(  104  > 


where  fiy"  tan  C-ig^a)]  *  n2“  0,2, 4, 6,...  for  the  symmetric  aodcs 

-  1,3,5 ,7. . .  fo 
so  dee  Gy(ui)  is  given  by 


and  n*  1,3, 5, 7...  for  the  antisymmetric  modes.  The  number  of  mixed 
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where 
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1  “  n/  (1+  11&  tan2  *)  6ji(tan  <P)*P 

»  A  6 


(  105  ) 
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tan  qj  •  «’*  i«,  .and  6y(tan  <?)  -  aCy/da'.  The  integral  I. can  be 
ovaluatod. by  contour  integration  yielding 


(1+2K)' 


(i-*0Ku+2xr+ff0} 
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C  107  ) 


The* bulk  nodes  have  both  attenuation  constants  replaced  by 
imaginary  Quantities,  *  iqj  and  q2  -  iq2-  For  large  *,  the  solu¬ 
tion  erf  Eq.  (  102  )  can  bo  written  as 


*2  "  <«VV-  *B 


(108) 


where  6g  «  tan  [gjr(a)tan(xq^  Lz)]  and  the  »a  are  non-negative 
integers.  For  tho  number  of  bulk  nodes  Gg(ui) ,  we  have 
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'  finally,  the  values  of  q£  for  the  pure  transverse  bulk  wavos 


are  specified  by 


q£  ”  nnz^ZLzr  n  “  °*1'2*3»4 <5* • ■ 


(  110  ) 


The  number  of '  transverse  modes  with  frequencies  less-  than  or  equal 


tO  u)  is 

„  /  3  2 

gt<«)  -  4r  f  -%  -  JIaL-T 

4  n“  l  leVt 


.  (Ill  ) 


The  total  number  of  modes  with  frequencies  less  than  os'  equal 
to  a  is  giv'm  by 


G(u»)  “  G^<a)  +  Gy(ui)  +  Gg<w)  +  G,j,(w) 


.(  112  ) 


The  part  of  G<«j)  involving  «r  gives  the  bulk  contribution  to  the 
'specific  heat,  while  the  part  involving  »  gives  the  surface  con¬ 
tribution.  Tife  total  number  of  modes  is  3N  whore  N  is  the  number 

of  atoms  in  the  crystal.  If  uiQ  is  the  maximum  frequency  of  the 
crystal,  theu 


3H  -  G<u>n)  . 


(  113  > 


Utilizing  the  above  expressions  for  the  contributions  to  G(e>)»  wo 
obtain  from  Sq.  (  113)  an  expression  for  s>*  in  the  fora 

SI 


“a  -  *nB  +  A®aS 


<  114  ) 


whore  is  the  maximum  frequency  of  the  bulk  crystal  given  by 

ct^8n-N/¥2*/3/|2+CK/a+2K)33/2  j  1/3  and  awmS  is  the  change  in  « 
produced  by  the  surface, 


Ao?«S  *  * 


n{Y0  -  I  ~  £K/4(l*3K)3l 
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where  A  is  the  area  of  the  crystal  and  n  ~  N/V. 

The  differential  frequency  distribution  gCiy>  is  related  to  G(uj) 
by  c(s/>  ~  dG(ui)/<lai.  Tae  internal  energy,  neglecting  the  zero-point 
contribution,  can  be  written  from  Eq. (  i55  )  as 
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froa  Sq^C  34  ),  the  specific  fcsat  becomes 
*V-  2  ftai/V 
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(  117  > 


Breaking  up  both  u>H  and  dG/djj  into  their  bulk  and  surface^  parts,  we 
can  write  the  surface  contribution  to  tho  apse if ic  heat,  Cvgf  in  the 


C  „  -  kAF(K)  n2/3  [D(x)-  E{x)i 


CU8  } 


where 


1/3  i~[K/(l+SK>3} 


{ 2+£S/ (1+2K) X  '  }  ' 
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and  «  «  haig.gAT  . 

In  tho  low  temperature  limit,  on®  finds  with  the  aid  of  the 
explicit  expressions  for  V  and  X  that 

O  •  * 
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^3  2e,-  3c  c  +3c 

cvs  -  3n«s>  fi  xV.;iiH'- 

b  ctV  4  t> 


AT 


(  120) 


where  c  is  the  speed  of  longitudinal  bull;  waves  specified  by 

■v 

c3«  (X+^i)/p0  and  C(3)  is  the  Itiemann  zeta  function  of  argument  3. 

Equation  (  120)  was  first  obtained  by  Dupuis,  Mazo,  and  Onsagar'  1 

who  used  a  contour  integral  representation  for  Cy, 

Tho  surface  specific  heat  is  seen  to  be  proportional  to  ife© 

•* 

surface  area  and  to  the  square  of  the  absolute  temperature  at  low 
temperatures.  For  fine  powders  at  sufficiently  low  temperatures, 
tho  surface  specific  heat  may  become  comparable  to  the  bulk  T3  coa- 
tribution. 


•  B.  'Discreto  Lattice  Theory. 

Tha  first  lattice  dynamical  theory  of  the  surface  specific  heat 

was  developed  by  Patterson^01  /  who  utilized  tho  Rosens tock-Kowoll^- 41  ^ 

model  involving  nearest-neighbor  central  and  non-central  interactions 

in  a  simple  c  ibic  lattice.  At  low  tersperaturos  kis  roC”tt  reduced 

_2 

-  to  a  quantity  proportional  to  AT'.  However,  the  lack  of  rotational 

•  invariance  of  tho  Rosenstock-Newell  model  makes  the  result  lack 
quantitative  significance. 

A  more  recent  lattlee-dynaaical  calculation  is  that  of  Uaradudin 
and  ffallis  who  treated  the  surface  as  a  perturbation  &nf  employed 
a  Green’s  function  method.  Tho  perturbation  involves  sefctiitg  to  zero 
the  force  constants  coupling  atoms  on  opposite  sides  of  the  plane 
defining  the  surface.  In  principle,  it  sight  also  involve  changes 
in  force  .constants  at  or  near  the  surface  but  not  crossing  iSs@  de¬ 
fining  plane.  Ike  latter  changes  s^glocted  is  tho  treatment  to 
follow,  *  ’  - 


i 


Jf4iuation<g4  )  tbe  specific  heat  nay  bo  ro-8xpress«d  as  a 


power  earios  in  ©xp(~3hm  ) 

u 

0.,  “  kP  S  bS  0^b)2  eap'(“.ipft»  )  .  {  121  ) 

v  B  n-1  s  s  * 


The  sun  over  the  normal  nodes  fiaa  bs  replaced  by  a  contour  Integral, 
so  that 
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(122  ) 


where 

F(z2)  v  S  [l/(z2n»2)3 

s 


(  123  > 


and  tbo  contour  encloses  the  positive  real  axis  and  thus  encloses 

2  • 
the  poles  at  -^g  of  F(z  ). 

fc 

The  frequencies  »_  are  those  of  the  crystal  with  a  pair  of  free 

a  * 

surfaces.  He  shall  also  need  the  frequencies  oj0g  of  the  perfect 

•  2 

periodic  crystal  and  the  function  FQ(z  }  obtained  by  replacing  «s 
by  (bCs  In  Eq. (123) o  The  surface  contribution  to  the  specific  heat 
can  then  be  expressed  as 
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z3e“R0ftz[ F(z2)~F0 (z2) j  az 


<  124  ) 


o  o 

By  expanding  ir.  inverse  powers  of  z,  one  finds  that  F(z  )-?o(z  )  is 
of  order  z  ,  so  that  the  contour  caa  be  taken  along  the  imaginary 


axis  and  closed  by  an  infinite  semicircle  lying  in  the  right-hand 
side  of  the  complex  plane.  Since  the  contribution  from  the  semi¬ 
circle  vanishes,  the  surface  specific  boat  can  bo  written  .as 
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(125) 


where  y  -  ~lz , 

v> 

/* 

J(n9fc)  -  f  .  y3sin<njifty}ft{/}<iy  , 
**o  ^ 

and  Q(y2)  -  y  (~y2)  -F(-y2)  . 
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(  126  ) 


(  127  ) 


In  the  low  teopcraturo  liailt,  ws  are  interested  in  the  asymptotic 
behavior  of  J (o$h)  as  n9ft  ■+  *  .  How  it  can  be  shown vX*^  that  fi(y2) 
has  as  its  only  singularity  a  logarithmic  dependence  on  Jy|  be 
7  -*  e  « 


fi(y  /  ~  -  0o  ln|y|  +  oClnjyJ)  , 


(  128  ) 


so  that  the  .dominant  terxa  ia  the  asymptotic  behavior  of  J(njft)  aa 
nph  “  is  give®  by  ^“°3^ 


-J(ngft)  ~  -  3«oo/Co9h>4^  0(0^)  . 


(129  ) 


Tha  low  temperature  limit  of  the  surface  specific  boat  then  becomes 


CvS  "  6v*0C(3)kB(k B1/hf  +  oCT2)  . 


(  130  ? 


The  problem  now  is  to  calculate  &(y  }  and  this  requires -  eon* 
sidoration  of  the  lattice  dynamics. 

hot  the  force  constant  matrices  for  tho  perfect  crystal  and 
for  tho  crystal  with  the  free  surface  bs  $0  and  respectively. 
The  equations  detcrjainins  tho  aorsal  sodas  for  these  two  cases  can 
oo  written  as 
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acid 


whore 


L  u  ~  o 


(L  ~.  l>h)  u  *>  o 


C.  131  ) 


<  132  ) 


h  ~  ® «  ~  io 


.  (  13?a) 

“  i  "  io  .  C  133b) 

M  is  the  (diagonal)  mass  matrix,  and  u  Is  the  complete  displacement 
vector  of  tho  atoms.  $L  is  the  perturbation  nat'rix  associated  with 
the  free  surface. 

Wo  now  introduce  the  Gretm’s  function  matrices  C  and  U  for  the 

*Sf  M 

perfect  and  perturbed  crystals,  respectively,  by 

« 

G  -  IT1  (  134  } 

<V  A*  . 

and  U  -  (L  -  6 L}”1  .  *<  135  ) 

*w  »w  «#  ( 

One  sees  from  Eqs. (134)  and  (1S5)  that  U  is  a  solution  of 
B  - G  +  G  6MJ 

M  M  r#  <V 

*G  +  G?G  C  13S  ) 

M  M  M  , 

where  the  matrix  T  satisfies 

‘  T  »  6L  +  4L  G  T  .  (137  ) 

w>*  n# 

2 

It  will  turn  out  that  n(y  )  can  bo  expressed  in  terms  of  T  and  that 
T  can  ha  aalculatod  exactly  for  certain  c-ses  of  interest, 

•  At  this  point,  we  find  it  convenient  to  express  the  elements 

of  if  in  terms  of  the  eigenvector  components  B, Ctaa)  and  eigenvalues 

2  ' 
aa  of  the  dynamical  matrix, 
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■shore  the  B_(4xa)  satisfy  the  normalization  condition 

}S 


E  B  Ctxtt)B  ,  Uxa)  ~  6SS.  .  (  139  ) 

4x0 

Utilizing  Eqs.  <123)  and  (138),  we  see  that  F(u>2)  can  be  written  as 

(  1^0  ) 


F(a2).  -  S  M  U  a(4K;4K;a-2)  . 

4xa  * 


2 

A  similar  expression  holds  for  ?0(ia  )  in  terns  of  G  -  namely, 


f0(b2)  ■£  m,  g  a<4x;4*;a,2)  , 

4xa  * 


<  141  ) 


so  that 
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,<o), 


The  eigenvectors  of  the  perfect  crystal  (txa).  can  be  taken  in 
the  fora 


b£°}Um)  -  .&<x|kj)ei^’5(4) 


(  143  ) 


whore  the  nods  index  s  now  stands  for  the  wave  vector  k  and  branch 
Index  The  elements  of  G  can  then  be  written  as 

2  1  „  '®n  4  j  k  j  }o  (k '  j  kj ) 

G^U.HU'k';*2)  -  - - — r  2 

aP  tsfv 


mK\,)s  kj 


x  e 


ixU)-x ({,')]  . 


(  14-4  ) 


The  matrix  T  is  conveniently  calculated  by  means  of  a  double 
Fourier  transform  defined  by  1 
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x  oiSl*S^>+iS2‘S^)  t<k1.t1;k2j2'<u!2)  .  (  145  ) 


Substitution  of  Sq. (14 5)  into  Eq«(137)  yields  the  integral  equation 
satisfied  by  t(kj  ;k'  j  *  jo?2) 

t(kj;k'3';»2)  -  V(teJ s**4' S«2)* 

M  M  fV 

V(kj  j-k. j,  ;<y2>  „ 

+  2  — - -  t<k3a,jk«j«j»  >  (146) 
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where 
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2  2 

In  tents  of  t(kl ;k’ j ’ ;»  } ,  the  function  0{y  )  can  be  expressed  as 


•  2'  ~  y  > 

fi(y)  "  -  S  * — x — o — — y 

kj  £y  +  »  <kj>32 


(  14a  } 


The  key  to  obtaining  ft(y  3  is  the  solution  of  tho  integral 
equation  (143).  la  general,  this  is  a  difficult  task,  but  it  can 
be  ftceoaplisfesd  for  special  casus,  ffe  new  discuss  one  such  case  - 
naasly,  a  sisspl®  cubic  lattice  with  nearest  and  next-nearest  neighbor 
central  forces.  •  The  potential  energy  for  this  social  esn  bo  written 
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where  the  prime  on  the  first  sum  means  tne  torn  for  4  -  i*  is 
omitted  -and  the  coupling  constants  t)_0(r  »)  may  be  taken  in  the 

Q-p  Cv 
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With  central  forces  Cauchy's  relation  c^2*°  c^^  must  be  satisfied  by 
the  elastic  constants.  We  also  impose  the  isotropy  condition 
Cj^"*  &12*  ^c44  wllicl1  l®ads  to  the  restriction 


9”  (aQ)  -  aD) 


(  151  ) 


where  &0  la  Che  lattice  constant.  In  the  long  wavelength  limit 

(which  is  all  we  need  for  the  low  temperature  specific  heat),  the 
2 

eigenvalues  ®  (kj)  are 
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(  152a) 


«2(k2)  -  a.'2(k3)  -  (c,./p)  i2  . 
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Too  eigenvector  o(kl)  is  given  by 
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The  other  two  eigenvectors  can  always  bs  eliminated  with  the  aid 
of  the  closure  relation 
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For  tho  model  just  discussed,  we  introduce  a  pair  of  (001) 
free  surfaces  by  setting  to  ZGro  all  interactions  coupling  atoms 


on  opposite  sSdos  of  the  plane  s  -  ia0*  A  total  of  five  "bonds" 
are  thus  cut  which  link  a  given  atom  in  the  plane  4  "  1  with  atoms 
in  the  plane  0.  The  perturbation  &L  is  specified  by 


6LaS<U,)"  <Pa9(rU,)[642064:t,l+  64a.064alj  *  1  *  l< 

6Lag ^^(tt*) 
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The  kernel  V(kj;k,j';u>  )  then  becomes 
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where  4  «  4-4,*  , 
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X  0(l/2)k-x(Z)  sin  J  k-x(t)  , 


(  15’  ) 


3 

and  the  volume  of  the  crystal  is  (La0>  ,  It  should  be  emphasized 
that  4  takes  ch  only  tho  five  values  corresponding  to  the  bonds 
broken  in  forming  tho  surface. 

Frca  JEq. (155)  one  sees  that  the  kernel  is  separable.  Equation 
(146)  can  then  be  solved  by  standard  nethoas  to  yield  the  result 
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t  ••V-v**—'.  «  »  ■«,  -  sf. 


wliero  I  is  the  5  x  5  unit  aatrlx  and  U(k  k  )  is  a  5  x  S  matrix 

**  x  y 

whose  elements  are 


o  v-(-kJ)v~, (kj) 
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It  is  convenient  to  introduce  a  new  5x5  aatrix  n(k  k  :u)2)  by 

^  x  y1  * 

&(ksV*2>  “  if^xV^^xV®2*3  \^txky  >  (  160  ) 

whore 

E«'(k*V  “  644'e<1  2)a°a^X+kyiy)  •  •  (  *61  ) 


Using  Eqs.(l57)>  (15S),  and  (ISO),  one  obtains  for  ft(y2)  the 


expression 


Q(y  >  -J£  S 


w-(kj)w",(ki) 


'Jj  -*3 — ~n-  ri~a(k  k  ;~y)] 

ui  kj  r?k+u,2(kj)32  w~~  *  y  u> 
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whore 


^~j)  ~7l X?  |k-x(D  , 


(  163  > 


Tho  singular  bebp-vlor  of  Q(y2)  as  y  i)  arises  froa  a  diver¬ 
gence  in  the  sum  over  k  is  Eq. (162)  associated  with  spall  k  values, 
in  fact,  it  is  thO  matrix  ;®2]  which  is  singular  as 

iS  **  °*  can  therefore  restrict  our  attention  to  snail  k  values, 
replace  the  bus  over  k  by  an  integral.  The  result  for  fj(y2}  is 
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whore  A  is  the  area  created  by  cutting  the  bonds, 
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and  c2  -  3c2  -  3c2  -  3*>»  <2*  a0)/aQp  . 


The  evaluation  of  the  elements  of  and  using  Eq.  (165)  la 

straightforward  but  tedious  and  is  given  in  detail  in  Ref. (102). 

2 

The  result  for  the  limit  of  fi(y  )  as  y  0  is 


n(y2)  -  -  (10/3)  ^njyj  +  o(tnjyj)  . 

SitcT 


(  166  ) 


Using  Eqs. (130)  and  (166)  we  r.btain  for  the  surface  specific  boat 


at  low  temperatures 


C  -  -  31T(kJ/h2)C(3X10/3cbAT2 


*  (  167  )  • 


teg  \  (g 

which  corresponds  to  the  results  of  Dupuis  ®t  alw/  and  of  Stratton 

*2  2 
when  c^  •*  3c  t  . 

At  this  point  it  might  be  well  to  discuss  experimental  attempts 

(1C3) 

to  observe  the  surface  specific  heat.  Lien  and  Phillips 
measured  the  specific  heat  of  SgO  powder  in  the  range  1.5  to  4aK  and 
found  a  contribution  of  the  form  C  •  0.163T2  ajcules/moie/deg.  The 
discrepancy  may  be  attributed  to  uncertainties  in  the  surface  area, 
the  Dcbyo  temperature,  and  the  true  elastic  constants  of  the  small 
particles.  Morrison  and  Patterson carried  out  measurements  or. 
KaCi  powder  in  the  range  10* K  to  room  temperature.  The  excess 
specific  ho&t  over  the  bulk  value  between  10° K  and  20* X  is  about 


1,5  tlxaos  tho  theoretical  surfaco  epocific  heat  of  Dupuis  ot  al'  , 

Similar  results  have  recently  boon  obtained  by  Barkraan,  Anderson, 

and  Brackott^®^  for  KaCl  powder.  However,  the  coefficient  multi- 
o 

plying  T  In  tho  surface  specific  heat  was  found  to  be  about  40  por 
cent  less  than  the  value  of  Uorrlson  and  Patterson  and  Is  therefore 
in  bettor  agreement  with  the  theoretical  result. 

Do  Sorbo  and  Nichols^0^  measured  the  specific  heat  of  lampblack, 
but  tho  excess  specific  heat  over  that  of  bulk  graphite-  seems  to  be 
duo  more  to  changes  in  the  elastic  constants than  to  surface 
effects. 

The  theoretical  discussion  so  far  presented  has  been  restricted' 

to  isotropic  materials.  Some  recent  work  is  applicable  to  non-iso*' 

tropic  crystals.  A  computer  calculation  of  the  surface  specific  beat 

for  a  (100)  surface  of  a  face-centered  cubic  lattice  with  Leonard- 

Jones  (6,12)  interactions  between  all  pairs  of  atoms  has  been  carried 

(  2SS 

out  by  Allen  and  Do  Wette'  '  using  Eq. (  94).  The  results  are  appli¬ 
cable  to  crystals  of  noble-gas  atoms  and  were  obtained  over  a  broad 
range  of  temperature.  At  low  temperatures,  the  expected  T2  dependence 
was  found;  at  higher  temperatures,  the  surface  specific  heat  was 
found  to  pass  through  an  inflection  point-  Indeed,  the  surface  specific 
heat  must  pass  through  a  maximum  and  approach  zero  as  the  temperature 
increases  since  tho  high  temperature  specific  heat  has  the  value  k0 
per  normal  mode,  independent  of  tho  presence  of  a  surface. 


.  This  behavior  has  been  found  by  Alien  and  De  who 

extended  the  calculations  Juat  discussed  to  the  (110)  and  (111) 
surfaces  of  noble  ga*  crystals  in  addition  to  the  (100)  surface. 


r 


•  -  ,  >-  /  ■.<•  1  .  W  '  .:.  t.  '  ,  k  *  *  * , 


ly.'* 
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Tho  surface  specific  heat  was  found  to  peak  in  the  vicinity  of. 
10 °X  for  Ne  and  at  slightly  higher  tempera turas  for  Ar,  Kr,  and 
Xe.  Allen  and  Do  Wette  also  calculated  the  surface  free  energy 
which  varies  linearly  with  temperature  at  high  T  and  approaches 
tho  zero-point  value  as  T  -•  0, 


Similar  behavior -has  boon  found  recently  by  Chen,  Alldrodgc, 

Do  Wetto,  and  Allot/109’’  who  have  calculated  surface  thermodynamic 
functions  for  MaCl  slabs  using  th3  rigid  ion  aoddl  involving 
Coulomb  interactions  and  noarest-neighbor  Bora-Mayor  interactions. 

The  results  for  the  surfaco  specific  heat  are  in  quite  good  agree¬ 
ment  with  experiment/105^  Tho  calculated  surface  specific  heat 
exhibits  a  maximum  at  45°K  which  is  close  to  the  experimental 
value  of  40°K.  % 

The  EosenstockrNewell  model/**1  ^  has  been  employed  by  Dobrzynski 
and  Leraan^110^  to  calculate  tho  surface  specific  heat  by  means  of  a 
phase  shift  method.  The  result  ia  the  low  temperature  limit  has 


the  fora 


CvS  -  (3/4rtc2 )  [1+2 (c44/cn)  ^C<3)  (k5/ft2)  AT2  (  *S8  > 


which  is  qualitatively  similar  to  the  result  of  Depuls  et  al_.  ' 

Cunningham^ 11 ^  has  calculated  tha  specific  heat  for  an 
anisotropic  (110)  surfaco  of  a  simple  cubic  lattice  using  two 
different  nearost-naighbor  forco. constant  models.  'One  model 
assumes  a  decoupling  of  the  motions  parallel  and  perpendicular  to 
the  surface;  the  other  is  the  Rosenstock-Xewoll  model.  The  Graen’a 
function  method ^1C2^  was  used  to  calculate  the  surface  specific 
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heat  which  was  found  to  bo  about  10  percent  lower  in  tho  isotropic 
»  • 

limit  than  the  correct  value.  This  discrepancy  may  bo  duo  to  the 
lack  .of  rotational  invariance  of  tho  models. 

The  surface  entropy  AS  of  fa-ije-centered,  body-centered  and 
diamond-type  cubic  crystals  has  been  investigated  by  Dobrzynski  and 
Friedel^**'^  utilizing  a  nearest-neighbor  contral-forco  model.  The 
(100),  (110),  and  (111)  surfaces  were  studied,  and  for  each  lattice 
it  was  found  thalj 

i 

*  ✓ 

AS{100)>  AS(110)>  AS(lll)  , 


C.  Adsorbed  Surface  Layers. 

Some  Investigations  have  recently  been  reported  concerning  the 
surface  specific  heat  of  crystals  having  adsorbed  atoms  on  the  sur¬ 
face.  Dobrzynski  and  Mills ^  have  used  a  Green's  function  method 
similar  to  that  described  above  to  study  tho  Rosenstock-Hewell  model 

of  a  simple  cubic  lattice  with  a  layer  of  adsorbed  atoms  on  a  (100) 

•  « 

surface.  The  central  and  non-central  force  constants  were  assumed 
to  be  equal,-  and  the  same  value  of  this  constant  was  used  in  the 
bulk  and  at  the  surface. 

Taking  the  masses  of  the  surface  and  bulk  atoms  to  be  M_  and  H, 

*> 

respectively,  Dobrzynski  and  Mills 
at  the  surface  gave  a  contribution 
temperatures  of  tho  fora 

ACvS  -  »s(i«  s/S’/5)n2a  k3 


at  tho  mass  deviation 
ipeciflc  heat  at- low 


(T/8jj) 
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whore  Us  is  the  number  of  surface  atoms  (mass  Mg),  a-  (Mg-M)/M, 
and  0O  is  the  bulk  Debye  temperature.  Since  this  contribution  is 
smaller  on  the  order  of  T/0d  than  the  ordinary  surface  specific 
heat,  it  will  bo  very  hard  to  observe  experimentally. 


Additional  calculations  of  surface  thermodynamic  functions  for 

n  13) 

adsorbed  particles  have  been  made  by  Allen,  Alldredge,  and  De  Wotts 
for  the  (100) ,  (110) ,  and  (111)  surfaces  of  a  face-centered  cubic 
lattice  w<th  Lennard-Jones  interactions.  The  changes  in  the  thorao-  . 
dynamic  'functions  due  to  the  change  in  mass  of  the  adsorbed  particles 
wore  calculated  numerically  as  functions  of  temperature. 


The  offset  of  a  step  on  the  surface  specific  heat,  entropy. 

and  fres  energy  has  been  wrxei  out  by  Jiasri,  Allan;  and  Dobrzynsfci 

for  a.ClQO]  stop  on  a  (001)  surface  of  a  simple  cubic  crystal  using 

< 

the  Eosenstock-Nawel  Model-  The  step  is  found  to  increase  the 
surf specific  heat  and  entropy,  but  decrease  the  free  energy. 


(114  > 
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'  4.  ATTENUATION  A NO  AMPLIFICATION  OF  SURFACE  WAVES 

Only  recently  have  either  the  experimental  or  theoretics! 

aspects  of  the  attenuation  of  surface  waves  received  detailed 

attention.  We  first  survey  the  theoretical  situation  and  then 

briefly  discuss  tho  available  experimental  results. 

A.  Theoretical  Aspects  of  Surface  Wave  Attenuation 

Various  mechanisms  can  contribute  to  the  damping  of  surface 

waves.  These  include  the  enharmonic  interaction  between  phonons, 

1 

tho  interaction,  of  the  surface  phonon  with  impurities  and  imper¬ 
fections  in  the  crystal,  3nd  the  interaction  of  surface  phonons 
with  conduction  electrons  in  metals  or  semiconductors. 

Since  long  wavelength  surface  waves  penetrate  deeply  Into 
*  the  material,  it  may  bo  expected  that  the-  damping  constant  for 
such  surface  waves  can  be  related  to  the  damping  constants  for 
bulk  waves.  Press  and  Hcaly^115^  have  considered  the  isotropic 
case  where  the  Rayleigh  wave  velocity  cR  is  specified  by 


„2  4  ,  2  .  ,  2  • . 


(l70) 


where  and  are  the  velocities  of  ccaprosslcnal  and  shear 
bulk  waves,  respectively.  Dissipation  can  be  incorporated  by 
lotting  each  velocity  become  complex: 


h"  vi~  ct{1  +  i5r 

Ct  -t  Yt  *  Ct(l  +  isJ  * 

*a  -  VB  •  °R(I  +  i6R>^. 


<1T1) 


If  one  takas  the  total  differential  ©?  5JQ»  (170)  and  selvas  for 


the  increoent  JR  in  tersas  of  and  c^,  one  f iada 


°a 


.  (172/ 


The  daapisg  constant  is  proportions!  to  She  iisaginpry  part  of  the 
wave  vector,  Writing 


V  £ 


Iff 


<173) 


and  using  Eqs,  (171),  wo  obtain  for  & 


a 


Ski 

2c 


Cl?  4) 


The  draping  constant  for  the  Rayleigh  waves  arR  can  new  be  written 
in  torus  of  the  danping  constants  #  and  tr+  for  cosprossional 

V 

and  shear  waves,  respectively,  as 


<17S  ) 


The  theoretical  prediction  contained  in  Eq.  (I??)  vre  found  by 
Press  and  Hoaly  to  be  in  good  agreessont  with  experimental  data 
for  plexiglass  sheets. 
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.  .The  procedure  of  adapting  the  theory  for  bulk  wavo  attenua¬ 

tion  to  surface  waves  has  been  pushed  further  by  He  Bride*  li6^  who 
.  applied  the  results  of  Maris  for  bulk  wavo  attenuation  by  cubic 

anhannonicity  to  the  surface  wave  problem.  Even  at  liquid  helium 
tooperatures,  the  wavelength  of  typical  thermal  phonor.  -»  is  such 
less  than  the  wavelength  of  the  ultrasonic  surface  phonon,  so  the 
principal  scattering  of  the  ultrasonic  surface  phonon  io  by  bulk 

phonons  and  not  b£  thermal  surface  phonons.  In  particular,  the 
.  ,  *  '  1 
dominant  process  is 

£  +  T  •*  T 

where  S  stands  for  the  ultrasonic  surface  phonon  and  T  stands  for 
a  thermal  transverse  bulk  phonon.  At  low  temperatures  whore 
ojt  »  1  (t  is  the  lifetime  of  the  T  phonons).  He  Bride  finds 


T4  t"1 


(  176  ) 


where  T  is  the  absolute  temperature.  This  result  is  pt  variance 
with  results  of  other  workers  to  be  discussed  later. 

A  somewhat  different  macroscopic  approach  has  been  followed 
by  Maris who  related  the  attenuation  of  ultrasonic  surface 
waves  to  the  bulk  viscosity,  heat  conduction,  and  thermal  expansion 
of  tho  crystal  in  the  limit  that  tfc^  mean  free  path  oi  the  theraal 
phonons  is  much  less  than  the  wave  length  of  tho  surface  wave,  The 
procedure  should  therefore  be  applicable  at  high  temperatures.  A 
similar  treatment  has  been  reported  by  King  and  Sheard^^ 
who  give  a  specific  expression  for  the  surface  wavo  damping  constant, 
•  The  displacement  vector  u  which  satisfies  tho  free  surface  boundary  ' 
conditions  can  bo  written  in  tho  fora  (cf.  Sq,  (8» 


1 

S 


& 


M 


c. 


-  S  Bj  Kj  exp{i(Hj*£  “  wt)J 

.J  ■—  X 


(  177  ) 


where  3.,  and  have  Cartesian  components  whose  analogues  in 

Eq.  (8)  are  ,  -iCj  and  xi,  xmf  ijcSj  i  respectively.  The 

viscosity  tonsor  components  ni1ki  enter  the  Bayleigh  dissipation 


function 


*  "  ^ijk-t  eij  ®kt 


<  178  ) 


whore  is  an  element  of  the  strain  tensor 


•ft*  %) 


(179  ) 


For  a  general  elastic  medium,  the  damping  constant  a  takes  the 


a  -  8.7  w2r,oJf/2pc|  dB/ca 


<180  ) 


where 


‘eif  ^  £ “  ”* ;;  *r -  ,  ( lei  ) 

j  j  »  Kj 1  B±  j  *  Q  j  »  '♦'ij) 

ij'£  “  ^  ®2X2  *  ®3X3^  53  x^=fixi»  End  &3  "  1*1  j-  For  on  isotropic 

medium,  there  are  only  two  independent,  non-zero  elements  of  the 
viscosity  tonsor  end  the  expression  for  the  damping  constant 
reduces  to  a  fora  equivalent  to  that  c?  Press  and  Eoaly. 

&  continuum  treatment  applicable  to  anisotropic  crystals  at 
low  temperatures  (®t  »  2)  has  been  presented  recently  by  King  and 
Shoard. *  20  X  Proper  account  ia  taken  of  the  free  surf see.  boundary 
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conditions.  The  displacement  voctor  u  which  satisfies  these 
boundary  conditions  can  be  written  in  the  form  given  by  Eq.  (177). 
Since  the  damping  constant  is  calculated  quantum  mechanically 
using  perturbation  theory,  it  is  convenient  to  quantize  the  dis¬ 
placement  field  by  writing 


u  -  £ 
~  J 


I  fcSy )  [*  *■’  v*p 

+  Sj\\+  exp(-ikj-r)] 


'{  182  ) 


where  and  bJ[+  are  the  annihilation  and  creation  operators  for 
the  mode  k,  cg  is  the  velocity  of  the  surface  wave,  S  is  the  surface 
area,  and  a  is  a  normalization  constant  which  can  be  calculated  from 
the  kinetic  energy  (half  the  total  energy)  and  is  given  by 


if  -  £  B  *3  K*K  /{qT*+q  ). 

J,  j*  ;i  J  J 


(  183  ) 


Sing  and  Sheard  neglect  the  effect  of  the  surface  on  the  bulk  codes. 
The  corresponding  quantization  condition  for  bulk  modes  of  wave 
vector  k  and  polarization  p  can  then  be  written  as 


/  \5  /  •**.■*■  f 

tt  -  S  - }  e.  U  «  ~  -  +  a.  0  ~  ~  \ 

~  k.p  \2pV».  J  \  / 


+  ik- r 


<  ^4  ) 


where  is  the  polarization  voctor  and  a^,  are  the  anniS na¬ 


tion  and  creation  operators. 


The  cubic  enharmonic  interaction  in  the  continuum  approx imat ion 


has  the  form 


H3  -  (1/3  1 
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whore  the  coefficients  arG  related  to  the  second  and  third 

order  elastic  constant*..  Considering  only  interactions  between  a 
surface  phonon  and  a  bulk  phonon  to  produce  another  bulk  phoi.cn, 

k  +  (£j,Pj)^S£  ^2,p2^  > 

Kins  and  Shcard  obtain  for  the  matrix  clement  of  Hg  for  this  process 

-  s(4?) l(irf  "(spif  '■ 
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a  [n^^+Dr 


(186  ) 


where  e, ,c2  are  the  bulk  mode  velocities,  ay,a^,  and  n&  are  ths 
mode  occupation  numbers,  and  ; 

'  w/3!)  JU  f  "uVi  <18?  > 

where  P  denotes  the  six  parnutations  of  the  pairs  lij) ,  (kt)„  (an). 
Cains  perturbation  theory,  the  transition  rate  and  thence  the 
relaxation  tiwe  Tg  for  surface  phonons  can  b®  calculated.  The 


89 


result  is 

A-  -  £ 

"S  plp2 


rrh  xJk,k„  „  _ 

1 1 - sV^ - Cn,-n,> 

JJS4 !7  p  cs  3  CjC2 
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£  - i - 

J  ia<Ij'i'klz~k22 
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where  rlj  and  tig  &i*e  the  mean  occupation  numbers  of  the  bulk  phonons, 
and  H>g,  la^,  tug  are  the  frequencies  of  surface  and  bulk  phonons, 
respectively.  It  should  be  noted  that  the  dominant  processes  atten¬ 
uating  surface  waves  are  energy  conserving  (cf.  He  Bride^*6^),  so 
no  relaxation  time  for  the  thermal  hulk  phonons  appears  in  £q,(188). 

The  integrals  in  Eq.  *.188)  may  bo  simplified  if  one  introduces 
the  ns*  variable  K  «  kg-jk^  anc*  uaos  polar  coordinates  kj  0,  p  for 
with  K  as  axis.  Only  small  K  are  important,  so  both  bulk  phonons 

Evaluating 

the  integral  over  k,  in  Eq. (1S8),  one  gets 
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can  be  taken  in  the  same  branch  t  **  k^  ,  v^v0 
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Xj“  hc^k./kgT  ,  ana  X  ~  K£/k  .  Xn  the  low  temperature  liait,- 
*  4^/15,  oo  ono  hae 
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The  evaluation  of  the  remaining  integrals  in  Eq.(189)  Is  difficult  f 

t 

and  has  been  carried  out  in  approximate  fashion  by  King  and  Sheard  f 

! 

far  surface  waves  propagating  in  the  X  direction  on  tho  Y  and  Z  j 

i 

surfaces  of  quartz.  The  results  for  the  'Y  surface  differ  by  less  i 

than  a  factor  of  two  from  the  experimental  data  of  Salzmann  et  al5121^  1 

A  lattice  dynamical  calculation  of  surface  wave  damping  for 
an  isotropic  cubic  crystal  has  been  carried  cut  by  Maradudin  and 
Mills using  a  Green's  function  method.  A  model  with  nearest 
and  next-nearest  neighbor  central  harmonic  interactions  and  nearest 
neighbor  cubic  anharaonic  interactions  was  employed.  Using  standard 
techniques  of  many-body  theory,  Maradudin  and  Mills  relate  the 
damping  constant  rc(u>) (proportional  to  the  reciprocal  of  the  lifetime) 
for  the  a-th  normal  mode  to  the  proper  self  energy  P  £«)  through 

,  J5 

the  equation 
» 

k  T  * 

Ajw)  +  i^(ui) - §~  Lim  PiVi€)  .  (  193  ) 
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whore  a  {•»>  is  the  frequsney  shift  due  to  the  anbarmonicitf  and  j 

S  { 

c  is  to  be  taken  equal  to  the  normal  mode  frequency  sjg.  * 

•  •  2 

In  terms  of  tho  Ore  .n’s  function  )  introduced 

in  our  discussion  c-f  surface  specific  heat,  P„  (w)  can  be  expressed 

is  lowest  order  of  perturbation  theory  as 
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wh ct0  tho  ladices  >;.K*  have  been  dropped,  M  is  the  mass  of  an  atom, 
^0rl°’2°'3^'^'3'“3^  *s  t*,e  cu^ic  enharmonic  coefficient  appearing  in 


el°'2a3  ~1^5&3 

the  cubic  inharmonic  Hamiltonian 


V  (1/6)  2  2  « 


i^3  £»  ‘w*  ( 195  ’ 


and  2n4kBT/h  with  4  an  integer.  At  this  point  Mnradudin  and 
MU1-  expand  in  ^arms  of  the  eigonfrequencies  wCO)  and  eigenvector! 


Sy. j  (t) ‘for  tho  perfect  periodic  crystal 


2  aW<ta)a(o>  <  ,  x 

£lJl  ^2j2  -1  1  -2J2 


U^~l'*l,~2^2'“2)  *  (  195a) 
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where  U«r>  is  given  b?  Eq.(143).  The  proper  self-energy 
becomes 
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where 
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X  1  £2^2  *^3^3^ 
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and 
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The  sura  over  n  in  Eq, (198)  can  be  evaluate u  if  use  is  made  of  the 


.  o 

spectral  representation  for  u(kj,k'j';ta  ), 


u(kj ,k'j’ 


a 

»2>  -/ 


u(kj ,k'J ' ;v) 


dv 


(■200  ) 


v  -  u» 


Using  the  theory  of  contour  integration,  one  finds 

„?  _  (v^lk-lt-Vj)  •  E*T  C8<'’l>*»<v2))/(v1.v2-l«t).  (  801  ) 


where  9(v)  -  4  coth(hv/2fc_T).  The  damping  constant  r  («)  now 

*>  a  5 

follows  from  Eos.  (193),  (198)  and  (202)  and  can  be  written  as 


rs(m) 


-T?-  2)  £  S  2  C„(k.J.)  c  (k_j  ' 

ft  %  -1~2^3  ~li52^3  ^iJ2^3 


(  202  ) 
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Turning  now  tr  the  case  where  the  mode  s  is  the  Rayleigh 
surface  mode,  we  consider  the  quantity  u(kj,k’j';v)  in  the  presence 
of  a  surface.  Since  the  wavelength  of  ultrasonic  surface  phonons 
is  long  compared  to  the  wavelength  of  thermal  surface  phonons  and 
the  penetration  distance  is  proportional  to  thS  wavelength ,  the 
interaction  of  the  ultrasonic  surface  phonon  with  thermal  hulk 
phonons  is  the  main  mechanism  for  damping.  If  the  effect  of  the 
surface  on  the  thermal  bulk  phonons  is  neglected,  Maradudint and 
Mills  show  that  the  spectral  density  u(kj ,k’ j * ;uj)  is  well  approximated 
by  that  for  the  perfect  crystal  without  the  surface,  u^Ocj  ;k*  j  ’  ,v) , 
given  by 


uo(ki  ,k’,1 '  ;©)**  A  (k+k' )6^y  [£[©+©(kj)]-6[©^(kj)]}/2©(kj).  (203) 

The  next  ingredient  required  is  Cg (kj)Cs|k'j ' )  for  the  Rayleigh 

•  2 
wave.  Now  from  the  definition  of  the  Green's  function  ), 

we  see  that 


8s<£")Bs(£’*> 
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2  2 
where  Res  fU  )}  stands  for  the  .residue  of  U  .Ul' lit  )  at 

©g*  op  —w  wp  ~A' 

the  simple  pole  ©'*©_.  From  Eqs. (138)  and  (196)  we  eoo  that 
s 


Res  [u(kJ,k'J’;08))  , 
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2 

The  double  Fourier-transformed  Green's  function  u(kj ,k * J 1 ;©  ).can 


be  written  as 


•A*-»  »■ 
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uCkj.k'j’ju,2)  -  uo<kj;k‘J1;lu2)  +  &u(kj;k'j';u)2) 


(  206  > 


2 

where  uo(kj  »£* J ’  ;u;  )  refers  to  the  perfect  crystal  and  is  giver,  by 


2  ^(k-tk'H.,, 


-jT(kJ) 


(  20?  ) 


2  ^ 
while  &u(kj,k!j! ;w  }  is  the  increment  in  u(kj ,k* j * ;,/)  produced  by 

the  surface  and  is  given  by 


2  t(kj  .k'j*  ;uj2) 

Au(kj,k'j*;/) - *-=— - ___ - 

'  '  [«  -v ^(kj)][U)2-lu2(k,j,)3  . 


(  208  ) 


with  t(kj,k'j';uj  )  specified  by  Eq.(158).  Since  uq  has  no  polo  at 

the  Rayloigh  wave  frequency,  only  &u  contributes  to  the  residue  in 

Eq.  (205).  Uaradudin  and  Mills  evalxiated  £u  for  the  nearest  and 

next-nearest  neighbor  central  force  model  of_a  simple  cubie  lattice 

employed  in  the  surface  specific  heat  calculation. 

The  final  ingredient  required  by  Eq. (208)  is  the  anharmonic 
> 

coefficient  V(kj  .k^ ,k2J2>.  Maradudin  and  Mills  argue  that  the 

nearest-neighbor  interactions  are  dominant  in  the  anharmonic  ‘•asms 

if  the  interaction  potential  qj(r)  drops  off  rapidly  with  increasing 
•  , 

r.  The  cubic  anharmonic  coefficient  in  general  involves  ^"’(r), 
cp"(r),  andtp‘(r).  The  third  derivative  sp”Hr)  generally  is  the 
most  important,  and  only  this  contribution  is  retained  by 
Maradudin  and  Mills.  Furthermore,  for  deeply  penetrating  Rayleigh 
waves,  the  effect  of  the  surface  on  the  anharmonic  coefficients 
gives  a  negligible  contribution  to  the  damping  constant. 

The  Rayloigh' wave  damping  constant  can  now  be  evaluated  for 
the  nearest  and  next-nearest  neighbor  central  force  model  with 


isotropy  in  the  long:  wavelength  limit,,  The  details  are  cumbersome 
and  can  be  found  in  Ref.(l22)*  Tho  result  in  the  low  temperature 
limit  is 


rR(»a)  -  C(4)  [0.0521  +  0.337  (  ^JW  *  *  <*0>]8  X 


U)p/kBT5 


~5~2~T~$ 
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whore  rQ  is  the  nearest  neighbor  distance  and  tho  Rayleigh  wave 

speed  is  given  by  Cp“  0.919c^.  Tho  result  contained  in  £q.  (209) 

4  i 123) 

has  the  same  dependence  found  by  Landau  and  Rumor v  for  tho 

damping  of  bulk  waves.  Maradudin  and  Mills  have  used  their  model 

to  calculate  the  ratio  of  the  damping  constants  for  Rayleigh  waves 

and  transverse' bulk  waves  and  find  that  the  Rayleigh  wave  damping 

is  much  larger  than  tho  bulk  wave  damping.  This  result  is  at  least 

in  parti  howover,  a  consoquoneo  of  tho  approximations  made  in  their 

treatment. 

*  ft 

Besides  anharmonicity ,  Rayleigh  waves  c<sn  bo  dampod  thx-ougb 
their  interaction  with  impurities  or  other  imperfections.  Steg 
and  Xlcmens^*2^  have  calculated  the  Rayleigh  \fuve  damping  constant 
for  impurity  scattering  using  tho  isotropic  continuum  model.  The 
surface  Wi.  .*e  displacements  wore  taken  to  have  tho  usual  fora  given 
by  Eqs. (id  ),  (15  ),  end  (16  )f  hut  the  bulk  wave  displacements 
were  not  chosen  so  as  to  £ul*y  satisfy  the  free  boundary  conditions. 
Toe  .result  for  the  inverse  scattering  time  when  tho  impurity  ie  a 
mass  defect  at  a  distance  z  from  tho  surface  can  bo  written  as 


90. 

2 

—  *  ~L2AY~4  u;5  [A1(h2>+  A2(kz)+  A3(Ha)]  (  2-10  ) 

1 R  8p  8  Cg 

whore  n  is  the  concentration  of  defects  per  unit  surface  area, 

Am  is  the  additional  mass  of  the  defect,  and  A^Ag.Ag,  and  B  are 
dimensionless  parameters.  The  terns  A^.Ag,  and  describe  the 
scattering  of  the  Rayleigh  wave  lsito  other  Rayleigh  waves,  into 
longitudinal  bulk  waves,  and  into  transverse  bulk  waves,  respectively, 
while  B  is  a  slowly  varying  function  *>t  Poisson's  ratio. 

Steg  and  Xlcmens . investigate!  kho  ratio  of  scattering  into  other 
Rayleigh  waves  to  the  scattering  into  bulk  waves.  They  found  that 
this  ratio  is  significantly  greater  than  unity  for  wavelengths 
greater  than  the  distance  z  of  the  impurity  from  the  surface,  hut 
becomes  very  small  when  the  wavelength  is  ouch  les3  than  z.  Equa- 

c  • 

tio’n  (  21(S  predicts  a  very  strong  ^  .dependence  on  frequency  for 
z-0  (defect  on  tho  surface).  This  suggests  that  the  successful 
propagation  of  very  high  frequency  surface  waves  (a>  a  10GHz)  may 
require  very  smooth  surfaces, 

B.  Experimental  Aspects  of  Surface  Wave  Attenuation 

An  excellent  review  of  experimental  techniques  used  in  the 
study  of  surface  elastic  vavos  has  recently  been  published  by 
Sransfeid  and  Salztsann.  A  somewhat  older  survey  is  contained 

in  Viktorov's  book. We  shall  give  only  a  brief  discussion. 

A  number  of  methods  have  been  employed  for  the  generation  of 
Rayleigh  waves.  A  mechanical  method  utilizes  a  wedge  or  prism  placed 
on  the  surface  to  bo  studied.  A  transducer  placed  on  an  oblique 
surface  of  t n+.  wedge  launches  longitudinal  elastic  waves  of  speed 


97. 


Cy  into  tho  wodg.*.  At  tho  interface  between  the  wedge  and  the 
solid  of  interost,  the  longitudinal  waves  are  partially  converted 
into  surface'  waves  of  speed  C  which  then  propagate  along  tho  surface 
of  tho  solid.  If  0  is  the  angle  between  the  oblique  surface  of  the 
wedge  and  the  surface  of  tho  solid,  then  sing  -  C^/C.  The  prise 
oust  therefore  bo  made  of  a  material  with  a  small  longitudinal  wave 
speed,  <  C,  such  as  a  plastic. 

Another  probedure  applicable  to  piezoelectric  materials  involves 
a  scries  of  linear  electrodes  deposited  with  a  spacing  d  on  tho 
surface?27!*  the  electrodes  are  all  connected  together,  one  has  a 
single-phase  array;  if  alternate  electrodes  are  connected  togother, 
one  has  an  alternate-phase  array.  Application  of  a  rf  electric 
field  to  the  electrodes  generates  surface  waves  through  the  electro¬ 
mechanical  coup'ling  of  the  piezoelectric.  The  wavelength  of  the 
Hayleigh  waves  generated  is  equal  to  the  electrode  spacing  d.  Non- 
piezoelectric  materials  can  be  handled  by  depositing. a  layer  of 
piezoelectric  material  on  part  of  the  surface  to  be  studied  and 
then  depositing  the  electrodes  on  the  piezoelectric  layer. 

The  single-phase  or  alternate-phase  electrode  schemes  can  also 
bo  employed  to  detect  surface  *ivcs.  Recently,  optical  techniques 
using  lasers  have  been  develo,  *>-<«'• ,28^  for  detecting  surface  waves. 

Relatively  few  experimental  results  are  available  on  the  damp¬ 
ing  of  surface  waves.  Salzmann,  Plieninger,  and  Dransfeld^121^ 
have  studied  Rayleigh  wave  damping  as  a  function  of  frequency  and 
temperature  for  several  surfaces  cf  quartz.  At  low  temperatures 

A 

(below  APS),  the  damping  was  found  to  exhibit  the  e>T  dependence 
calculated  by  Maradudin  and  Mills. At  high  temperatures 
(above  60°  K)  where  tar  «  1,  the  damping  is  independent  of  T  anti 


i 


Sm&s;  * 


i 
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varies  as  w  .  Tho  iu  dependence  can  be  understood  in  terras  of  the 

viscosity  coefficient  theory  of  JJaris^ 118  ^  and  of  Kin**  and  Sheard^  11S^ 

The  actual  experimental  data  lie  somewhat  above  the  values  predicted 

(  1 2°  1 

by  these  theories,  but  improved  experiments  by  Budroau  and  Carr'  ' 

have  removed  this  discrepancy  as  well  as  the  discrepancy  with  tho 
froquoncy-independent  results  of  Daniel  and  Do  Klerk. ^130^ 


C.  Amplification  of  Surface  Haves 

Recent  interest  in  surface  wave  devices  has  focused  attention 
or.  the  amplification  of  surface  waves,  particularly  in  piezoelectric 
materials;  For  6uch  materials  it  is  possible  to  couple  drifting 
current  carriers  to  the  surface  wave  and  thus  transfer  energy  from 
tho  carriers  to  tho  surface  wave. 

In  practice,  one  can  avetd  screening  effects  in  the  piezo¬ 
electric  by  having  tho  carriers  in  a  thin  zlzb  of.  semiconductor  - 
placed  adjacent  to  the  piezoelectric.  This  case  has  beon  discussed 
theoretically  by  Gulyaev  and  Pustovoitf  131  ^  The  attenuation  or 
amplification  of  the  surface  wave  is  determined  by  tho  imaginary 
part  of  the  wave  vector  k, 
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where 
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F(u; , k)*~l ,  h  is  tho  thickness  cf  tho  slab,  Q  is  the  coupling  constant, 
rD  is  the  Debye  radius,  v  is  tho  carrier  collision  frequency,  <u 
is  the  frequency,  and  £2  are  the  dielectric  constants  of  semi¬ 
conductor  and  piezoelectric,  v T  and  ore  the  thermal  and  drift 

volocitios  of  tho  carriers  (assumed  non-dogonornto;  if  aogonorote, 
vT~v?/«/3  whore  vy  is  the  Fermi  velocity)  ,  and  cg  is  tho  velocity 
of  tho  surface  wtfvo.  From  £qs.  •(  211  ),  we  see  that  amplification 

•i 

occurs  if  v^  >  Cg. 

A  theorotical  discussion  of  amplification  of  Rayleigh  waves  by 

S 

an  electron  beam  passing  near  the  surface  has  been  given  by  Gulyaev 
and  2il'boraan$  ^32  ^  Viktorov^  ^  and  Kaliski^  ^  have 

discussed  the  particular  case  of  CdS,  Tho  effect  of  a  magnetic 
field  has  been  treated  by  Bers  and  Burke.  w‘  ' 

Tho  experimental  observation  of  surface  wave  amplification 

was  first  reported  by  White  and  Voltaer^  135  ^  and  by  Vas'kova  and 

(  137  ) 

Viktorov'  for  CdS.  The  procedure  involving  a  semiconductor 

slab  adjacent  to  the  piezoelectric  has  been  utilized  recently  by 

(  13  8  y 

Collins,  Dakin,  Quate,  and  Shaw'  '  to  study  amplification  of 
surface  waves  on  LlKbQg.  Tho  amplification  of  Blousteln-Gulyaev 
waves  propagating  in  the  [110]  direction  on  tho  (110)  surface  and  of 
Rayleigh  waves  propagating  in  the  [110]  diroction  on  the  (001) 
surfaco  has  just  boon  reported  by  Ludvik  and  Quata? 
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r.i AN  SQUAPE  D I E PLA C Ei'.E NTS  AND  VELOCITIES  0?  SURTACK  ATOMS 
Introduction 


A. 

The  mean  squnro  displacements  of  surface  atoms  can  be  studied 
oxporimontally  using  low-onorgy  oloctron  diffraction  (LEED)  or  tho 

Mbssbauar  effect.  Since  only  a  few  radioactive  nuclei  such  as 
57  119 

Fo  and  Sn  arc  suitable  for  MtSssbauer  investigations,  wo  shall 
devote  our  primary  attention  to  LEED. 

Let  us  consider  tho.  diffraction  of  lew-energy  electrons  from  a 
crystal.  As  the  temperature  of  tho  crystal  increasej,  the  increas¬ 
ing  vibrational  amplitudes  of  tho  atoms  cause  tho  diffraction  spots 
to  broaden  in  width  and  docrcaso  in  poak  intensity.  The  decrease  in 
intensity  of  a  Bragg  peak  is  described  by  the  Dobye-Waller  factor, 
exp  [-2M (lx)],  where 

«Ux>  -  £<[Q  *  u(ix)J2>,  .  (212  ) 


u(vx)  is  tho  displacement  of  atom  lx,  and  §  i's  a  wave  vector  to  be 
specified.  In  tho  kinematic  or  single  scattoring  approximation 
tho  intensity  of  scattered  electrons  from  a  crystal  with  one  atom 
par  unit  coll  can  be  expressed  ^ 

2  _  _  iQ-[R<0-*U‘)3 
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°  °  l  4’  ~  ~ 


X  0 


4<(S-Ch(£)  -H^‘>3)2> 


(213  ) 


i 


101. 


whore  I  is  tho  incident  intensity,  a  is  the  Thomson  factor,  f  is 
<-•  w  o 

tho  atoalc  scattering  factor,  or(4)  is  a  transnission  factor,  and 
Q  “  K*  **  X  with  X  and  X'  the  wove  vectors  of  the  incident  and  scat- 

^  4**  f*’  « 

torod  oloctron.  Tho  angular  brackets  in  Eqs.  (  212  )  and  (  213  ) 
denote  a  thermal  average.  The  factor  involving  tho  thermal  average 

in  Etj,  (  213  )  can  be  rewritten  as 
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Tho  quantity  I1  is  tho  Law  scattering  and  I2  is  ths  first-order 
thermal  diffuse  scattering. 
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From  £q.  (  215  )  we  see  that  the  moan  square  displacement  ol 


an  atom  in  the  crystal  can  bo  extracted  from  the  logarithm  of  the 


peak  intensity  of  a  Laue  scattoring  spot.  If  very  low  energy  electrons 


(typically  less  than  SO  volts)  are  used,  the  scattering  is  mainly  from 


the  surface  layer,  and  the  mean  square  displacement  then  refers  to 


a  surface  atom.  It  is  convenient  to  discuss  the  mean  square  dis¬ 


placement  in  torms  of  an  effective  Debye  temperature  S  which  is 


dofinod  by  the  rolation 


<u2>  -  <3h2T/tnkBe2)§(e/T) 


(  218  ) 


where 


4(x)  -  D1(x)  +  \x 


(  219) 


and  D^Cx)  is  the  Debye  function,. 

D> <5C)  ~h  f  ^ 

K  °  -1 


(  220  ) 


At  high  temperatures,  §(8/T)  *•  1,  and  <u^>  is  proportional  to  the 


tomporature  T. 


The  discussion  Just  given  has  been  based  on  the  kinematic  or 


single  scattering  approximation.  Since  very  low  energy  electrons 


aro  strongly  scattered  by  tho  crystal,  one  must  consider  the  effects 
of  multiple  scattering.  This  has  boon  done  by  Duke  and  Daranore^*^  ^ 


who  show  that  the  lattice  vibrations  enter  tho  thoory  in  tho  fora 


of  a  Debyo-Waller  type  factor  multiplying  tho  electron-ion  core 


scattering  amplitude.  The  multiple  scattering  loads  to  products  of 
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Dobye-Waller  type  factors,  so  the  peak  intensity  may  decay  more 
rapidly  with  increasing  temperature  than  would  occur  for  single 
scattering.  The  extraction  of  the  moan  square  displacements  of 
surface  atoms  f^pm  LEED  data  is  a  complicated  problem  when  multiple 
scattering  is  significant.  ’  . 


An  attempt  to  assess  the  importance  of  multiple  scattering 

•  t 

effects  on  the  extraction  of  surface  Debye  temperatures  from  LKED 
data  has  very  recontly  been  initiated  by  Jepsen,  Marcus,  and 
Jons.^142^  They  consider  the  (HI)  surface  of  silver  and  calcu¬ 
late  the  LEED  spectra  for  a  givon  value  of  9  from  dynamical  theory 
and  then  evaluate  the  results  using  kinematic  theory  to  re-extract 
the  value  of  0.  The  procedure  yields  values  close  to  the  surface 
Debye  temperature  for  electron  energies  around  40  -  50  eV  and 
values  close  to  the  bulk  Debye  temperature  for  energies  abovo 
140  eV.  Below  40  eV,  the  values  of  e.turn  upward  and  become 
larger  than  the  surface  Debye  temperature  by  10  to  15  percent. 


We  now  discuss  a  numbor  of  oxporimcptal  results  which  have  boon 
analyzed  using  kinematic  theory.  The  surface  Debye  tempera tares  so 
obtained  must  bo  treated  with  some  reservation  until  the  role  of 
multiple  scattering  effects  is  fully  clarified. 


B.  Experimental  Results 


Experimental  observations  of  surface  atom  mean  square  dis¬ 
placements  wore  first  made  by  Zamsha  and  Kalashnikov ^  on 
silver.  Subsequent  work  was  carried  out  by  Menzel-Kopp  and  Monzel^ 
on  copper  surfaces  and  by  MacRae  and  Germer^14£*  ^  on  nickel  surfaces. 
These  early  investigations  demonstrated  that  surface  atoms  typically 
havo  larger  moan  square  displacements  than  bulk  atoms.  The  effective 
Debye  temperature  for  surface  atoms  is  therefore  smaller  than  that 
for  bulk  atoms. 

The  availability  of  commercial  LEED  apparatus  during  the  last 
decade  has  stimulated  considerable  work  or,  surface  Dcbyo-Kaller  . 
factors.  Of  particular  interest  is  the  investigation  of  various 
surfaces  of  nickel  by  MacRae^1**6  ^  which  rovoaled  that  the  moan 
square  displacement  of  a  surface  atom  in  a  cubic  crystal  is  aniso¬ 
tropic  even  though  symmetry  requires  that  tho  bulk  atom- mean  square 
displacement  be  isotropic.  UacRae's  data  a-e  shown  in  Fig.  12 
for  the  (110)  surface,,  The  effective  Debye  temperatures  were  found 
to  bo  220°K  for  the  [110]  direction  (perpendicular  to  tho  surface) 
and  330°K  and  220°K  in  the  [110]  and  [001]  directions,  respectively 
(both  parallel  to  the  surface) .  All  three  surface  Debye  temper¬ 
atures  are  iess  than  the  bulk  value  of  3S0°K. 

Additional  experimental  results  are  available  for  other  f.c.c. 
metals  at  room  temperature  and  above.  Jone?,  McKinney,  and  ffebb^147^ 
studied  the  (111)  surface  of  silver  and  found  approximate  isotropy 
of  the  surface  mean  square  displacements.  This  may  be  compared 
with  tho  early  results  of  Zamsha  and  Kalashnikov^*43  ^  for  tho  (001) 
surface  who  found  a  strong  anisotrdpy  with  the  tangential  aoaa 


square  displacement  larger  than  the  normal  component.  Additional 
investigations of  silver  have  been  reported  by  Corotte, 
Ducros,  and  Macall  for  the  (111)  surface  and  by  Rovida,  Torrini, 
and  Zanazzl  for  the  (001)  and  (111)  surfaces.  The  latter  workers 
found  a  parallel  wean  square  displacement  at  the  (111)  surface 
which  is  nearly  equal  to  the  perpendicular  mean  square  displace¬ 
ment.  Normal  components  of  surface  mean  square  displacements 
have  been  measured  for  various  surfaces  of  platinum  by  Lyon  pnd 
Sossorjai'  and  for  palladium  and  lead  by  Goodman,  Farrell  and 
Somorj&i .  Some  additional  work  on  nickel  has  been  reported 

by  Andorsson  and  Kasemo. Copper  has  been  investigated  by 
Woodruff  and  Seah^152^  and  by  Reid^153^  and  iridium  by  Goodman. 
Both  Woodruff  and  Seah  and  Reid  discuss  difficulties  in  obtaining 
surface  mean  square  displacements  from  LEED  data.  A  summary  of 
surface  Debye  temperatures  for  fee  metals  is  given  in  Table  1. 
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Table  1.  Surface  Debye  temperatures  (normal  component 

of  displacement)  for  fee  metals. 


Kotnl 

Surface 

SS(°K) 

es(°K) 

Reference 

Hi 

(110) 

220 

390 

146 

M 

(001) 

280 

415 

i.51 

Ag 

(111) 

155 

225 

147. 

Ag 

(111) 

145 

225 

148 

Cu 

•(111) 

244 

334 

152 

cu 

(001) 

230 

334 

153 

Pb 

(111) 

44 

SO 

150 

Pb 

(110) 

25 

00 

154 

Pt  ‘ 

(001) 

118 

234 

.149 

Pt 

(110) 

107 

234 

149  ‘ 

Pt 

(111) 

111 

234 

149 

Pd'  . 

(001) 

181 

274 

,150 

Pd 

(111) 

145 

274 

150 

Ir 

•  (001) 

218 

285 

154' 

i.  =  : 

:^x:) . 
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VK?." 
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Somewhat  less  exoerimontal  data  is  available  on  bcc  metals 


Aldag  and  Stern^10^  and  Baudoing^^have  studied  the  (110)  su  .ce 
of  tungsten,  while  Estrup^57^  has  studied  the  (001)  surface.  Tno 
(001)  surfaco  has  been  investigated  for  niobium  by  Tabor  and 
Wilson^*5®  ^  and  for  chromium  and  molybdenum  by  Tabor,  Wilson,  and 
Bastow.^**®  ^  The  foregoing  work  is  restricted  to  room  temperature 
and  abovo.  Kaplan  and  Somorjai^6®  ^  have  made  the  first  investiga¬ 
tion  of  surface  Debye-Valler  factors  at  lo*?  tempera turos.  The  case 
considered  was  the  (110)  surfaco  of  chromium  at  temperatures  down 
to  100°K.  More  detailed  discussion  of  this  work  will  bo  ^ivsn  later. 
A  summary  of  surface  Debye  tempera turos  for  bcc  metals  is  given  in 
Table  2. 


Table  2.  Surfaco  Debye  temperatures  (normal  component 


• 

of  displacement) 

for  bcc 

metals. 

Metal 

Surface 

es(°30 

eB(°K) 

Reference 

W 

(110) 

200 

280 

155,156 

W 

(001) 

183 

280 

157 

Hb 

(001) 

106. 

231 

158 

Cr 

(001) 

175 

*  440 

159  . 

Cr 

(110) 

333 

450 

160 

Mo 

(001) 

239 

380' 

159 

Relatively  little  experimental  work  on  surface  Debye-Waller 
factors  has  been  carried  out  for  insulating  crystals.  However,  an 
interesting  investigation  has  recently,  boon  reported  on  solid  xenon 
by  Ignatjevs  at  al.^^^  who  found  that  tho  low  energy  electron 
scattering  from’  this  material  is  essentially  kinematic.  The  experi¬ 
mental  values  of  the  surface  Debye-Waller  factor  agree  quite  well 

t/ 


with  theoretical  values  calculated  for  fee  crystals  by  Wallis 

i  JL 62  )  (  04  s 

ot  al.'  “  "and  by  Allen  a«l  Getfobte  -  *  .  Semiconductors  have 

faicn  investigated  by  Nesterenko,  Borodkin,  and  Snitko-'  ^ 

who  found  surfaca-to-bulk  aeon  square  displacement  ratios  (normal 
components)  of  1.5  for  GaAs  (110)  and  2*2  for  silicon 

The  setjiuetsl  graphite  has  been  investigated  by  Alfcinet, 
Siberian,  end  Bien?ait^1S‘'^ who  obtained  a  surface  Dobye  temperature 
of  690°K  for  the  XGO.l)  surface.,  This  is  to  bo  compared  with  tho 
bulk  value  of  800°K.  A  model  calculation  gave  a  value  for  8g  6x 
720-?30°K. 

An  alternative  way  of  obtaining  information  about  the  mean 

square  displacements  of  surface  aloss  is  through  measurements  of  , 

the  recoilless  fraction  in  tfce>  HBssbauei*  effect.  The  recoilless 

fraction  is  .proportional  to  the  Dobyo-Wallor  factor  and  is  conso- 

quontly  temperature  dependent,  Ordinarily,  the  radioactive  nucleus 

is  an  impurity,  but  for  a  fe&  nctsrials  such  ss  iron  and  tin,  one 

can  got  inforaation  about  tho  host  crystal  atoms.  Results  have 

bean  reported  for  57?d  on  various  surfaces.  Ailen^3'65  ^  found  a 
•  57 

larger  teeaa  square  displacement  for  Fe  at  the  C  <00 >  and  (111) 
surfaces  silicon  than  in  the  bulk.  Qualitatively  similar  re¬ 
sults  were  obtained  by  Burton  and  Godwin^168  ^  for  57Fe  on  poly- 


crystalline  tungsten.  The  latter  workers,  however,  found  the 
opposite  effect  for  Fa  on  the  (100)  surface  of  silver.  This  re¬ 
sult  is  in  interesting  contrast  to  the  L5SED  results  of  Zsmsha  and 
Kslushnikov^ 143  ^  an#  .Jones  et  al.^i47  *  for  gbre  silver  surfaces. 


57 

tungsloii.  A  posslblo  explanation  is  that  the  Ko  a  tons  lio  abovo 

t  1 

tho  Bifriac©  layor  rathor  than  in  it,  but  furthor  experimental  and 

theoretical  work  is  required  to  clarify  this  situation.  It  is  clear 
that  a  major  problem  in  using  the  IHJssbauer  effect  to  study  surfaces 
is  tho  lack  of  precise  information  as  to  the  location  of  the  radio- 
active  nucleus  relative  to  the  surface  layor. 

C.  Theoretical  Results 

Qualitatively,  one  can  visualize  several  effects  arising  from 
a  free  surface.  The  Kean  square  displacement  of  an  atom  is  deter¬ 
mined  by  the  interatomic  forces  acting  upon  it.  A  surface  atom 
is  acted  upon  by  fever  neighbors  than  an  interior  atom,  so  this 
in  itsolf  will  generally  cause  tho  mean  square  displacement  to  be 
larger  for  the  surface  atoi^68i‘urthoraoro ,  the  interatomic  forces 
aay  have  different  coupling  constants  at  the  surface  than  in  the 
bulk.  Tfse  environment  of  a  surface  atom  has  different  symmetry 

f* 

than  that  of  an  Anterior  atom,  and  this  may  load  vo  an  anisotropy 
of  tho  moan  ■square  displacements  at  the  surface  evon  in  cubic 
crystals. 

To  calculate  the  moan  square  displacements  we  assume  initially 

that  displacement  components  satisfy  harmonic  equations  of 

^  • 

me  lion 

~  2  ^(txit'xOu  U'x*)  £221) 

*  *  ~  i’x’p  P  ~  % 

where  u  (tx)  is  tho  or-th  Cartesian  component  of  displacement  of  the 

Cf  *** 

x-tfe  atom  is  the  i-th  unit  cell,  a  is  the  nass  of  tbo  r-t'tx  kind  of 

X 

*soa  and  tbs  $  are  the  harmonic  coupling  coefficients. 

-  Cfji  'fi 

Asking  the  transformation  •  - 


where  the  D  Alx.U'x,')  are  the  elements  of  the  dynamical  matrix 
erg  ~  ~ 

dofined  by 

DtfS(£x;£,x,)  “  *aSiC£x:£'x,)  ■  (224) 


Tho  presence  of  the  free  surface  is  reflected  in  the  values  of 
appropriate  elements  of  either  the  coupling-constant  matrix  or  the 
dynamical  matrix. 

The  dynamical  matrix  can  be  diagonalized  by  means  of  the  nor¬ 


mal  coordinate  transformation 


vJ£x'>  ~  S  V(~x)Qp 


<225  ) 


where  e  (£x)  is  f-xor-th  component  of  the  eigenvector  of  the  dynacr- 

PQf  ** 

ical  matrix  for  the  p-th  normal  code  and  Q  is  the  normal  coordinato 

*  >.  * 

for  that  mode.  Then  *  * 


<u2Ux»  -  <l/«x)  T)  !opa<ixj-  ( |Qp|) 
where  wo  have  used  the  fact  that 

<QpV>  “  <iQPl2>  W  •• 


<  22s  :i 


(  227  ) 


How  from  tho  quantum  statistical  mechanics  of  a  harmonic  oscillator. 


<|QP!2>  -  €(«9/u,i 


(228  > 


wbero  u>  is  the  frequency  of  normal  node  p  and  6(®_)  is  tho  mean 
P 

energy  of  the  mode  given  by 


C(aO  “  $  ha  coth(ha  /2kRT) 
P  P  " 


P  p  P  _ 

The  mean  square  displacement  component  can  now  be  written  as 


(229  > 


(230  ) 


( 140}  • 

Alternatively,  using  a  well-known  theorem  of  matrices'  '  and  Eq. 
(229  >  wo  can  rewrite  Eq.  <230  )  directly  in  terms  of  the  dynamical 


matrix 


-  <u*Ux)>  -  (ft/2mx) [D~Scoth(ftD*/2k3T) • 

At  high  and  low  temperatures,  we  have  the  limiting  results 

<**<w>  -  .  t  >  e0 

tufa»  .  (X^.y)lU-ihtta^  .  T  -  0°K  . 


C  231  ) 


(232a) 


(232b  ) 


For  a  monatomic  crystal,  tbo  mean  square  displacement  is  independent 

_1  < 

of  the  atomic  mas3  at  high  temperatures  and  varies  as  m  a  at  low 
temperatures. 

In  actual  calculations,  it  is  convenient  to  introduce  periodic 
boundary  conditions  in  "the  two  directions  (say,  the  1  and  2  direc-. 
tione)  parallel  to  the  surface  and  an  associated  two-dimensional 
wave  wrier  q  **  (q.  ,q_> .  The  eigenvectors  then  have  the  form 

*  *+  i  St  V 

VW  <233> 

ti.s>  siPChsT'  of  unit  cells  An  a  surface  layer  and  j  iden- 

0  j-  • 

r t'm>  oorrs i  lor  given  q^.q^.  A  reduced  dynamical  matrix 

can  bo  cufxr-e-i  «e  i-saents  are  specified  by 

i  • 

•  i[?i  ^*19(^2”^  3 

LD(q,q2)3  -  S  D  (U;VX‘)8  1  ^  3  . 

3  J  A  ^  (234  ) 

2 

The  alternative  expressions  for  (u  (4x) >  ,ean  now  be  rewritten  as. 

U  ^  j 

<u*(tx>>  -  (l/mxK6)  ss  x  CUCqJl/u*2^)  (235  ) 

3  ^ 


<u2U*)>  "  (ft/2a  N  )  S  [D**  (q)coth(AD^  (q)  /2ItRT}]  .  (236  ) 

3  -  •  h**‘h** 


If  L  is  the  number  of  layers  in  the  crystal,  the  sine  of  tho  reduced 
dynamical  matrix  is  3rL  x  3rb  whore  r  is  the  number  of  atoms  per 
unit  coll,  whereas  the  original  dynamical  matrix  is  3rLh's  x  3rLKg. 

For  a  crystal  20  atomic  layers  on  each  edge,  the  reduced  dynamical 
matrix  is  smaller  by  a  factor  of  400  than  the  unreduced  dynamical 
matrix.  One  can  save  considerable  computer  time  by  working  with 
the  reduced  matrix.  , 

ffo  now  discuss  specific  calculations.  An  early  theoretical 
discussion  of  surface  moan  square  displacements  was  given  by 
Kalashnikov^*®®'  who  used  the  Debye  isotropic  continuum  model. 

Re  found  an  anisotropy  of  tho  surface  mean  square  displacements 
in  qualitative  agreement  with  the  experimental  results  of  Zaasha 
and  Kalashnikov^143^  on  silver.  Recently,  Dennis  and  Hubor^170^ 
havo  given  theoretical  expressions  based  on  the  continuum  model. 

The  first  lattice-dynamical  calculation  was  reported ‘by  Rich^171^ 
who  used  a  nearest-neighbor  model  tp  study  tho  (001)  surface  of  a 
simplo  cubic  lattice  anti  found  a  10-30%  larger  mean  square  displace¬ 
ment  for  a  surface  atom  than  for  an  interior  atom  over  a  temperature. 
Tho  anisotropy  of  tho  mean  square  displacement  at  tho  surface  was 
not  revealed  by  Rich's  calculation  because  of  tho  special  character 
of  his  model.  Additional  calculations  based  on  a  noarost-noighbor 
model  havo  boon  carried  out  by  Corciovei  and  coworkers J17^  by 
Cel7,^173^  and  by  Kotheri  and  Siugal.^  174  ^  The  latter  authors 

considered  finite  aicrocrystallitos  and  calculated  the  mean 
square  displacements  for  atoms  on  surfaces,  edges,  and  corners. 

As  one  might  surmise,  the  mean  square  displacements  increase 
successively  from  surface  to  edge  to  comer.  A  related  problem- 
th*  scan  equor*  displacements  of  atoms  on  steps  on  surfaces  - 
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has  boon  attacked  by  Masri,  Allan,  and  DobrzynskiVii*  'for 
a  [100]  step  on  a  (001)  surface  of  a  simple  ‘cubic  crystal  using 
the  Roser.stock-Nowoll  model.  They  find  that  the  ratio  of  moan 
equaro  displacements  for  an  atom  on  the  stop  edge  and  an  atom 
on  the  flat  surface  is  ~  1.3  at  temperatures  coraparnble  to  the 
Dobye  temperature. 

iiaradudin  and  Melngailis* l75)  have  studied  tho  (001)  surfaco 
of  n  simple  cubic  lattice  with  nearest  and  next-nearest  neighbor 
central  forces  which  is  isotropic  in  the  long  wavolcngth  limit. 
Their  procedure  is  based  on  the  Green's  function  formalism  of 
Saction  3  „p  the  surface  is  created  as  a  perturbation  on  an 

Oht.T’-ise  perfect  periodic  crystal.  The  mean  square  displacement 
component  of  an  atom  at  lattice  site  t,  (only  one  atom  por  unit  coll 
is  assumed)  is  given  by 

<«?<£>>  ■"  “(V71^  S  S  oJWejvy)  x 
®  “  kj  k*j«  *  ~  ®  ~ 

A/  AW 


i(k+k')*xU>  &  ...  ,  ‘  2x 

a - ~~  h  u(kj,k’j’;-u)“) 

^  ^  u 


(237  ) 


where  u(£4  .**4’ is  given  by  Eq.  (206)  and  un  -  (2nnkgT/h) 
with  n  an  integer.  In  the  high  temperature  limit,  Eq,  (237)  re¬ 


duces  to 


<u“(<,»  -  (k  T/Nra)  S  2  o  (kj)e  (k!^')  x 
°  B  kj  k'j*  a  ~  *  ~  J 

•  x  _  t(k3,k'3»;°)  j  ^  ^ 

whore  t<kj,k'j*;0)  is  given  by  Eq.  (158)  with  »2  -  0.  Using 
Eq.  (  238  ) ,  K&radudts  and  We Inga ills  have  evaluated  the  mean  square 
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displacements  for  their  model  numerically  on  a  computer.  In  addi¬ 
tion  to  the  features  found  by  Rich,^17*^  thoir  results  indicate 
that  at  tho  surface  the  Debye- Haller  factor  is  anisotropic  with 
the  component  of  mean  square  displacement  perpendicular  to  tho 
surface  significantly  larger  than  tho  parallel  component.  The  mean 
square  displacements  approach  the  bulk  value  very  closely  at  five 
or  so  atomic  layers  from  tho  surface. 

Calculations  based  on  tho  high  temperature  limit  of  Eqs.  (235  ) 
and  (  236)  havo  been  carried  out  by  Clark,  Herman,  and  Hall is ^17®  ^ 
for  a  noarest- neighbor  central-forco  model  of  a  face  centered  cubic 
lattice.  Tho  (001),  (110),  and  (111)  surfaces  were  studied.  It 
was  found  that  tho  normal  component  of  mean  square  displacement 
is  about  double  tho  bulk  value  for  all  throe  surfaces.  The  (110) 
surface  roveals  an  anisotropy  in  the  tangential  components  not 
found  in  the  other  surfaces.  This  anisotropy  together  with  tho 
rapid  fall  off  of  the  mean  square  displacement  with  increasing 
distance  from  the  surface  ore  illustrated  in  Fig.  13.  The 
results  for  this  surface  are  summarized  in  Table  3  and  compared 
with  tho  experimental  results  of  UacRao^  146  ^  for  nickel. 

•  Table  3 .  Theoretical  and  experimental  mean  square 
displacement  components  (relative  units) 
for  atoms  at  a  (110)  surface  and  in  the 


bulk  of 

a  fee 

crystal. 

[no] 

[Ho] 

[001] 

Bulk 

Theoretical 

0.80 

0.64 

0.86 

0.40 

Experimental 

1.41 

0.63 

1.41 

0.45 

The  agreement  between  theory  and  experiment  Ls  qualitatively  satis¬ 
factory,  but  quantitatively  leaven  something  to  be  desired.  Sinca 


tho  nearost-neighbor  model  gives  a  reasonable  fit  to  the  phonon 
dispersion  curves,  the  discrepancy  probably  lies  elsewhere.  One 
possibility  is  multiple  scattering  effects.  Another  is  that  the 
.force  constants  at  the  surface  are  different  from  those  in  the 
bulk..  This  possibility  has  been  Investigated  by  Wallis,  Clark, 
and  Herman^  *77?  who  found  that  agreement  of  theory  with  MacRao's 
data  can  be  significantly  improved  by  changing  .to  surface  force 
constants  on  the  order  of  50  per  cent.  Further  work  is.  required, 
however,  to  establish  whether  such  changes  are  meaningful.  As 
mentioned  in  Sec.  2,  ,  Clark  et  alv  1  have  calculates  changes 
of  25-30  psr  cent  in  the  surface  force  constants  for  a  rsodol  of 
iron  (bcc  lattice)  with  nearest  and  next-nearest  noighbor  Lennard- 
Jones  interactions,  so  relatively  large  changes  at  the  surface  do 
soca  possible. 

*  ( 178  1 

Recently,  Rovida,  Torrini,  and  Zanazzi'  7  have  used  a 
Mcrse  interaction  potential  to  investigate  the  (110)  surfaco  of 
nickel.  They  allow  for  surface  relaxation  and  thus  take  into 
account  force  constant  changes  near  the  surface.  Taking  into 

account  surfaco  relaxation  leads  to  significant  improvement  in 
the  agreement  between  the  calculated  mean  square  displacements 
and  the  experimental  results  of  MacHae. 

At  temperatures  sufficiently  small  compared  to  tbo  Debye 
tempera ture ,  the  mean  square  displacement  of  an  atom  should  ex¬ 
hibit  quantum  mechanical  zero-point  notion  effects.  This  behavior 
has  been  revealed  in  calculations  by  Wallis  et  aiv  *  based  on 
Sqs.  (235  )  and  (236  )  for  the  (110)  surface  of  a  faco-eentered 
cubic  lattice  with  nearest- ns igkbor  central  forces.  Additional 
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calculations  of  the  temperature  dependence  of  tho  :  ?an  square  dis¬ 


placements  for  a  face— centered  cubic  lattice  have  been  carried  cut 
by  Allen  and  DeWette^®  }  who  treated  the  (100),  (110),  and  (111) 


surfaces  using  a  nodol  involving  Lennar?!- Jones  interactions'  be¬ 


tween  all  pairs  of  atoms.  In  a  subsequent  paper,  these  ?uthorsf“  ^ 


point  out  that  tho  ratio  of  surface  to  bulk  mean  squaro  displace¬ 


ments  is  a  nearly  universal  function  of  T/e^.  Their  results  are 


showt  in  Fig.  14  .  The  larger  value  of  this  ratio  at  the  higher 


temperatures  may  be  attributed  to  either  the  heavier  weighting  of 


surface  modes  at  high  temperatures  thhn  low  temperatures  or  to  tho 


earlier  onset  of  zero-point  notion  in  the  bulk  than  at  the  sur¬ 


face  ns  the  temperature  is  decreased. 


In  the  work  just  mentioned,  soma  account  of  inharmonic  effects 


was  included  through  the  use  of  the  qausiharnonic -approximation. 


These  effects  tend  to  increase  the  surface- to- bulk  ratio  through 


decreases  in  the  surface  force  constants.  A  acre  complete  treat¬ 


ment  of  inharmonic  effects  is  contained  in  tho  molecular  dynamics 


calculation  of  Allen,  DeWette,  and  Rahman^18®  ^  who  solved  the 


equations  of  motion  on  a  computer  for  atoms  interacting  with  a 


Lennard-Jones  potential  in  a  fee  lattice  with  (100) ,  (110) ,  Or 


(111)  surfacos.  The  surface  aaan  square  displacements  at  half 


the  melting  temperature  were  found  to  be  1(K20  par  cent  larger 


than  those  calculated  using  lattice  dynamics. 


The  calculations  described  so  far  have  either  boon  restricted 


to  nearest- no ighbor  interactions  (Rich^17*^)  or  have  involved 


central  forces  so  that  Cauchy’s  relation  on  the  elastic  constants 


C12  "  C44’  l^.-ds.  A  more  general  codel  consisting  of  nearest  and 
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next-neareat  neighbor  central  forces  plus  angle-bonding  interac¬ 
tions  involving  pairs  of  nearest  neighbors  has  been  employed  by 
Wallis,  Clark,  and  Herman^ 181  ^  to  the  (001)  surface  of  a  fee 
lattice  in  the  high  temperature  limit.  The  force  constants  can 
be  chosen  to  give  an  exact  fit  to  tho  elastic  constants  for  a 
given  fee  material,  and  Cauchy's  relation  need  not  bo  obeyed. 

The  results  for  the  ratio  of  surface- to-bulk  mean  square  displace¬ 
ments  are  remarkably  independent  c_  the  material  for  tho  normal 

i 

component,  ranging  only  botweon  2.01  and  2.03  for  Al,  Pb,  Pt,  Pd, 
Ag,  Au,  Cu,  and  Ni.  The  ratio  for  the  parallel  component  shows 

core  variation,  ranging  from  1.21  for  Au  to  1.51  for  Hi.  The  re¬ 
sults  for  the  normal  component  ratio  arc.  in  reasonable  agreement 
with  experiment  (cf .  Table  1  )  for  the  (001)  surface  of.  Cu  and 
Ni,  but  are  only  about  one-half  the  experimental  values  for  Pt 
and  Pd.  The  latter  discrepancy  may  indicate  weak  binding  of  th© 
surface  transformations  can  take  place  in  Pt^  149  ^  and  Pd^50  \ 

• 

A  model  similar  to  that  described  above  but  including  third 
neighbor  central  forces  has  been  used  by  Wallis  and  Cheng ^ 
to  calculate  tbs  temperature  dependence  of  the  mean  square  dis¬ 
placements  of  atoms  on  the  (001)  and  (110)  surfaces  of  chromium. 

* 

The  result  for  the  surface-to-bulk  ratio  of  tho  normal  components 

at  the  (110)  surface  is  in  qualitative  agreement  with  the  experi- 

(  160) 

mental  results  of  Kaplan  and  Soaorjai  ,  but  quantitative 

agreement  will  require  inclusion  of  multiple  scattering  effects 
or  temperature- dependent  changes  in  the  surface  force  constants. 

The  result  for  the  ratio  at  the  (001)  surface  in  the  high  tempera- 

* 

lure  regime  is  significantly  less  than  that  iousig  experimentally 


by  labor,  Wilson,  and  Basto**  Horo  work  is  required,  both 

experimentally  and  theoretically,  on  this  problem. 

An  interesting  approximate  procedure  for  calculating  mean 
square  displacements  has  been  proposed  by  Hasri  and  Dobrzynskl. 

The  dynamical  matrix  D  can  be  written  as  the  sum  of  the  diagonal 
part  d  and  the  off-diagonal  part  R: 

D  -  d  +  R  (239  ) 

If  Eq.  (  239  )  is  substituted  into  Eq.  <  231 )  and  the  result  ex¬ 
panded  in  a  matrix  power  series  in  R,  the  result  for  tho  first  few 
terms  in  tho  expansion  of  the  mean  square  displacement  As 

<uf  C4x)>  «*«  (ft/2a  )  ]  d“^  coth  a  + 

**  A  }  ^  . 

(  -tXCr.iXOf 

*  J -*  [*  i  2  °°a  ~  *  I  * 

3 

+ -|(ft/kgT)2  d""^  coth  a. sinh"2a  +  •••( 

-Uxcr.-Uo  •  j(240  ) 

whore  t  “  fcd^/2kQT.  Since  it  is  very  easy  to  calculate  the  square 
root  and  other  functions  of  a  diagonal  matrix  and  since  only  non- 
negative  integer  powers  of  the  non- diagonal  matrix  R  are  required, 
the  right-hand  side  of  Eq.  (  240)  is  readily  evaluated.  Masri 
and  Dobrzynski  treated  tho  (001)  surface  of  a  simple  cubic  lattice 
using  Rich’s  model ^  ^  and  obtained  results  accurate  to  about  20 

per  cent.  In  his  thesis.  Masri^  134 ^  has  employed  this  procedure 
to  calculate  tho  mean  square  displacements  as  functions  of  tempera¬ 
ture  for  atoms  on  various  surfaces  of  a  number  of  face-contered  and 


body-contorod  cubic  metals.  This  work  has  recently  oppeorod  in 
papers  by  Maori  and  Dobrzynski^  '  on  the  face-centered  cose 
and  by  DobrzynBki  and  MnBri^  186  *  on  the  bady-centored  case. 

/  1  0*9  \ 

The  procedure  has  also  been  applied  by  Thcjsten  and  Dobrzynski  * 

to  the  calculation  of  mean  square  displacements  of  atoms  on  the 
(100)  and  (111)  surfaces  of  diamond,  silicon,  and  germanium. 

The  effect  of  removing  rotational  invariance  from  the  force- 
constant  model  has  been  studied  by  Musser^  188  ^  for  the  (001)  sur¬ 
face  of  a  fee  crystal.  The  non-rotationally  invariant  non-central 
forces  were  used.  Rotational  invariance  was  then  imposed 
by  modifying  the  surface  force  constants  using  the  procedure  of 
Ludwig  and  Longoler. ^  Musser  found  that  the  nox-rotationally 

invariant  model  gives  larger  mean  square  displacements  thm  the 
rotat innally  invariant  model.  For  the  perpendicular  component  of 
Pd,  the  increase  is  rather  sizeable,  ~3?%,  while  for  Ni»  it  is 
much  less,  ~€&.  Musser  exploited  matrix  partitioning  methods  te . 
reduce  the  size  of  the  matrices  to  be  inverted. 

Theoretical  studies  of  mean  square  displacements  of  adsorbed 

< 

189) 

atoms  have  been  carried  out  by  several  workers.  Dlouha  has 

utilized  the  Rosenstock-Kewell^^modol  to  study  thin  films  on  the 
(100)  surface  of  a  sioplo-cubic  lattice.  Allen,  Alldrod-o,  and 
Do  Wotte^^  have  employed  the  LennardrJones  interaction  to  study 
the  (100),  (110),  and  (111)  surfaces  off  a  monatomic,  face-centered 
cubic  crystal.  Relaxation  of  the  atoms  near  the  surface  was  taken 
intc  account.  They  find  that  the  mean  square  displacement  of  the 
adsorbate  &toa  is  enhanced  or  depressed  according  to  whether  the 

adsorbate  atomic  sans  is  smaller  or  larger  than  the  substrate 

•  ... 


atomic  mass,  respectively.  Tho  effect  incrcasos  as  the  temperature 
is  lowered. 

Recently,  Theoten,  Dobrzynski,  and  Doraange^®0^  have  studied 
the  <001)  surface  of  Ni  with  a  (2x2)  adsorbed  layer  where  the 
absorbod  atoms  are  S,  Na,  Cs,  or  0.  They  calculated  tho  mean 
square  displacements  parallel  and  perpendicular  to  the  surface 
for  tho  adsorbed  atoms  and  compared  the  results  with  experimental 
data.  They  find  that  the  adatom  particularly  affects  the  "extra 
spots"— o.g. ,  {£$) — but  not  the  ordinary  spots.  For  S  on  Hi  (001) 
Theeten  et  al.  conclude  that  the  Ni-S  force  constant  is  larger 
than  the  Ni-Ni  nearest  neighbor  forco  constant. 

The  theoretical  work  discussed  so  far  has  been  restricted 

to  non-ionic  crystals.  Recently,  Chen,  Alldredge,  Do  Wette,  and 

Allen  *  have  reported  calculations  for  the  (100)  surface  of 

NaCl  using  the  rigid  ion  model.  Interplanar  relaxation  near  the 

surface  was  allowod,  but  was  assumed  to  be  tho  same  for  cations 

and  anions.  Tho  dopondencos  of  the  mean  square  displacements  on 

distanco  from  the  surface  and  on  temperature  are  qualitatively 

the  same  for  cations  and  anions  and  also  the  same  as  for  non-ionic 
» 

crystals.  However,  the  magnitudes  of  the  surface-to-bulk  mean 
square  displacement  ratios  are  about  20  per  cent  smaller  than 
for  non-ionic  materials. 
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D.  iloan  Square  Velocities  of  Surface  Atoms 


Closely  related  to  moan  square  displacements  aro  mean  square 


velocities.  The  mean  square  velocity  of  an  atom  is  in  principle 

{ 182} 

measurable  through  the  second  order  Doppler  shift  in  the 


lltissbauer  effect.  The  energy  2  of  the  emitted  gamma  ray  is  shifted 


by  an  amount  A£  given  by 


AS  -  E«u2>/2c2} 


(241  ) 


2  i 

whore  (Cs  }  is  the  mean  square  velocity  of  the  emitting  nucleus. 


By  experimentally  measuring  both  E  and  AS,  one  can  obtain  an  experi¬ 


mental  value  for  <u  }.  Such  measurements  can  be  used  to  obtain 


information  about  the  coupling  constants  between  the  active  atom  and 


its  neighbor f*®3  ^  The  procedure  has  the  significant  limitation  that  It 


57  119 

is. restricted  to  suitable  radioactive  atoms  such  as  Fe  and  Sn. 


Theoretical  expressions  for  the  mean  square  velocity  can  be 


obtained  in  n  fashion  similar  to  that  for  the  mean  square  dis¬ 


placement.  In  a  normal  mode  of  vibration  of  frequency  uip,  the 


volociry  and  displacement  components  of  atom  -tx  aro  related  by  the 


equation 


*  Ux)  -  ia-  u  (Ax)- 

Of  P  Of  ^ 


(242  ) 


Consequently,  the  expression  for  (tii(tx) >  analogous  to  £q.  {  230  ). 


taken  the  form 


<trUsO>  -  a/aj  2  K  Ux)  j2  ttoJ 

«  ~  X  _  ’  so  ~  1  p 


{  243  ) 
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Transforming  this  oquatior.  using  tho  theory  of  matrices,  we  obtain 

<^Ux)>  -  (ft/Sta^LO*  COthUD^T))^^  .  (244)  | 


Equations  analogous  to  Eqs.  (  233)  and  (236  )  follow  Iron  tho  saso 
approach . 

It  is  instructive  to  considor  the  expansion'  of  £q,  (244  ) 
appropriate  to  high  temperatures: 


«#£*>>*  r  + 

X 


_ u _ _  rjyi  +  .... 

12fnxkT  LUJtxo,ixcr 


(245  ) 


Wo  note  that  in  the  extreme  high  temperature  limit,  the  mean  square 
velocity  has  the  oquipartifcion  value  kT/o^  which  is  independent  of 
the  position  of  the  atom  relative  to  the  surface.  The  correction 
term  in  1/T  is  proportional  to  the  dynamical  matrix  0,  sq  expert-  , 
mental  measurement  of  this  term  will  give  the  force  constants  at 
the  surfaco  directly.  The  primary  obstacle  is  the  extreme  difficulty 
of  making  such  experimental  measurements  ■'■>r  surfaco  atoms. 

At  lew  temperatures,  the  expansion  of  the  we.-n  square  velocity  starts 


>  “  <•*»,> 


(246  ) 


The  xoro-point  velocity  consequently  will  depend  on  tho  proximity 
of  the  atom  to  the  surface  through  the  D*  factor  in  £q.  (248). 

Tho  situation  is  illustrated  in  Eig.  15  for  a  one-dimensional 
monatomic  lattice  with  free  ands.^®®  ^  The  moan  square  velocity 
is  smaller  for  the  end  atoms  than  for  the  intorior  stoss  an  a 
consequence  of  the  smaller  nunbor  of  interactions  involving  tho  end 
atoms  compared  to  interior  a toys. 


•Maradudin  and  Molngnili.u**76  ^  hnvo  developed  oxprossiono 
for  the  neon  equaro  velocity  of  a  surface  atom  in  n  Bravais  lattice 
using  a  Green's  function  approach;  The  result  is 


<d2U)>  -  {keT/Nn)  £  £  o(kj)-e(k‘J*)  x 

tt  k'J*  ~~  ~~ 


i(k+k’>  *x(-t)  *  o 

5(kj)u)(k’j‘>o  ~~  ~~  £  u(kj,J  j’l-u)2) 

^  n~-*9  ~  * 


'v  247  ) 


•shore  utkJ.k’J5 ;-u>“)  is  given  by  ISq.  (*J06  >  with  -  (2nnk0T/h) 
and  n  an  integer.  High  temperature  axpansioss  such  as  £q.  '(  245  ) 
can  bo  derived  from  Eq.  (247) . 

Detailed  calculations  of  mean  squaro  velocities  for  atoms  or! 
a  (110)  surface  of  a  fee  crystal  have  been  carried  out  by  "jftllls 

et  al*162-  )-  using  a  nearest  neighbor  central  force  model. 

> 

Anisotropy  shows  up  only  at  low  temperatures,  sinco  each  component 

• 

has  the  e^uipartition  value  kT/a  at  high  temperatures.  The  smaller 
number  of  boitfb  which  connect  a  surface  atom  to  its  neighbors  than 
connect  a  bulk  atom  with  its  neighbors  leads  to  a  smaller  moan 
squaro  volecity  for  tho  surface  atom  at  low  temperatures , 

Very  few  experimental  investigations  of  the  socond-ordor 

f.- 

Doppler  shift  for  surface  atoms  have  been  made.  Allen*  165  ^ 

<r 

studied  the  temperature  dependence  of  this  shift  for  57?e  nuclei 

t 

both  in  the  bulk  and  on  the  (100)  and  ^1X1)  surfaces  of. silicon  ' 
and  found  that  a  surface  mean  square  velocity  roughly  five  per  cant 
loss  than  in  the  bulk.  Similar  qualitative  results  have  been  ob~  • 

tained  by  Burton**66  ^  for  577«  os  tungsten  and  by  Godwin*194  * 

S7  ^ 

for  ?o  on  silvar. 
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€.  INELASTIC  SCATTERING  FROM  CRYSTAL  SURFACES  - 
A.  Introduction 

A  very  attractive  method  for  determining  surface  phonon  onergios 
is  tho  inelastic  scattering  of  particles  such  as  electrons  or  helium 
atoms  from  crystal  surfaces.  Helium  atoms  would  be  particularly  suit¬ 
able,  ainco  they  are  strongly  scattered  by  the  surface  and  since 
their  onergy  is  on  the  order  of  millivolts  -  i.e.,  in  the  range  of 
phonon  energies  -  when  their  DaBroglie  wave  length  is  on  the  order 
of  a  typical  lattice  spacing.  Unfortunately,  the  experimental  diffi¬ 
culties  are  rather  severo.  Neutrons  are  of  little  use  because  they 
penetrate  dooply  .into  tbo  crystal.  Electrons  do  not  havo  tho  favor¬ 
able  energy  -  DeBroglio  wave  length  relationship  of  helium  atoms  but 
are  experimentally  more  tractable,  and  some  success  has  been  achieved 
by  their  use. 

B.  Theoretical  Development 

Formal  treatments  of  the  inelastic  scattering  of  particles  from 

{  197  j 

crystal  surfaces  have  been  given  by  Cabrera,  Celli,  and  Hanson, 
Cabrera,  Celli,  Goodman,  and  .Hanson \  198  \  Goodman/199  *  Roundy  and 

Mills^200  *  Lucas  and  Sunjic/201  * ,  and  Evans  and  Mills/202  * 
Cabrera  ot  a/  197  ^  employed  the  T-natrix  formalism  to  derive  an  ex¬ 
pression  for  tho  one-phonon  reflection  coefficient  which  was  then 
evaluate*/ 19S  ^  using  a  Xorse  interaction  potential.  Specific 
calculations  ware  made  for  the  cases  of  helium  incident  normally  on 
silicon^  203  ^  and  on  Ll/  204  *  with  the  aid  of  a  Debyo  continuum 
modal  for  the  lattice  vibrations.  The  differential  reflection  co¬ 
efficient  as  s  function  of  emission  angle  exhibits  a  continuous 
variation. up  to  a  cut-off  angle  and  then  a  dolta-function  peak  split 
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,  off  b>r  a  finite  gap.  The  continuous  part  la  associated  with  the 

creation  of  bulk  and  mixed  phonons,  while  the  delta-function  is 

associated  with  tho  creation  of  surface  phonons  (Rayleigh  vavos) . 

Those  calculations  havo  been  extended  by  Goodman^  '  to 

tho  case  of  He  atoms  incident  on  silver.  The  results  are  found  to 

bo  in  reasonable  agreement  with  experimental  data  of  Subbarao  and 
(  206  ) 

Hiller.  One-phonon  inolastic  scattering  has  boon  observed 

f  207  \ 

oy  Fishor  and  81odsoov  '  for  tho  He  atoms  incident  on  tho  (001) 

surfaco  of  LiK.  It  should  bo  mentioned  that  there  is  -an  extensive 
literature  on  the  scattering  of  atoms  from  solid  surfaces  which 
has  been  loviewed  recently  by  Goodman.^  208  ^  Additional  references 

include'  the  work  of  Chambers  and  Kinzer-^208  ^  and  of  Trilling^210  ^ 
who  treat  tho  solid  as  a  semi-infinite  elastic  continuum.  The  dis¬ 
crete  atomic  character  of  tho  solid  has  been  takon  Into  account  by 

,  ..  {  211  }  (212  "j 

Araand  and  by  Beeby  and  Dobr&ynsRi .  ' 

(  200  ) 

Rcundy  and  Hills  have  developed  a  Green's  function  theory 

of  ono-phonon  inelastic  scattering  of  low-energy  electrons  from  . 
crystal  surfaces  where  the  electron  interacts  with  tho  atoms  of  the 
crystal  through  short  range  forces.  The  interaction  potential  onorgy 
of  tho  electron  with  the  atom  at  tho  position  R (•£.)  may  bo  expanded 
in  tho  form 

iq*(r-R(t)) 

vl  <£  "  SCO)  -  (1/0)  £  v,  (q)o -  (  248  ) 

z  q  z 

whoro  it  is  assumed  that  the  surface  is  in  tho  xy  plane-,  the  sub¬ 
script  4^  takes  into  account  that  .the  potential  for  an  atom  near 
tho  surface  may  bo  different  than  for  an  atom  in  tho  bulk,  and  n 
is  a  quantisation  volume.  The  incident  and  scattered  electrons  have 
energies  eo  and  «s  and  wavs  functions.  f0  and  which  can  bo  written 
in  the  ferns 


*V' 


*  '  4*w  - 


*0<£a  .»C  "  e  ~S  JVo)(ka°)  »fc0;£°,z) 

♦s(J«*2)  "  8  ~C  ~Vs)<jc<s>,es;r  ,z> 


(248a  ) 


(248b  ) 


whore  r  is  the  two-dimensional  position  vector  parallol  to  .the 
~a 

surfaco,  and  k^3^-  aro  two-dimensional  wave  vectors  parallel 

~a  ~i 

to  the  surface  for  the  incident  and  scattered  electron,  and  the- 
functions  ,z)  are  periodic  in  rg .  Outside  the 

crystal,  Roundy  and  Kills  write  ^s(£g,s)  in  the  fora 


.  ,  \  n-lr  i(iS.(s)  +  £»>'£«  ikz{&>)5 

,z)  -  n  5  S  e  a  o 


5  »j&) 


(249  ) 


whore  a(k|s)  ,jgu>  is  a  Fourier  amplitude,  kz(j?fi)  “  [Zn2s  *  -( }«b f -*-£3 5 2 j ^ » 
and  gj  is  a  two-dinonsional  reciprocal  lattice  vector  associated 
with  the  surface. 

In  order  to  describe  tha  phonons  emitted  or  absorbed  in  the 
scattering  process,  it  is  convenient  to  introduce  the  normal  co¬ 
ordinate  transformation 
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where  u (f.)  is  the  displacement  of  the  atom  at  site  l  from  equilibrium, 

i,*. 

u?p  and  j[p(£)  are  the  frequency  and  eigenvector  for  tho  p*  normal 

mode  of  tho  crystal  with  the  free  surface,  and  at  aro  the  phonon 

P  p 

annihilation  and  creation  operators,  and  M  in  tho  atomic  mass.  Tha 
translational  periodicity  parallel  to  tho  surface  can  be  exploited  by 
writing 
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whore  K  is  the  number  of  atoms  in  the  surfacr  layer  and  j  specifies 
8 

the  normal  nodes  for  given  value  of  tho  wave  vector  kg  parallel  to 
the  surface.  The  sun  over  p  in  %q.  (250  )  now  bocoaes  a  sun  over 

k  and  J.  Ono  effect  of  the  theraal  vibrations  is  to  re normalize 

~g  i 

tho  interaction  potential  throu  ‘b  a  Debye-Wallor  factor 

.2. 


Js>  -  v^(s)exp[-J((S*«<£>)  >T3 


(253  ) 


rhero  denotes  a  theranl  average  over  a  canonical  ensemble. 

The  scattering  cross  section  involves  an  electronic  matrix 
element.  A  convenient  measure. of  this  matrix  element  is  provided 
by  the  quantity  . 
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where  Gg  is  a  two-diracnsiotial  reciprocal  lattice  vector  that  places 
kf°3  -  +  Gg  in  the  first  Brillouin  zone,  X  is  defined  by  its 

components  Kj  -  kg°^  -  kg s ^  +  Gg ,  Xs  *»  q2,  and  the  range  of  Integra 
over  z  signifies  that  tho  crystal  lice  in  the  region  z  <  0, 


4 


Let  us  restrict  our  attention  to  phonon  emission.  The  equations 


of  conservation  of  energy  and  of  wavo  vector  component  parallel  to 
the  surface  take  the  forms 

«s  “  eo  "  (255  ) 

and 

feiS>  "  *i0)  "  +  £s  <256  ) 


The  specification  of  the  directions  of  the  incident  and  scattered 
electron  is  completed  by  the  equations 

k^o)  -  ~  [2neo  -  <kSo))2]*  (257  ) 

and 

kgS*(g|)  **  t2acs  -  {kj}s*  +  £a)23^  (258  ) 


A  significant  simplification  of  tho  analysis  occurs  if  one  noglccts. 
the  energy  of  the  phonon  created  relative  to  tho.  energies  of  tho  inci¬ 
dent  or  scattered  electron.  Then  a  Rayleigh  wavo  and  a  bulk  phonon 
with  the  samo  value  of  ke  will  scatter  an  electron  into  the  same 
differential  element  of  solid  angle  d  £}  .  •  • 

Tho  differential  cross  section  for  one-phonon  emission  can  he 
expressed  in  the  fora 
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where  k'  *  is  tho  magnitude  of  tho  scattered  wave  vector,  «  -  c  -  ce, 
a(ai)  is  tho  Bose-Elnstein  occupation  probability,  ,  ?  sa)  is 

tho  phonon  spectral  density  function. 
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*4  4  ,(••*•)  “  B.  (••*•)»*  (•••♦)*  (  280  > 

2  2'  v  t  , 

2  2* 

and  Q  -  kg  +  2  Q  (4  )  with  Q_(4„)  specified  by 

Q2(v  -  i  -5T-  log  n  (••••>  .  (  261  ) 

« 

The  phonon  spectral  density  function  is  related  to  the  phonon  Green’s 
function  Uffg0s.i;4s6E,;ai)  j,y 
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C  283  ) 


For  a  particular  direction  of  the  jcattered  electron  (5c *  fixed) 
aad  scattering  only  froa  tho  surface  layer  of  a  terns  ?  ono  can  arita  t« 
a  good  approximation 
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(  264  > 


where  K  is  independent  of  the  energy  leas  ©.  Houndy  and  llills  have 
.  *.  * 

evaluated "the  a,  dependence  of  the  cross  suction  for  a  (100)  surface 
of  a  fee  crystal  with  nearest  neighbor  central  forceB.  They  made 
very  sinple  approximations  for  the  periodic  parts  of  the  electron 
wave  functions  and  the  Fourier- transformed  potential  v(q) .  The  re¬ 
sults  arc  shown  in  Fig.  16,  Tho  delta  function  peak  is  associated 
with  the  excitation  of  Rayleigh 'waves,  while  tho  broad  band  at 
higher  energy  loss  is  associated  with  excitation  of  bulk  nodes. 

Thus,  the  surface  node  dispersion  curve  can  is  principle  be  deter¬ 
mined  Iron  energy  loss  xeasuroaonta  for  various  kj.  To  date,  however, 


no  one  has  succeeded  in  carrying  out  such  measurements  for  P.ayleigh 


Wo  now  turn  to  the  case  where  the  scattered  electron  interacts 
with  the  crystal  through  a  macroscopic  electric  field.  This  caso 
concerns  ionic  crystals  or  covalent  crystals  with  8urface"indnced 
offectlvo  charges  and  is  directly  comparable  to  experiment. 

A  classical  macroscopic  theory  of  energy  loss  by  elections  due 
to  excitation  of  optical  phonons  has  been  given  by  Fujiwara  and 
Oh  taka  who  adapted  the  treatment  of  Ritchie  for  plasmons. 

Cousidor  a  slab  of  an  ionic  crystal  between  the  planes  z  «*  0  and 
z  »  a  and  an  electron  with  velocity  v  incident  on  the  slab.  The  . 
energy  loss  per  unit  path  length  is  simply  given  by 


~  -  -  •£.  y.E 
dx  v  ~  ~m 


(  265  ) 


where  E_  is  the  electric  field  due  to  the  medium  alone.  The  total 
«*»  . 

field  £  is  given  by  £  -  -  vtp  where  the  potential  $  satisfies  the 


equations 


€(»)  V  tp  -  4uu  6  (r  -  vt)  ,  ‘  0  i 


(  266a) 
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.  *  <? 


4 ne 6 (v  "  vt)  elsewhere 


(  266b) 


and  g{a;)  is  the  frequency  dependent  dielectric  constant  given  by 
Eq.  {84  ).  Equations  (  266  )  are  to  be  supplemented  by  the  condi¬ 
tions  of  continuity  of  the  normal  component  of  the  electric  displace¬ 
ment  and  the  tangential  components  of  the  electric  field  at  the 
boundaries  of  the  slab. 

It  is  convoniont  to  introduce  the  Fourier  transform 
‘P(rl»Ex*t)  ~  (1/2n>  <p(rE,kA,o)e  \  (  2g7  } 
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vhoi-o  rj  and  r  are  coordinates  parallel  and  perpendicular  to  v  , 
respectively,  and  kA  is  a  two- dimensional  wave  vector  perpendicular 
to  v  .  Equations  (  266  )  can  now  be  solved  together  with  the  boundary 
conditions  and  the  energy  loss  computed  using  Eq.  (265)  .  For  normally 
incident  electrons,  tho  probability  of  loss  of  energy  ft tu  can  be 


written  as 


P(k^,  »}  -  -  U/fta)  dr9(e/v)v.Em|ri  „  yt 


(  268  ) 
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x  2<g~  l)cos<ma/v)  +  (£-  IJ^xpH^a)  +  U“  ^)oxp(kAa). 


(£  -  l)“exp(-kAa)  -  (£  +  D^expO^a) 


(  269  ) 


•  1 
j 


The  first  term  on  the  right  band  side  of  Eq.  (269)  corresponds  to 
excitation  of  bulk  LQ  phonons  and  the  second  term  to  excitation  of 
surface  optical  phonons.  Setting  the  denominator  of  tho  second. tera 
to  zero  gives  exactly  the  dispersion  relation  for  the  Fuchs-Kliewer 
surface  optical  phonons,  Eq.  (83  ),  in  tho  unretarded  limit. 

.  v/ 

Lucas  and  Sunjic  * '  have  given  a  seniclassical  treatnont  of 
the  problem  for  polar  crystals.  The  electron  is  treated  as  a  classi¬ 
cal  point  particle  which  moves  on  a  fixed  trajectory  at  constant  speed. 
W©  have  soon  in  Section  2  that  the  Fuchs-Kliower  typo  surface  optical 
node  has  an  electric  field  which  ponotrafcee  outside  the  crystal*  The 
electron  couples  to  this  electric  field -and  drives  the  lattice.  The 
interaction  energy  between  the  electron  and  tho  surface  sods  in  the 
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non-rotarded  limit  is  simply  ecp(r)  whore  c?(r)is  tho  electrostatic 
potential  sot  up  by  the  surface  modo.  The  probability  of  oxcitation 
of  one  or  more  surface  optical  phonons  is  then  calculatod  quantum 
nochanically.  For  a  slab  of  thicknoss  2a,  Lucas  and  Sunjic  find 
that  tho  probability  of  the  electron  losing  an  energy  h<u  is 


/»  .  .  r  -iu>  t  -ia_t 

dt  expC  /  dk.(Q+e  *  +  Q_e  )] 

— u# 


C  27^  ) 


where  PQ  is  the  no-loss  probability,  tuA  are  solutions  of  the  Fuchs- 
K1 tower  dispersion  relation  for  surface  optical  phonons, 


sinh(2ka)  “  ±  2€(u») /C€2(is)  -  1]> 


(  271  ) 


*'si 


k  "  |kg|,  and  €(<u)  Is  the  dioloctric  constant  givon  by  2q.  (84  ). 

Tho  quantities  are  specified  by 

q  (ka)  .  0~2ka  - 1_ - 
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where  (1  J'ija  the  ion  plasma  frequency  and  tho  vectors  v  -  (vs,v  )  and 
Pv/  •  ***  ** 

v*  •*  (vi oro  tho  volocitios  (rosolvod  into  components  parallel 

~  -**A 

and  perpendicular  to  the  surface)  of  tho  incident  and  scattered 
electron . 

In  general,  Eqs.  <  270  )  and  (  272  )  describe  a  Poisson 

distribution  for  the  excitation  of  n  phonons.  For  the  case' of 
specular  reflection  from  a  very  thick  slab,  the  probability  for 
the  excitation  of  ono  phonon  beceaetv 
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where  G(|vcj,vA)  is  a  relatively  slowly  varying  function  of  its 
arguments .  From  £q.  (  273  )  wo  see  that  tho  probability  varies 
with  energy  E  of  the  incident  electron  as  E  ^  to  a  good  approximation. 
This  dependence  is  in  agreement  with  tho  experimental  results  of 
Ibach^215^  on  ZriO.  For  i  review  of  this  topic,  see  the  recent 
article  by  Lucas  and  Sunjic*2'1'6^  in  this  Journal. 

A  quantum  mechanical  treatment  of  inelastic  electron  scatter- 

ing  from  surface's  of  polar  crystals  and 'covalent  crystals  with 

•  -  (  202  ^ 
surface-induced  charges  has  been  givan  by  Evans  and  Mills:  ' 

For  polar*  crystals,  their  results  reduce  to  those  of  Lucas  and 

Sunjic  in  the  specular  case.  Fox*  c<  /alent  crystals,  the  surface 

optical  modes  are  always  localized  within  a  few  atomic  layers  of 

the  surface.  Evans  and  Mills  find  thqt  the  angular  dependence  of. 

the  one-phonon  scattering  crocks  section  is  less  for  tho  covalent 

crystals  than  for  the  ionic  crystals.  Furthermore,  tho  energy 

dependence  of  tho  cross  section  for  the  covalent  case  is  more 

complicated  than  tho  E  2  dependence  for  tho  ionic  case.  These 

differences  con  in  principle  be  used  to  identify  the  typo  of  surface 

node  involved  in  the  scattering.  .Recently,  Evans  and  Mills^217^ 
have  extended  their  work  to  auliiphonon  processes, 

C.  Experimental  Results 

Anong  the  earliest  observations  of  surface  vibrational  modes 
by  inelastic  electron  scattering  is  that  of  BoerccJi,  Geiger,  and 
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Stickolv  '  on  LiF,  Those  workers  used  relatively  high  energy 
electrons  (25  KeV) ,  but  wore  able  to  achieve  sufficient  resolution 
to  identify  loss  peaks  due  to  excitation  of  phonons.  They  found 


a  peak  in  the  energy  loss  spectrum  in  the  region  bo tween  the  long 
wavelength  TO  and  LO  phonon  frequencies  and  identified  this  peak 
as  duo  to  excitation  of  Fuchs-Kliewer  typo  surface  optical  phonons. 
This  identification  has  been  confirmed  by  Fujivara  and  Ohtaka.^  213  ^ 
Surface  modes  of  a  different  type  have  been  observed  by  Propst 
and  Piper  who  Studied  hydrogen  and  other  molecular  species 


adsorbed  on  tungsten  using  low-energy  electrons.  Loss  peaks  were 


observed  corresponding  to  the  excitation  of  vibrational  nodes 


localized  in  the  vicinity  of  the  adsorbed  species. 


Further  extensions  of  the  inelastic  scattering  technique  have 
been  achieved  by  Ibach^  21S  ^  who  observed  surface  optical  modes 


in  ZnO  using  electrons  in  the  range  1  -  100  oV  and  resolution 


~  10  -  20  moV.  For  the  (1100)  surface,  Ibach  observed  loss  peaks 


associated  with  the  excitation  of  0,  l,  2,  3,  and  4  surface  optical 


phonons.  The  intensities  W(o-n)  of  the  peaks  are  consistent  with 


the  Poisson  distribution  predicted  by  quantum  mechanics  for  an 


oscillator  driven  by  a  tine-dependent  force  K(t) 


S(tm»)  -  ~~jG|2V^ 


<274  ) S 


G  -  i/(2n/wi)^  f*  dt  e^KUJdt 


<275  ) 


This  result  is  also  in  agreement  with  the  more  elaborate  theory  of 


Lucas  and  Sunjic/  201  ^ 
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In  n  subscquont  paper,  Ibach/  ^  studied  the  inelastic 
scattering  o t  low-onorgy  electrons  from  tho  (111)  surface  of 
silicon  and  found  a  peak  at  55  meV,  which  is  to  be  compared  to  tha 
Ra’-'an  frequency  of  65  raeV.  The  intensity  of  the  loss  peak  is 
relatively  small  compared  to  tho  case  of  ZnO,  but  is  larger  than 
one  would  expect  for  scattering  from  localized  atomic  potentials. 
Ibacb  concludes  that  his  scattering  is  due  to  a  macroscopic  electric 
field  d  wived  f r'  a  -.  iv.^ting  surface-induced  offoctivo  charges  and 
estimates  thu  charg:  to  be  on  the  order  of  0.1  e. 

As  i  jntioned  earlier,  Inelastic  neutron  scattering  is  not 
particularly  favorable  for  surface  studies  of  lattice  dynamics;, 
Uowovor ,  mon*.ion  should  be  made  of  tho  work  by  Rioder  and  Horl^  ^ 
who  studied  fine  powders  of  MgO,  Using  the  incoherent  approximation, 
they  find  deviations  of  tho  frequency  distribution  from  that  of 
bulk  crystals  which  suggest  the  presence  of  both  acoustical  surface 
modes  and  optical  surface  modes  (both  Lifshitz  *  and  Fuchs-Xliawer 
types) . 
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7.  THERMAL  DIFFUSE  SCATTER (KG 


A.  Introduction 

Thermal  diffu»j  scr.ttorlns  is  a  form  of  inelastic  scnttorir.g 
which  is  quito  familiar  from  x^ray  diffraction.  Not  only  do  the 
thermal  vibrations  of  the  atoms  lead  to  a  gain  or  loss  of  energy 
by  the  incident  boaia,  but  also  they  lead  to  a  broadening  or  dif- 
fusencss  to  the  diffraction  spots.  Ordinarily,  in  discussing 
thermal  diffuse  scattering,  one  ignores  the  change  in  energy  of 
the  incident  beam.  This  is  satisfactory,  ovon  for  lon—onergy 
electrons.  If  one  focuses  on  the  diffuceness  of  the  diffractioh 
spots  and  not  on ‘energy  loss  or  gain,  since  phonon  energies  are  . 
/«  0.01  eV  and  therefore' nsall  compared-  to  typical  incident  elec¬ 
tron  energies. 

S.  Theoretical  Development 

A  theoretical  oppression  for  tha  intensity -of  first-order 
thermal  diffuse  scattering  has  already  been  exhibited  in  2q.  (217) 
for  a  Brava is  crystal.  In  tenas  of  tha  correlation  function  de¬ 


fined  by 
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ye  can  rewrite  as 
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If  we  utilise  the  normal  coordinate  transformation  given  by  £q, 
(225  ).  then  the  correlation  function  can  be  rewritten  as 

Vi’i't  -  <VW  2  £S-sj<£>KS-&,(i')W(p)/*|'  (271 


fir  v-  ■■  ' .  >  -  *  -i^i Lj 


’-V: 

v4".  Vv ' 


-  ^§11 

*•<*'  Pi 


•I 


137. 


where  £(p)  is  the  moan  energy  of  mode  p.  Taking  the  crystal  sur¬ 
face  to  bo  parallel  to  the  1-  and  2-directions,  we  use  Eq.  (233) 
and  obtain  after  some  manipulation 
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where  wo  have  replaced  k^p'  by  £,  q  •  (ij.qg),  Q  *  (Q.  .Q 3>  j  Ms) 

~B  ^  * 

is  unity  if  its  argument  is  a  two-dimensional  reciprocal  lattice 
vector  G  and  zero  otherwise,  and 


BT(t3,t3,$«  (i/l4^S,®*^S»J*'t3)KS‘o(q.d#t./>3C5(aJ>/a2<qJ)3* 
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(280)' 


Applying  the  previously  mentioned  theorem  of  matrices,  we  can 
rewrite  ,'q)  as 

fcD*(q)| 
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where  D(q)  is  defined  in  Eq.  (234  ).  Both  Eqs.  (280  )  and  (251)' 

AT 

can  ho  conveniently  used  for  numerical  calculations.  We  note  a 

« 

rather  general  result  that  integration  of  Ig  over  and  Q3  yields 
a  quantity  which  is  a  linear  combination  of  mean  square  displace¬ 
ments  [«f.  Eq.  (235  }}.  I*  the  scattering  is  entirely  from  the 
•surface  layer,  tho  integrated  Ig  is  simply  proportional  to  the 
surface  atom  mean  square  displacement. 

The  thermal  diffuse  scattering  intensity  can  5>a  expressed  ia 
terms  of  Graen‘s  functions  by  uaihg  tlr»  relations 
. .  €<p)  -  h«d(P-n  +  £> 
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np+J  -  (u»pkBT/fi)^S^{a)p+«^>*‘1  , 


(284) 


whero  wn  -  2nnkgT/ft  with  n  an  integer.  .Then  C7U,l!)  can  be  re¬ 


written  as 


-  -  V  J  Vs  J.  Vi'£v”n> 


(285) 


-  “  *BT  Sq  QaQg  Uej»<fc*'j0)  as  T  - 

ftp 
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wheire 


}  “  <1/M>  ?  2  2 “ 


P  ttl  “ID. 


The  surface  can  be  treated  as  a  perturbation  on  an  otherwise  pcr- 


■  ;H  .■ :  ./*  • 


feet  crystal.  If  the  perfect  crystal  eigenfrequonoies  and  eigen¬ 
vectors  axe  ui(kj )  and  e(kj),  then 


rj_,a,4'j»2)  -  u/nvo  s  s  e*(kj)eB<k'jO  x 

03  kJ  Of  -  P  ~ 

—«  /V 


ik.HU)+ik'.R'(4) 

^  /V  A/  <V 


u(k.1,lc'j‘',B^) 

A#  A/ 


whore  uCkj  ,k'j m2)  is  given  by  Eq.  (206). 

A*  #«* 


Wallis  ard  Maradudin'222'  have  evaluated  Ig  with  the  aid  of 
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Eqs,  (277)  and  (285  )  for  tho  Maradudin-Velngailis  isotropic  tsodel'* 


•of  a  sc  crystal.  Both  high  temperatures  and  long  wavelengths  wore 
assumed.  For  scattering  only  fro::  tbs  ‘surface  layer  of  a toes ,  the 
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2  2  2 

whore  q  r  f  q^,  a  is  the  lattice  constant, 

b  -  I  «  (knT/Mc2)c  fV(0)e-2«(0) 
o  s  a  2  o  o 

and  c0  is  the  speed  of  long  wavelength  transverse  acoustic 


:290> 


Note  that  the  Intensity  is  proportional  to  the  nunber  of  surface 

atoms  and  to  the  absolute  temperature.  An  interesting  feature  is 

the  inverse  proportionality  of  intensity  to  the  distance  q  of  tho 

scattering  vector  from  the  nearest  reciprocal  lattice  rod.  t  Bulk 

scattering,  in  contrast,  varies  as  1/q  .  Wallis  and  Maradudin 

also  found  that  the  iso->£nt.ensity  curves  for  the  surface  and  bulk 

cases  are  frequently  different.  This  is  illustrated  in  Fig.  17. 

Another  calculation  based  on  the  .isotropic  model  is  that  of 
(223) 

Huber  who  considered  a  slab  of  a  general  isotropic  material 

,  {98) 

following  -the  approach  of  Stratton.  •  Huber's  result  for  tho 
intensity  of  thermal  diffuse  scattering  can  be  written  as 


I2  -  ^  E'rlQi-GfclrtlsJ  cos  2«?  +  VqJ  sin 
Gg 


vah 


(291) 


where  the  scattering  vector  Q  a  K'-Z  -  (QB ,Q  ),  U,  V,  apd  W  arc 

fv  »v  r*  rwH  J, 


0 

1 

functions  < 

1  . 

waves,  p  i 

0 

Jr  ; 
s  5 1 

8iIoV^* 

minimum  value  of  ISj-Gj I  is  to  be  counted.  Qualitatively,  Huber's 
result  is  eiwilar  to  that  of  Wallis  and  Maradudin,  but  quantita¬ 
tively  It  is  larger  by  about  a  factor  of  2.  The  '-yigin  of  this 
discrepancy  has  hot  been  pinned  down. 

A  numerical  calculation'  of  the  thermal  diffuse  scattering 
from  the  (10G),  (110),  and  (111)  surfaces  of  argon  kae  been  carried 
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out  by  DeWette  ond  Allen 


(224) 


using  a  henna rd-Jones  potential, 


Tho  oigenfrequencies  and  eigenvectors  for  a  slab  11  atonic  layers 


thick  wore  calculated  and  then  used  in  Eq.  (  280  )  to  obtain  lg. 


Tho  results  are  not  valid  in  tho  immediate  vicinity  of  the 


reciprocal  lattice  rods,  but  give  a  better  picture  of  tho  situation 


far  from  the  reciprocal  lattice  rods  than  does  the  continuum  ap¬ 


proximation.  These  calculations  have  recently  been  extended  by 


(225  5 

Xesmodel,  DeWette,  and  Allen  who  have  treated  in  detail  the 


(100) ,  (110) ,  and  (111)  surfaces  of  Xe.  They  included  contribu¬ 
tions  to  the  scattering  from  tho  first- five  layers  of  the  nine- 


layer  slabs  employed  and  assumed  a  geometrical  decrease  in  the  trans¬ 


mission  factor  with  distance  from  the  surface.  Kesmodel  et  al.  found 


significant  interference  effects  between  layers  leading  to  different 


intensity  distributions  around  different  reciprocal  lattice  rods. 


Thoy  also  found  an  ''extra"  spot’ near  the  (04)  rod  for  the  (100) 


surface. 


C.  Experimental  Results 

Early  observations  of  surface  thermal  diffuse* scattering  were 
made  by  Aldag  and  Stern ^  155  ^  who  studied  the  (110)  surface  of 


tungsten  with  electrons  in  the  230  -  700  eV  range.  They  observed 
an  anisotropy  in  tho  scattering  vhich  thoy  attributed  to  aniso-  . 


tropic  surface  modes. 

A  detailed  study  of  thermal  diffuse  scattering  from  the  (111) 

(  226) 

surface  of  silver  has  been  reported  by  RcKinnoy,  Jones,  and  ^obb; 


Those  authors  verified  c  number  of  the  theoretical  predictions  for 


surface  thermal  diffuse  scattering.  Thus,  they  found  tho  predicted 


inverse  dependence  of  intensity  on,  distance  from  tho  nearest  recip¬ 


rocal  lattice  rod.  Furthermore,  the  ratio  of  diffuse  intensity 
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absolute  toaporaturo  and 


to  p6ak  intensity  is  proportional  to  tho 
to  tho  square  of  tho  magnitude  9x  tho  scattering  vector.  Finally, 
tho  ratio  of  integrated  diffuse  and  peak  intensities  is  approxi¬ 
mately  proportional  to  the  exponent  in  the  Dobye-Waller  factor. 

In  a  subsequent  paper,  Barnes,  Lagally,  end  ffebb(227)  dis~ 
cuesod  experiments  using  tho  (111)  surface  of  nickol.  They  ob¬ 
served,  in  addition  to  the  one-phonon  scattering,  a  uniform  back¬ 
ground  which  is  a  function  of  energy  and  temperature.  Barne?  et  al 
present  Arguments  that  this  background  is  due  to  multiphonon 
scattering  -  i.e.,  to  higher  order  terms  in  the  expansion  in  Eq. 
(213).  In  the  6ase  of  tho  (111)  surface  of  silver,  tho  multi¬ 
phonon  scattering  is  relatively  small.  Lagally,  Ngoc,  and  W228> 
have  'extended  this  work  and  have  show  that  kinematic  intensities 
can  bo  recovered  from  LBED  data  by  averaging  tho  LEBd  intensities  ‘ 
at  constant  momentum  transfer.  This  procedure  could  prove  useful 

ia  determinations  of  surface  structures  and  surface  Debyo-Waller 
factors* 
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8.  SURFACE  EFFECTS  ON  OPTICAL  PROPERTIES 


In  this  section  we  shall  discuss  the  effect  of  freo  surfaces 


on  optical  proportios  associated  with  lattice  vibrations,  In 
particular,  we  shall  discuss  infrared  absorption  or  reflection 


and  Raaan  scattering. 


A.  Infrared  Absorption 


(i).  Theoretical  Results 


Since  infrared  wavelengths  are  large  compared  to  lattice 


spacings  in  crystals  one  can  treat  the  wave  vector  of  the  infrarod 


adiation  as  negligible.  Undor  these  circumstances,  the  infrared 
absorption  coefficient  of  a  crystal  lattice  is  proportional  to 
tho  quantity*229* 


Uw)  -SEv^J  tf(R)H<R)?i(R)dR|2  x 


x  &(£j  ~  Ej  +  ha) 


(292) 


where  R  stands  for  the  set -of  miclehr  coordinates,  f*  and  aro 
the  vibrational  wavo  functions  for  the  initial  and  final  states. 


and'Ef  aro  the  corresponding  energy  eigenvalues.,  is  the 
Bolttaann  occupancy  factor,  and  M(B)  is  the  electric  dipole  aojsont. 

Tho  dipolo  ssosent  can  bo  expanded  in  power  series  in  the 
displacements  of  tho  ions  frea  equilibriua  ti(ty) : 
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«Tha  constant  term  «i  is  of  no  interest  for  infrared  absorption, 
oa 

Tho  terras  linear  and  quadratic  in  the  displacements  are,  rospcc- 
*  »  * 
tivoly,  the  first  order  and  second  order  terms,  tfe  shall  restrict 

our  discussion  to  the  first  order  terms,  which  are  important  in 

ionic  crystals  and  are  subject  to  strict  selection  rules  if  the 

crystal  is  poricdic. 

Th3  coefficients  H  „(&t>  play  the  role  of  effective  ionic 
ep  ~ 

charges.  Infinitesimal  translational  invariance  imposes  tho 

i 

•  condition 

£  U„aUx)  “  0  ..  (294) 

Ik  ?  ~  . 


In  a  periodic  crystal,  ll  (tx)  does  not  depend  on  l,  bo 

Qt  p 


S  H  (x)  -  0  •  (295) 

x  a* 

It)  many  cassc,  tho  tensor  li  Ax)  can  be  taken  to  be  diagonal 

crp 

V*5  “  V*  ■  <29£)  • 

where  *  is  a  scaler  effective  charge  for  ions  for  type  x.  Again 

4  # 

cons.tderi&g  a  periodic  crystal,  one  can  expand  in  terms  of  tho 
noraiil  coordinates  Q(kj)  instead  of  the  displacements  uUx) 


“«<S>  -*<»*£  ya>«;y>  +  •• 

S3 


whore 


-1  ik*R(o)U) 

-  WO  5  B  fxJjjJ)*,  ~~  ~ 


(297) 


Q<kj>  (298) 


and  f^(xjkj}  is  specified  by  Eq.  (143).  It  can  be  shown <3) 
that  the  coefficient  II  (kj)  vanishes  unless  k  •*  0  and  J  refers  to 


'v**  -t*,^ 


144. 


the  transvorso  optical  branch.'  This  result  is  a  consequence  of 
periodicity;  when  a  free  surface  is  present  other  modes  can  give 
a  non-vanishing  dfipolo  monent  and  infrared  absorption. 

A  theoretical  investigation  of  the  effect  of  free  surfaces 

on  tho  infrared  absorption  of  an  ionic  crystal  was  carried  out  by 

/  230  % 

Rosenstock'  '  who  simplified  the  calculations  by  retaining 

only  short-range  forces  in  tho  lattice  dynamics.  Tho  first  order 
dipole  monent  was  evaluated  with  the  aid  of  Eqs.  (293), 

(  295  ),  end  <296  ).  Tho  use  of  Eqs.  <  295  )  and  (296) 

is  not  strictly  correct,  but  should  bo  all  right  for  qualitative 
estimates.  For  a  linear  chain  of  N  ions  of  alternating  charge  but 
equal  mass  and  nearest  neighbor  interactions,  Rosenstock  found  for 
tho  first-order  dipole  aohent  associated  with  normal  code  s, 

Mg  •"  ouo/cos(<ps/2)  ,  s  odd  .  *  (299) 

“0  ,  s  even 

where  <p8  •  sn/K,  s  -  0,  1,  2,  3,  . .K-l,  and  uQ  is  an  arbitrary 
anplitude.  Tho  dipole  jeoreont  has  a  raaximun  for  s  *•  ?f-l  (the  mode 
which  has  tho  naxinun  frequency  wn  and  which  would  bo  the  long-wave 
length  optical  mode  if  tho  ‘lattice  were  made  diatomic)  and  decreases 
in  magnitude  as  s  decreases.  The  important  point  is  th2t  the  dipolo 
Btoaent.,  and  hance  tho  infrared  absorption,  are  non-zero  for  s  *  N-l. 
Thus  tho  froe  surfaces  load  to  absorption  in  frequency  regions  not 
allowed  with  periodic  boundary  conditions.  Taking  into  account 
the  frequency  distribution,  tho  absorption  coefficient  turns  out 
to  bo  proportional,  to  {»_  -  w)“3/2. 


-r*^-*^®*  <**.  vr~«r-»f  /«-»->  ^wlN. . ,«  ,^_  ^  ^ 
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•  |,  * 

*  ,  la  o  second  paper,  Rosonstock*  231  *  treated  two-  and 

throo-diiaonsionni  crystals.  The  Coulomb  interactions  wore  ignored 
•  in  tho  lattice  dynamics  by  ubc  of  tho  Rosonstock-Howoll  modol*41^ 

The  main  difference  from  tho  one-dimensional  coho''  is  tho  appoarnnee 
of  critical  points  in  the  frequency  distribution.  Since  the  free 
surfaces  lead  to  a  continuous  distribution  of  infrared  absorption 
over  the  normal  modos,  tho  critical  points  produce  subsidiary  peaks 
superposed  on  a  continuous  background  absorption.  This  is  in  sharp 
contrast  with  tho  single  delta-function  absorption  predicted  by 
poriodic  boundary  conditions, 

Rosenstock  has  subsequently  included  tho  Coulomb  interactions 
in *hi8  analysis.  For  tho  one-dimensional  casof232  ^  he  finds  that 
tho  absorption  coefficient  varies  as  t-log^-w)  ]^(u»  -u>)“3/2  rather 
than  as  (u^-ui)  without  Coulo&b  interactions.  For  three- 
dimensions,  *  ho  has  calculated  tho  frequencies  of  subsidiary 
peaks  using  a  model  of  NaCl-type  lattices  involving  Coulomb  inter¬ 
actions  and  nearest~n8ighbor  repulsive  interactions.  Rosenstock 
found  that  the  experimental  values  of  the  fundamental  lattice  ab¬ 
sorption  frequency  for  many  alkali  halides  av3  fitted  better  by  a 
subsidiary  frequency  than  by  tho  classical  limiting  frequency.  How¬ 
ever,  it  is  not  clear  that  tho  model  used  is  sufficiently  precise 
to  removo  all  doubt  about  such  an  identification. 

Recently.  (Bass  and  Rosenstock* 234-235>  havo  taken  into 
account  the  difforemt  masses  which  the  ions  of  different  charges 
generally  have  in  an  ionic  crystal.  They  furthermore  investigated 
,  dependence  of  the  infrared  absorption  in  various  frequency 

regions  on  the  volume  and  surface  area  of  the  crystal.  ThoF  found 
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that  the  absorption  in  a  frequency  interval  containing  tho  funda- 
montal  lattice  frequency  is  proportional  to  the  volume,  but  tho 
absorption  in  othor  frequency  regions  is 'proportional  to  tho  surfaco 
area.  Thus,  to  experimentally  observe  tho  absorption  induced  by 
the  free  surfaces  away  from  the  fundamental,  one  should  do  experi¬ 
ments  on  samples  consisting  of  many  thin  layers  to  enhance  the 
surface-to-voluoo  ratio.  Tho  experiments  should  also  be  dono  at 
low  temperatures  to  avoid  complications  due  to  anhariaonic  broaden¬ 
ing  of  the  fundamental  absorption. 

Since  tho  fundamental  absorption  of  infrared  radiation 
by  a  crystal  is  primarily  restricted  to  small  wave  vectors,  it  is 
reasonable  to  apply  macroscopic  theory  in  tho  frequency  region 
near  tho  fundamental.  This  has  boon  done  by  Englman  and  Ruppln^7^ 
who  have  discussod  tho  optical  properties  of  finite  ionic  crystals. 
Those  authors  have  recently  published  a  review^  71  \'of  the 
situation,  r.o  shall  take  up  only  selected  topics  and  refer  tho 
reader  to  tho  Kuppin-Englman  review  article  for  further  details. 

Vo  now  consider  the  optical  proparties  of  finite  crystals. 
For  spheres,  cne  can  utilize  tho  Hio  theory $  236  ^  Tho  absorption 
cross  section  (in  units  of  nil-,  R  tho  radius)  is  given  by 

08  “  ofm1  m+1)(~Bo  V}*JS  ~  VIbJ2)  (300) 
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For  small  but  finite  R,  tho  l  »  1  modo  frequency  has  corrections 
specified  by  tho  eqvtations 


4>a)  .  go  *  2eM(1^ 

uT  €w  +  2€y(l+b) 


(306) 


<307) 


In  Fig.  18  ,  the  extinction  cross  section  calculated  for  a  1  pm 

radius  sphere  of  KaCl  is  plotted  as  a  function  of  frequency.  Tho 
peak  due  to  tho  A  •*  1  mode  is  clearly  evident  between  the  long  wave¬ 
length  transverse  and  longitudinal  optical  mode-  frequencies.  As 
the  size  of  the  sphere  increases,  additional  peaks  appear  at  tho 
l  -  2,  3, ....  surface)  modo  frequencies  and  at  bulk  modo  frequencies 
lying  Just  below  Wj,. 

Englnan  and  Ruppin  have  also  considered  absorption  by 
cylinders.  For  very  narrow  cylinders  (radius  ~  Ipa)  and  electric 
vector  perpendicular  to  the  axis  of  the  cylinder,  there  is  a  peak 
in  extinction  at  ,a  surfaco  modo  with  frequoncy  specified  by 


“so  •“?<eo+€«>/<^+«I! '?  ♦  <3CS) 

For  parallel  polarization,  the  extinction  peaks  at  u;^.  A  related 
problem,  tho  absorption  by  needle-shaped  crystals,  has  been  dis¬ 
cussed  by  Rosenstock. ^237^ 

A  particularly  interesting  case  is  the  thin  slab.  Ws  have 
aeon  in  Section  2  that  a  slab  exhibits  surface  modes  of  the  Fuchs- 
Kliewer  type  with  frequoncy  specified  by  Eq.  <  S3  ).  These  nodes. 


however,  aro  rionradiativa  and  do  not  interact  directly  with  elec¬ 
tromagnetic  radiation,  although  coupling  can  bo  achieved  by 
special  meano  to  be  described  later.  The  optical  properties  of  a 
slab  of  an  alkali  halide  have  been  discussed  by  Berrenan. '  ' 

For  radiation  incidsut  normally  on  the  slab,  there  is  the  usual 
absorption  peak  at  m-j,.  For  non-normal  incidence,  Berreman  ob¬ 
served  an  additional  absorption  peak  at  ui^  if  P-polarization  is 
used  (electric  vector  parallel  to  the  plane  of  incidence).  That 
coupling  should  occur  in  this  situation  is  clear,  since  the  elec¬ 
tric  vector  has  a  component  normal  to  the  surface  and  therefore 
parallel  to  the  macroscopic  electric  fi?ld  of  the  longitudinal 
optical  phonon.  A  detailed  theoretical  analysis  of  the  optical 
properties  of  an  ionic  slab  ha3  been  given  by  Fuchs,  Kliower,  and 
Pardee^^  using  a  macroscopic  approach  and  by  Jones  and  Fuchs^®5  ^ 
uaing  a  microscopic  approach. 

In  a  second  paper,  Berrenan^240^  analyzed  the  effects  of 
microscopic  pits  and  domes  on  tho  infrared  reflectance.  He  found 
that  these  defects  introduce  structure  into  the  reflectivity  .just 
boic<» 

In  tho  foregoing  discussion,  the  boundary  conditions  have 
been  handled  satisfactorily  either  for  the  short-range  forces 
(Bosenstocb)  or  tho  long-range  forces  but  not  for  both.  Tong  and 
Haradudin^  30^  and  Grimm,  ilaradudin,  and  Tong*-24**  have  tackled 
the  problem  of  handling  tho  boundary  conditions  with  both  types  of 
forces.  The  imaginary  part  of  the  dielectric  response  tensor 
(rhich  relates  the  polarization  to  the  g stomal  electric  field) 
can  b®  writ tea  a? 


•VS-—'-  -m  jW“W?  T.:*;;  •  ‘.f ‘ 


C3 

€®>(.u)  -  (2n/ftn}(opft*-l>  f  dfc  e“iwt(M  ,  (t)M  (0)>  <309  ) 
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whsre  n  is  tho  volume  of  the  crystal,  0  ~  1/kgT,  M^(t)  is  given  by 
£q.  (293  )  with  the  displacements  u^(tx)  taken  to  be  time  dependent, 
and  the  angular  brackets  denote  an  average  over  a  canonical  ensemble. 
By  generalising  the  normal  coordinate  transformation  specified  by 
Eqs.  (2S0  )  and  (251  )  to  two  atoms  per  unit  cell  and.  sriting  tae 
time  dependence  of  the  creation  and  annihilation  operators  ar 


a^(kj|,t)  “  aj{kii,0)exp£“itu.i(ka)t3 


(310  > 


a^f  (~k5  ,t)  -  a^  (*-k5  ,0)exp[.i^  (kj)tJ  , 
one  can  roduce  the  expression  for  €jv  («)  to  the  fora 
<21  nV  *E5W"U>i'P)3  * 
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whoro 
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see  that  only  k^-O  modes  contribute  to  ^v'(u>)  and  hence  to  tho 
absorption.  The  tetragonal  sysaaotry  of  the  slab  is  reflected  in 
^(e)  by  the  fact  that  fa)  +  €^5(o). 

.These  calculations  have  been  extended  by  Grimm,  Maradudin, 
and  Tong^^  who  included  the  real  part  of  €^v<sb).  For  a  100-layer 
slab  of  KaCl,  they  find  that  the  diagonal  components  can  be  repre¬ 


sented  by 
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where  ui0  *  10  sec  and  y  Is  the  damping  constant  (.-0.02tao  at 
roots  temperature).  The  poles  of  £  (ta)  In  the  limit  y*o  are  at 
oi~  “  2.491ui  and  at  a  transverse  surface  mode  (s~.„  «  2.41R»  .  This 
surface  mode  is  localized  within  a  few  atomic  layers  of  the  sur¬ 
face  and  is  distinct  from  the  Fuchs-Kllewer  mode.  The  zero  of 
ezz(w)  in  the  limit  r*o  is  at  u,'L  **  5.837w0.  For  an  anglo  of  inci¬ 
dence  of  30°,  the  absorption  coefficient  is  plotted  against  fre¬ 
quency  in  Fig.  19.  Peaks  are  clearly  evident  at  «TSjj»  s»T,  and 
As  mentioned  previously,  there  is  no  direct  coupling  of  the  radia¬ 
tion  tc  the  Fuchs- Kliewor  mode  and  hence  no  peak  in  absorption  at 
its  frequency. 

<iil  Experimental  Results 

A  great  many  workers  have  made  infrared  absorption 
measurements  of  very  f ino  powders  of  various  materials.  He  con¬ 
sider  only  a  few  examples  and  refer  the  reader  to  the  review 

1 7i ■) 

article  of  Ruppin  and  Englraan  7  for  further  references.  Typi¬ 
cally,  the  powders  are  suspended  in  ainqral  oil,  polyithylene,  or 
KBr  disks  and  are  found  to  give  peak  absorption  at  a  frequency 
higher  than  <s_,.  Assuming  that  the  particles  can  -be  approximated 
by  spheres,  these  results  are  entirely  consistent  with  the 
Englaan-Ruppin  result  that  the  peak  absorption  of  very  small 
spheres  occurs  at  the  surface  mode  frequency  given  by  Eq.  (89) 

„  rather  than  at  bj^.  Thus,  for  CaO;  MeBevitt  and  Baun^2*2^  *  ob¬ 
served  peaks  at  400  cbT*  and  290  ceT",  These  are  close  to  the 
calculated  surface  mode  frequency  of  402  cs*1  and  the  bulk 


transverse  frequency  of  289  cm  ,  Sira liar  results  have  boon  ob¬ 
tained  for  NiO,  MgOf  and  a  variety  of  n ore  complicated  crystals. 

Martin  has  studied  fine  powders  of  KCl^*53^  and  KaCl^344^ 


:  ‘ :.-.!  r 
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and  observed  absorption  in  the  region  between  asT  and  su^.  The 
peak  in  absorption  occurs  at  a  somewhat  lower  frequency  than  &y. 
This  may  be  duo  to  absorption  by  £>1  nodes  or  to  the  r.on-spherical 
shape  of  the  crystals.  Genzel  and  Mart  in  have  recently  ten  no 
infrared  measurements  on  finely  powdered  CdO  containing  free  car¬ 


f  . 
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riers.  One  may  therefore  expect  coupled  surface  optical  phor.on- 

w  • 

surface  plasaon  modes.  Tho  experiments  indeed  reveal  two  reso¬ 
nances  in  the  absorption  coefficient  as  well  as  an  antlresonancq 
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which  is  also  predicted  by  theory. 

Systems  of  particles  which  can  be  approximated  by  cyl¬ 
inders  rather  than  spheres  have  also  been  investigated  experiraent- 
.  ally.  Pultz  and  Hertl^^  studied  very  thin  fibres  of  p-SiC  and 

observed  peaks  at  av,  and  at  the  surface  mode  as  predicted  by  the 

(2  *7) 

Englraan-Ruppin  theory.  Hass'*'  *  carried  out  an  interesting  study 
of  tho  transmittance  of  thin  films  of  NaCl  and  of  I»iF  deposited 
on  diffraction  gratings.  For  radiation  with  electric  vector 
polarized  parallel  to  the  grooves  the  transmittance  had  a  minimum 
at  cT  whereas  with  polarization  perpendicular  to  the  gjeoves,  tho 
minimum  moved  to  a  higher  frequency  whose  value  corresponds  to 
tho  surface  mode  in  an  elliptic  cylinder  with  axial  ratio  of  30. 

Ke  have  alluded  to  Berreman’s  work  oik  the  effect  of  pits 


py  ■  •;  • 


and  domes  on  the  optical  properties  of  a  crystal.  Another  typo  of 
surface  defect  which  must  bo  considered  is  an  impurity  atom  ad¬ 
sorbed  oa  the  surface.  Some  years  ago  Plisfcin  and  Eischens^8^ 
studied  tho  infrared  absorption  associated  with  hydjrogea  adsorbed 
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on  plhtinura.  Th sy  observed  two  absorption  peaks  which  they 

attributed  to  two  different  types  of  bonding.  The  vibrational 

origin  of  th8  peaks  was  confirmed  by  the  isotopic  shift  when  hy- 

di  jgcn  was  replaced  by  deuterium.  In  another  paper,  Eischens  and 
(243) 

Pliskin  investigated  carbon  monoxide  adsorbed  on  platinum  and 
found  two  absorption  bands  at  476  cm-i  and  2000  cs~*.  If  it  is 
assumed  that  the  CO  molecule  is  bound  to  a  surface  Pt  atom  by  at 
C-Pt  bond,  then  the  force  constants  for  the  C-Pt  and  C-0  bonds 

I 

/Cg\ 

can  be  determined  using  the  theoretical  analyses  of  Grimley'  ' 

(38) 

and  of  Masri  and  Armand. 

It  has  been  pointed  out  that  radiation  doos  not  interact 
directly  with  the  Fuchs-Kliewer  surface  inodes  of  an  Isolated  ionic- 
slab  with  perfectly  planar  surfaces.  This  is  a  result  of  the  dis¬ 
persion  curve  for  the  Fuehs-Kliewor  modes  lying  to  the  right  of 
the  light  lino  ui-kc.  However,  by  suitable  experimental  tricks, 
it  is  possible  to  make  those  modes  accessible  to  experimental 
study. 

One  procedure  is  to  study  the  infrared  reflectance  from 
a  sample  upon  which  a  grating  has  been  ruled..  If  the  spacing  of 
the  grooves  is  d,  then  the  wave  vector  component  of  the  radiation 
parallel  to  the  surface  is  augmented  by  2nn/d  where  n  is  an 
integer.  ThU3 , 

kj  -  (w/c)  sin. 9  +  2nn/d  (316) 

where  6  is  the  angle  of  incidence.  The  reflectance  curve  exhibits 
dip's  at  frequencies  specified  by  the  surface  polariton  dispersion 
eurvo  and  Eq.  (316  ).  From  the  observed  reflectance  dips,  an 


experimental  dispersion  curve  can  bo  deduced.  This  technique  has 
been  employed  by  Uarscball,  Fischer  and  Queisser^23^  to  determine 
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the  dispersion  curve  for  surface  plasaons  in  n-InSb  and  by  Ander¬ 
son,  Alexander,  and  Boll^5^  to  investigate  the  interaction  of 
surface  piasmons  with  optical  phonons,  also  in  n-InSb. 

An  alternative  and  perhaps  more  convenient  procedure 

£252  (252 

is  attenuated  total  reflection.  Ruppin  suggested  that 

this  technique  would  be  useful  in  the  experimental  investigation 
of  nonradiative  surface  polaritons.  The  incrementation  of  the 
wave  vector  is  achieved  by  placing  a  prism  adjacent  to  the  sample 
.  of  interest  and  separated  from  the  sample  by  a  few  microns.  In 
the  absence  of  the  sample,  tho  radiation  is  totally  internally 
reflected  in  the  prism.  In  the  presence  of  tho  sample,  the  ex¬ 
ponential  tail  of  radiation  penetrating  outside  the  sample  couples 
with  the  surface  polaritons  in  the  sample.  Coupling  is  possible 
because  the  effective  wave  yector  of  the  radiation  parallel  to  the 
surface  is  given  by 

tcj  “  (<a/c)  n  sin  o  (317  ) 

where  n  is  tho  refractive  index  of  the  prism  and  a  is  tho  angle  of 
incidence  in  the  prism.  Since  n>l,  wave  vectors  to  tho  right  of 
the  light  line  in  vacuum  can  be  achieved.  This  technique  has  re¬ 
cently  been  exploited  by  Harschall  and  Fischev^25^  who  determined 
the  dispersion  curve  for  surface  optical  phonons  in  CaP.  An 
independent  investigation  ‘~ang  tho  satt&  lines  is  that  of  Sryksin, 
Gerbshtoln,  and  Mlrlin*255*  on  KaCl  and  CaFg. 
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B.  Raman  Scattering 

In  Raman  scattering,  tho  frequency  of  tho  scattered  radiation 
is  down-shifted  (Stokes)  or  up-shifted  (anti-Stokes)  from  t'bc 


155. 

frequency  of  the  incident  radiation.  The  advent  of  high  intensity 
lasers  has  stimulated  ar.  enormous,  upsurge  of  experimental  investi¬ 
gations  of  tho  Raman  effect  in  hulk  crystals.  However,  the  use 
of  the  Raman  effect  to  study  surface  phonons  is  in  a  very  rudi¬ 
mentary  3tate. 

Raman  scattering  can  be  ascribed  to  the  modulation  of  the 
electronic  polarizability  by  the  vibrations  of  the  ions.  Tho  in¬ 
tensity  of  Raman  scattering  at  frequency  ta  is  proportional  tb  the 
fai 

quantity v ' 

i«Y.»^<“>  "5/i<V,lPoYivXVlP3-xlv/>5  6(w-®0-*VV'>  (318) 
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where  v,  v *  designate  the  vibrational  quantum  numbers  of  the' crys¬ 
tal,  u)0  is  the  incident  frequency,  ftai  i  is  the  change  in  vibra¬ 
tional  onorgy  of  the  crystal  and  P  is  tho  electronic  polariz¬ 
ability  tensor.  The  polarizability  can  be  expanded  in  power  series 
in  the  ionic  displacements.  Kaon  one  transforms  to  normal  coor¬ 
dinates  Q(kj),  one  finds  for  a  bulk  (poriodic.)  crystal  that  only 
optical  modes  with  k«Q  are  active  in  tho  Raman  effect.  The  polar¬ 
izability  expansions  enn  then  bo  writ ion  an 

P«g  -pSJ  +SPo3(j)Q(oj)  +  ...  (319) 

w 

where  j  runs  over  the  optical  branches.  It  turns  out  from  symmetry 
considerations  that  crystals  such  as  NaCl  with  every  atom  at  a 
center  of  inversion  do  not  exhibit  a  first-order  Raman  effect. 

However,  crystals  such  as  diamond  which  have  centers  of  inversion 
midway  betwoen  pairs  of  atoms  do  show  a  first-order  Raman  effect. 

Tho  creation  of  a  free  surface  can  result  in  the  elimination 
of  centers  of  inversion  symmetry  and  lead  to  surface-induced  Raman 


scattering,  A  theoretical  discussion  has  been  given  by  Ruppin 
and  Engloan^56^  who  treated  tho  ease  of  tho  incident  beam  perpen¬ 
dicular  to  a  thin  slab  and  tho  scattered  beam,  at  an  angle  6  to 
the.  normal.  Conservation  of  energy  and  momentum  components  parallel 
to  the  surface  yield  the  conditions  for  Stokes  scattering 

Vj,  “  +  ®  (320) 

“  kja  -  jkgjsin  8  (321) 

where  the  subscripts  1  and  2  refer  to  the  incident  and  scattered 

radiation  and  the  unsubscripud  quantities  refer  to  the  phonon. 

From  Eqs.  (320  )  and  (321  ),  one  finds 

cqB 

8ln  6  *  5^5lqi7  ♦  {322> 

This  relation  between  qB  and  6  can  bo  used  in  principle  to  deter¬ 
mine  the  dispersion  curve  iu(qj).  . 

According  to  Rupp in  and  Englnan,  tho  Fuchs-Klicwcr  modes  of  a 
-lab  of  NaCl-type  crystal  will  not  exhibit  first-order  Raman 
scattering,  since  tho  displacement  patterns  of  thqso  surface  modes 
are  similar  to  those  of  tho  bulk  longitudinal  and  transverso  opti¬ 
cal  modes  which  are  Raman  Inactive.  For  crystals  which  are  Raman 
active  in  the  bulk,  however,  the  Fucbs-Xllewer  surface  modes  should 
be  Raman  active. 

A  calculation  which  emphasizes. the  effect  of  tho  surface  on 

(2571 

the  short-range  interactions  is  that  of  Litzman  and  Cely  'who 

analyzed  the  Rosonstock-Kewell  model  of  a  NaCl-type  crystal  with 

•  {258} 

a  (OGi)  surface.  Following  tho  approach  of  Loudon,'  .  Litzman 
and  Cely  show  that  the  transition  probability  of  first  order  Raman 
scattering  into  unit  solid  angle  can  be  expressed  as 
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where  the  subscripts  o  and  s  refer  to  the  incident  and.  scattered 

I 

radiation,  o  is  the  polarization  vector  of  the  radiation, 
is  a  third-order  electronic  matrix  element,  n  is  a  Bose-Einstein 
factor,  the- lower  and  upper  lines  refer  to  Stokes  and  antiStokos 
scattering,  and  UoX (4x,4 'x* ,ujZ)  is  the  Green's  function  specified 
by  Eq.  (  136).  After  evaluating  the  Green's  function  for  a  mon¬ 
atomic  Rosenstock-Xewell  model,  Litznan  and  Cely  find  a  continuous  . 
band  of  scattering  which  peaks  at  a  frequency  below  the  maximum 

frequency  of  the  lattlco.  Unfortunately,  no  surface  nodes  exist 

fjsg) 

for  the  model  considered.  In  a  recent  publication,  Kawabatav 

>  * 

has  given  a  theoretical  treatment  of  light  scattering  by  surface 
phonons  (Rayleigh  waves)  in  metals. 

C.  Brillouin  Scattering 

The  inelastic  scattering  of  light  associated  with  the  emission 
or  absorption  of  acoustic  phonons  is  termed  Brillouin  scattering 
and  can  be  viewed  as  arising  from  the  Doppler  effect.  The  frequency 
shift  is  where  a>a  is  the  acoustic. mode  frequency.  Experimentally, 
one  can  drive  a  surface  acoustic  wav®  with  an  external  power  source. 
Tho  acoustic  wav©  then  serves  as  a  diffraction  grating  so  that  dif¬ 
fracted  beans  appear  in  too  scattered  radiation.  Tor  the  at21  order, 
sin  @a  -  sic  60  +  (a X/A)  (324) 
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§.  SURFACE  THERMAL  EXPANSION 

tfe  have  seen  in  Section  2  that  the  equilibrium  spacings  of 
tho  atomic  layers  close  to  a  free  surface  are  in  general  <1  if  lor- 
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ent  from  the  bulk  value.  One  may  also  anticipate  that  the  chango 
in  equilibrium  spacing  with  temperature— i.e. , -the  thermal  expan¬ 
sion-will  also  be  different  near  a  surface  from  that  in  the  bulk. 

Experimental  information  concerning  surface  thermal  expansion 
can  be  obtained  ih  principle  from  the  temperature  shifts  of  i^ragg 


peaks  in  low  energy  electron  diffraction  if  tbo  kinematic  approxi¬ 
mation  is  valid.  Multiple  scattering  effects  will  tend  to  obscure 
the  extraction  of  surface  thermal  expansion  values  from  the  experi¬ 
mental  data. 
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A  discussion  of  experimental' and  theoretical  aspects  of  surface 
thermal  expansion  has  been  g Iven  by  Wilson  and  Bastow. ^260^ ,  The 
positions  of  normal  incidence  diffraction  peaks  are  given  in  the 


kinematic  approximation  by 
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(326) 


where  Q  is  the  scattering  vector,  1  is  the  electron  wavelength,  n 
is  an  integer,  and  a.  is  the  spacing  between  equivalent  layers  in 
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tho  direction  normal  to  tho  surface.  The  electron  accelerating 
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potential  Bp  and  the  inner  potential  are  related  to  X  by 


where  E *  “  B_  +  V. . 
p  p  i 


UA)  ~  [150.4/Ep.(eV}}* 

The  thermal  expansion  is  then  given 


(327) 
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For  low  values  of  E  ,  the  scattering  will  bo  mainly  from  the 

r 

surface  and  tho  surface  thermal  expansion  will  be  obtained. 
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Experimental  results  have  been  reported  by  various  workers, 
Celatt,  Legally,  and  Webb^261^  found  values  of  «s/«b  («*b  is  the 
bulk  value)  on  the  order  of  one  to 'two  for  Ag  and  Ni.  However, 
temperature-dependent  asymmetries  of  the  Bragg  peaks  due  10 
absorption  and  tha  variation  of  the  Debyu- Waller  factor  with  dis¬ 
tance  from  the  surface  complicate  the  interpretation  of  the  data. 
Wilson  and  Bastov/2®®^  have  studied  the  (100)  surfaces  of  Cr  and 
Mo.  They  find  values  of  0^/0^  in  the  range  between  two  and  three. 
Ignatievs ^262*  has  investigated  the  (111)  surface  of  Xe  and  found 
Og/ofjj  valuos  of  about  four  or  five. 

A  simple  theoretical  treatment- of  surface,  thermal  expansion 
has  been  given  by  Wilson  and  Bastow^2®*^  based  on  the  Griineiscn 
expression  for  thermal  expansion.  They  find 


o  -  (kjjY/a^Xft/kgej))4 


(329) 


where  y  is  the  Grunoiscn  constant,  W  is  the  atomic  mass,  aQ  is 
the  lattice  spacing,  and  9D  is  the  Debye,  temperature.  If  y  is  the 
same  for  surface  layers  as  for  bulk  layers , . then 
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This  relation  seems  to  be  in  rough  agreement  with  the  available 
experimental  data. 

In  a  more  complete  theory,  one  would  minimize  tho  Helmholtz 
free  energy  with  respect  to  variations  in  the  various  layer 
spacings.  This  has  been  done  by  Allea^263*  and  by  Kenner  and 
Allen. ^264^  For  the  change  in  thermal  expansion  relative  to  the 
bulk  value,  they  obtain  tbe  expression  for  tho  layer 
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whore  5  is  the  potential  energy  matrix,  Wp  is  the  p  normal  mode 
frequency,  Yj(u>p)  is  the  Grtineisen  parameter  given  by 

?!<»>  -  -  .  (332) 

ai  is  the  ith  lat'tice  spacing,  and  aQi  is  a  reference  latticp  spac- 
.  ing.  Kenner  and  Alien  have  carried  out  detailed  calculations  for 
the  (100)  and  (111)  surfaces  of  Ar,  Kr,  and  Xo  using  Lennard-Jones 
interactions.  They  find  a  peak  in  e  /»b  at  low  temperatures  assoc¬ 
iated  with  the  dispersion  of.  the  surface  modes.  Their  results  for 
'  Xe  are  significantly  lower  than  the  experimental  results  of 
Ignatjevs. 

An  alternative  formulation  of  surface  thermal  expansion  has 

recently  been  given  by  Dobrzynski  and  Maradudln. ^265^  They  evaluate 

the  Helmholtz  free  energy  including  enharmonic  terms  with  tho  aid 

of  many- body  theoretic  techniques.  Specific  calculations  of  tho 

surface  thermal  expansion  for  the  (100)  surface  of  or-iroa  were 

(32) 

made  using  the  model  of  Clark  et  al.  They  found  that  «s/»b 

is  somovhat  greater  than  unity  at  high  temperatures ,  but  increases 
rapidly  at  temperatures  below  5GPK. 
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dashed  lias  and  a  genera lifted  -surface  "save  u>  the  right 
of  this  lino, 

0.  Frequercy  versus  vavc:  vector  curves  Sor  a  soonatcaic  fee 
crystal  vitfe  two  (Ill)  surfaces.  The  surface  Staves  are 
denoted  by  S4.  <&}  5.1-layer  crystal;  (b)  Brillouir. 

A 

zono;  ic}  21- layer  crystal:  (<S>  21-i&yer  crystal  v?ith 
outeraost  layers  consisting  of  light  adsorbed  atoms, 
Adatom-host  atom  r&ss  ratio  1:5.  (After  Allen  et  ai„ , 
Ref.  43), 

10.  Frequencies  of  transverse  optical  surface  modes  versus 
wave  vector  Sot  &  (001)  surface  of  LiF.  (a)  high  fre¬ 
quency  sanies;  (b)  low  frequency  modes.  {After  Kliewer 
ami  Fuchs,  .  S§.) 

11.  Fruquenc i*o  versus  «av@  vector  for  a  15-layer  crystal 

of  H&C6  with  tvo  {001)  surfaces.  The  surface  modes  are 
labeled  by  S,  ana  missed  sodas  by  (After  Chen, 

*6S  at. ,  &a£.  87), 

12.  Intensity  versus  temperature  for  electrons  scattered 
from,  the  £110)  surface  of  Hi.  (a)  close  to  the. [0011 
dinectioa;  (b)  close  to  the  {1103  direction.  (Aftor 
Mae?.ae,  Ref.  146.) 

IS,  Moan  square  displacements  versus  atomic  layer  index  for 
the  (11.0}  surface  of  Hi.  a  is  the  nearest  neighbor  force 
constant. 

14.  *reffiperhiir6  dependence  of  the  surface- to- bulk  mean 
square  displacement  ratios  for  the  (100) ,  (110) ,  and 
•  (111)  surfaces  of  a  sana tonic  fee  crystal.  {After 
Alien  and  Ivetta .  Ref.  i?S), 
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Fig.  IS.  iienn  square  velocity  versus  atomic  position  for  a  mona¬ 
tomic  linear  chain  of  20  atoms  at  0 9X, 

Fig.  16.  Enorgy-loss  spectrum  for  a  particular  value  of  its.  The 
2 

lino  at  u>  “1  corresponds  to  the  surface  mode  and  tbo 
broad  band  to  tho  buiff  modoa.  (After  lioondy  and  slille. 
Ref.  200}. 

Fig.  17.  Iso-intensity  curvos  in  the  q, ,  plane  for  surface 

.  scattering  and  for  bulk  scattering.  In  each  case,  tho 
situations  Q3  «  Qj  and  Q.,  »  Q1(  are  presented. 

Fig.  18.  Extinction  cross-section  in  units  ot  the  geometrical 
cross-section  versus  frequency  for  a  small  spherical 
.  NaCi  crystal  of  radius  ljiss,  (After  Rupp  in  and  Englsan, 
Ref.  71). 

Fig,  19.  Absorption  coefficient  versus  frequency  for  a  crystal 
of  .laCt  100  layers  thick.  (After  Grits  ot  ai,.  Ref. 


to'jtJUct 


ATOMIC  LAYER  INDEX 


<u2>,(liO) 


MEAN  SQUARE  VELOCITY 


